Научная статья на тему 'Математическое обеспечение метода контроля технического состояния гидросистем воздушных судов по параметрам рабочей жидкости'

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

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

Аннотация научной статьи по математике, автор научной работы — Ковалев М. А.

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

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

Похожие темы научных работ по математике , автор научной работы — Ковалев М. А.

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

SOFTWARE OF THE QUALITY MONITORING OF THE TECHNICAL CONDITION OF HYDROSYSTEMS OF AIR COURTS ON PARAMETERS OF THE WORKING LIQUID

Article is devoted to development of theoretical base of a method of definition of a technical condition of hydraulic systems on the basis of the analysis of parameters of a working liquid. The generalized statistical model of the control of serviceability of hydrosystems on parameters of a working liquid is resulted. The mathematical model of a condition of a working liquid of the hydrosystems, describing dynamics of change of parameters of a liquid and their mutual influence is received. Both models are submitted in vector and tensor kinds. The methodology of construction and algorithm of work of the monitoring system of a technical condition of hydrosystems on parameters of a working liquid is specified.

Текст научной работы на тему «Математическое обеспечение метода контроля технического состояния гидросистем воздушных судов по параметрам рабочей жидкости»

УДК 62.001.4

МАТЕМАТИЧЕСКОЕ ОБЕСПЕЧЕНИЕ МЕТОДА КОНТРОЛЯ ТЕХНИЧЕСКОГО СОСТОЯНИЯ ГИДРОСИСТЕМ ВОЗДУШНЫХ СУДОВ ПО ПАРАМЕТРАМ РАБОЧЕЙ ЖИДКОСТИ

© 2008 М.А. Ковалев Самарский государственный аэрокосмический университет

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

В последнее время большое внимание специалистов уделяется проблеме повышения надежности бортовых систем воздушных судов (ВС). Одной из наиболее важных из них, техническое состояние которой непосредственно влияет на безопасность полетов, является гидравлическая система (ГС).

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

При определения технического состояния ГС по параметрам РЖ используется до-пусковый метод контроля. Этот подход предполагает оценивание диагностических признаков и сопоставление полученных оценок с допустимыми значениями.

Допустим, что для определения технического состояния ГС необходимо произвести оценивание т признаков, характеризующих состояние РЖ, которые образуют вектор оцениваемых величин Х=[ х, х2 ..., хт]Т. В общем случае величины х. являются взаимо-

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

Геометрически состояние РЖ в фиксированный момент времени I можно охарактеризовать некоторой точкой в т - мерном пространстве, соответствующей оконечности вектора Х(1). Эта точка в процессе эксплуатации занимает различные положения, которые в совокупности представляют область £ возможных состояний РЖ. Причем метод контроля технического состояния ГС по параметрам рабочей жидкости предполагает, что работоспособному состоянию ГС будут соответствовать определенные значения признаков х То есть на каждый признак задаются допуски х'. и (или) х ". , характеризующие минимальное и (или) максимальное значения, соответствующие работоспособному состоянию ГС. Тогда работоспособному состоянию ГС будут соответствовать определенная область 5р т - мерного пространства значений оцениваемых величин РЖ ( 5 р е 5 ), ограниченная поверхностью, положение которой задается значениями допусков на эти величины. До тех пор пока рабочая точка находится в пределах этой ограниченной области, ГС может считаться работоспособной.

Иными словами, геометрически в т -мерном пространстве область возможных

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

Х(г) - а = 0, (1)

где а - вектор, описывающий поверхность рабочей области S.

Уравнение (1) является обобщенной статистической моделью контроля работоспособности ГС на основе анализа параметров РЖ. При построении модели для конкретной системы контроля наиболее важной и сложной является задача оптимизации выбора признаков, составляющих вектор X и подлежащих оцениванию. Предпочтение следует отдавать тем признакам, которые имеют наибольшую диагностическую ценность, позволяют обнаружить неисправности на возможно более ранних стадиях их развития и которые удобно определять, измерять в процессе эксплуатации и обрабатывать при анализе.

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

Очевидно, что наименование и количество анализируемых параметров будут определяться спектром задач, решаемых системой контроля. Так способность ГС выполнять свою основную функцию (передача гидро-

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

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

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

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

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

Для исследования динамики изменения параметров РЖ и их взаимного влияния целесообразно построить математическую модель, описывающую состояние РЖ ГС.

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

У, ^) = /г У^ У2,•••, Уг-^ У¡+1,■■■, Уп ). (2)

Соответственно вектор параметров Y(t) можно описать выражением

Y ^) =

" У1(^ )" /1(t, У 2 , Уз^^ Уп )

У2(t) /2(t, Уl, Уз^^ Уп )

у, (t) = (t, Уl,■, У, Уг'+1, ^ ", Уп )

_уп (t)_ _ Л (t, Уl, У2 , • " ", Уп-1)

или Y(t)=F(t,Y).

(3)

Для того чтобы учесть динамику изменения параметров и их взаимное влияние воспользуемся представлением параметров у (^ путем их разложения в ряд Тейлора.

Пусть вектор Y(t) содержит лишь один параметр у(0. Тогда значение этого параметра в момент времени t+Дt можно определить на основе выражения

у(г + ^) = у(г) + Ш At +1 Дt 2 +

Л 2 ёг

1 Л3у(0 . 3

+--¿±zAt3 + ■

6 ё?

(4)

В случае, когда анализу подлежат п параметров РЖ динамика изменения параметра убудет описываться уравнением [3]

1 2 1 3

у,. ^ + At) = у,. (t) + й/г + - ё2 / + - ё3/ + •, (5)

2 6

где /г = /г Уl, у 2 , Уг-1, Уг+1,■, уп ),

ёк/ =

д

— At + Е

я* ¿—1

дt

1=1

Х\

Ау 1

/г (6)

- условное обозначение полного дифференциала к - го порядка функции /. Его надо понимать так: сначала "возводим в к - тую

степень" многочлен д^ ^ + Е

удУ]

Ау,

так

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

/ = д/- At + Е К Ау.

ё 2/ =

а2

^2

1=1

1=1 ыу. 'д2/

чду;

Ay

+ 2 ^ At х

дt

хЕ

1=1

кдУ]

л

(

Ay.

+ 2Е

1=1

д/

(

д/

Л

При этом

ёу,

= 0

д AУj £ Ыу- Ayl

дУ} I=1 ^

I * 1

в силу независимости

/1(1у1,у2,...,у1_1,у1+1,...у) оту,.

С учетом (3) уравнение (5) можно записать в векторном виде

к

1 2

У (Х + ДХ) = У (Х) + dF (Х, У ) + - d2 F(Г, 7 ) + 1

+ - d3 F (Х, У ) + ... 6

где

(7)

—кГ = [—к^ —kf2 ... —7„ ]г. (8)

Уравнения (6) - (8) описывают процесс изменения во времени вектора параметров РЖ, а также учитывают взаимное влияние параметров. Однако в правой части выражения (7) присутствует бесконечное число слагаемых. Поэтому следующим шагом при построении математической модели состояния РЖ является определение количества слагаемых для каждого из параметров, образующих вектор У, обеспечивающих решение поставленных задач с заданной точностью. С учетом осуществления этого шага уравнения (6) - (8) можно назвать математической моделью состояния РЖ ГС. На основе ее анализа можно определить величины, характеризующие состояние РЖ и подлежащие оцениванию в ходе контроля технического состояния ГС, т.е. построить вектор оцениваемых признаков X(Х). Источниками информации о значениях этих признаков являются датчики. Кроме того, вектор Х(Х) может содержать некоторые бинарные величины, характеризующие наличие или отсутствие сигналов от различного рода сигнализаторов (реле), установленных в ГС. Например, на большинстве самолетов в ГС устанавливаются сигнализаторы давления, срабатывающие при снижении давления в ГС ниже заданного значения, что сигнализирует об отказе ГС.

Таким образом, метод контроля технического состояния ГС по параметрам РЖ с использованием модели (6) - (8) предполагает следующую методику построения и функционирования системы контроля:

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

- определяем наименование и количество параметров РЖ, подлежащих оцениванию (формируем вектор параметров РЖ У);

- строим математическую модель состояния РЖ на основе уравнений (6) - (8), определив при этом число слагаемых в правой части выражений, описывающих каждый из параметров вектора У;

- формируем вектор оцениваемых величин X и с учетом допустимых значений диагностических признаков на основе обобщенной модели (1) синтезируем статистическую модель контроля работоспособности анализируемой ГС;

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

Для наглядности построим статистическую модель контроля работоспособности одноконтурной ГС, в которой оценивается только давление РЖ в линии нагнетания р(Х). Будем считать, что в данной ГС помимо датчика давления РЖ установлен также сигнализатор давления. Такая система контроля работоспособности ГС используется на многих отечественных самолетах. В этом случае математическая модель состояния РЖ ГС будет описываться уравнением (4), в котором

у(х)=р(х).

Исходя из условий конкретной задачи и выбрав интервал ДХ<</, количество слагаемых в правой части выражения (4) можно существенно сократить, что значительно упростит данное уравнение. Так, например, если ограничиться тремя слагаемыми, то математическая модель изменения давления примет вид приближенного равенства р (Х + Д Х) « р (Х) + • А Х + у —'Д12. (9)

Анализ выражения (15) показывает, что значение давления РЖ ГС в момент времени Х+ДХ будет определяться значением этого параметра в момент времени Х , величиной и—решала Д=Х, а та—рр(значениями скорости -—

и

ускорения

—х2

изменения давления.

С учетом выражения (9) статистическая модель контроля работоспособности ГС будет описываться следующей системой уравнений

р» < р(х) < рв; , —р(Х)

р» <^тт<рв;

—Х

и —2р(Х)

р» < Л ;

рс = 0.

(10)

Г> ))>>>>

Здесь р р р», рв, рн и рв - нижние и верхние допустимые значения (границы) дав-

ления, а также скорости и ускорения изменения давление соответственно. Последнее уравнение в системе (10) характеризует ситуацию, когда сигнал на выходе сигнализатора давления отсутствует. При срабатывании сигнализатора р==1.

Модель (10) предполагает, что объект -ГС - работоспособен, если выполняются все условия, составляющие систему Если же хотя бы одно из условий не выполняется, то объект неработоспособен.

Уравнения (6) - (8) в векторной форме записи удобны для представления только несложной одноконтурной ГС. Вместе с тем ГС ВС являются многоконтурными и состоят из нескольких отдельных ГС, каждая из которых обслуживает своих потребителей. Например, на тяжелых самолетах используют три или четыре самостоятельные, независимые друг от друга гидросистемы. Это позволяет повысить суммарную мощность бортовой гидросистемы, а также увеличить ее надежность за счет многократного резервирования функциональных подсистем [2]. Причем используются различные схемы резервирования. Они предполагают либо одновременное подключение функциональных подсистем к нескольким отдельным ГС ("горячее" резервирование), либо переключение с отказавшей отдельной ГС на исправную ("холодное" резервирование). Тем самым отдельные ГС, составляющие бортовую ГС, действуя независимо, взаимно дополняют друг друга.

Многоконтурное построение ГС и сложные схемы резервирования значительно усложняют задачу анализа функционирования ГС. Для решения этой задачи целесообразно использовать представление ГС на основе математического аппарата тензорного анализа сетей [4]. Такое представление позволяет описать с необходимым уровнем детализации многоконтурные ГС и учесть резервирование любого типа и уровня сложности. Соответственно при таком описании ГС и математическая модель состояния РЖ (6) - (8) также должна быть представлена в тензорном виде. Для решения этой задачи воспользуемся индексами и обозначениями, применяемыми в тензорном анализе [4]. Тогда математическая модель состояния РЖ ГС будет иметь следующий вид

Ау,. '-Аук'+

(11)

у'(t + Аt) = у'(t) + -Ау/+И)

+ О. ■ Ау.'• Аук'■ Ау,'+ ...,

где I = 1...п; .,к,1 = 0...п; у - тензор валентности один (вектор) параметров РЖ; у = у1 ... уп ]Т - расширенный тензор (вектор) параметров, в котором роль нулевого члена выполняет время I; Ау = [А* Ау ... Ауп]Т - расширенный тензор (вектор) приращений параметров; 3 М, D - тензоры.

Тензор валентности два (матрица)

з; =

ду ду,

/ = 1...п; . = 0...п, (12)

содержит п(п+1) элементов. Для параметра с номером а (столбец ) он описывается вектором

V дуа дуа

за =

д ду

ОУп

Тензор валентности три (геометрический эквивалент - параллелепипед)

м;к =

1 д2 у

2 су,'дук

I = 1...п; к = 0...п, (13)

содержит п матриц с числом элементов (п+1)2. Тензор И. для параметра с номером а - это матрица вида

И. =1 t 2

а д у 2а д у 2а д у

дгдг д^ Здуп

2а д у 2а д у 2а д у

дуд ЗУдУ ^Су,

д2 уа д2 уа д2 уа

дУпдУ: суп дУ,

о;

Тензор валентности четыре

_ 1 д3 у'

г, / = 1...п; к, I = 0...п,(14)

6 дУ.'дУк'дУ,

содержит п(п+1)3 элементов.

Таким образом, уравнения (11) - (14) описывают математическую модель состояния РЖ ГС в тензорном виде. В рассмотренном варианте записи полученные уравнения также как и выражения (6) - (8) можно применить для описания параметров РЖ несложной одноконтурной ГС. Математическая мо-

дель состояния РЖ сложной ГС, состоящей из м> отдельных ГС, должна описывать п параметров каждой из м> систем. Получить такую модель можно, если в выражении (11) каждый параметр уг в тензореуг, г=1...п, заменить тензором (вектором), описывающим значение этого параметра во всех м? отдельных ГС. Такой тензор, элементами которого являются также тензоры, называют мульти-тензором [4]. Обозначим его как у'И. Здесь 1=1...п - номер параметра, И = 1...ж - номер отдельной ГС в составе сложной ГС.

Если допустить, что параметры РЖ одной из ГС не влияют на параметры РЖ другой ГС, то мультитензорное представление ГС не приведет к существенному усложнению математической модели. Тогда в правой части выражения (11) вектор (тензор) Ay' необходимо заменить матрицей (мультитензо-ром валентности два) AyJ', столбцами которой являются векторы Ay' , составленные для каждой из отдельных ГС, а ^ М и D - муль-титензорами третьей, четвертой и пятой валентности соответственно, т.е. выражение (11) примет следующий вид:

угИ ^ + At) = угИ ^) + Sfh ■ Ayhj '+МгкИ х

х AyJ 'Лу] '+Dfklh ■ AyhJ 'Лу] '^у^ '+ • (15)

При этом в случае анализа отдельной ГС с номером а необходимо в уравнении (15) переменную И заменить на постоянную а, что можно рассматривать как понижение размерности всех тензоров в этом уравнении. Тогда валентность величин в уравнении (15) будет равна размерности аналогичных величин выражения (11), а вычислительная сложность задачи при переходе к мультитензорам не изменится. Фактически для того, чтобы описать состояние РЖ сложной ГС необходимо решить задачу (11) - (14) для каждой из отдельных ГС.

Попытка учесть взаимную зависимость параметров различных ГС приведет к существенному усложнению модели (11) - (14), поскольку при этом значительно возрастет валентность тензоров в правой части уравнения (11), а, следовательно, и вычислительная сложность задачи.

Следует отметить, что в реальных ГС параметры РЖ отдельных ГС можно считать независимыми. Поэтому выражение (15) с достаточной степенью достоверности описывает состояние РЖ многоконтурной ГС.

С учетом выражения (1) обобщенную статистическую модель контроля работоспособности многоконтурной ГС по параметрам РЖ можно описать следующим уравнением хгИ - а = 0, (16)

где г=1...т, И=1...ж; хгИ и а - мультитензоры одинаковой размерности и валентности.

Таким образом, статистическая модель (16) позволяет описать процесс контроля работоспособности сложных ГС по параметрам РЖ. Эта модель предполагает следующую методику построения системы контроля.

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

2. Для выбранных параметров РЖ синтезируется математическая модель, описывающая состояния РЖ ГС. Для этого используются выражения в векторной (6) - (8) или тензорной (11) - (15) формах записи, в которых исходя из требуемой точности решения задачи определяется количество слагаемых в правой части каждого из уравнений.

3. На основе математической модели состояния РЖ формируется вектор (тензор) величин, подлежащих оцениванию в ходе определения работоспособности ГС. Этот вектор (тензор) будет содержать параметры РЖ и производные от них величины, характеризующие динамику изменения параметров и их взаимное влияние.

4. Используя выражения (1) в векторной или (16) в тензорной формах записи строится статистическая модель контроля работоспособности анализируемой ГС по параметрам РЖ, в основе которой лежит допусковый принцип. Она определяет алгоритм работы и структуру системы контроля.

СПИСОК ЛИТЕРАТУРЫ

1. Башта Т.М., Бабанская В.Д., Головко Ю. С. и др. Надежность гидравлических систем воздушных судов / Под ред. Т.М.Башты. М.: Транспорт, 1986.

2. Акопов М.Г., Бекасов В.И., Долгушев В.Г. и др. Системы оборудования летательных

аппаратов: / Под ред. Матвеенко А.М., Бекасова В.И. М.: Машиностроение, 2005.

3. Выгодский М.Я. Справочник по высшей

математике. М.: Астрель; АСТ, 2002.

4. Крон Г. Тензорный анализ сетей / Под редакцией Л.Т.Кузина, П.Г.Кузнецова. М.: Советское радио, 1978.

SOFTWARE OF THE QUALITY MONITORING OF THE TECHNICAL CONDITION OF HYDROSYSTEMS OF AIR COURTS ON PARAMETERS OF THE WORKING LIQUID

© 2008 M.A. Kovalev Samara State Aerospace University

Article is devoted to development of theoretical base of a method of definition of a technical condition of hydraulic systems on the basis of the analysis of parameters of a working liquid. The generalized statistical model of the control of serviceability of hydrosystems on parameters of a working liquid is resulted. The mathematical model of a condition of a working liquid of the hydrosystems, describing dynamics of change of parameters of a liquid and their mutual influence is received. Both models are submitted in vector and tensor kinds. The methodology of construction and algorithm of work of the monitoring system of a technical condition of hydrosystems on parameters of a working liquid is specified.

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