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

УДК 629.7.064.3

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

© 2011 А. М. Гареев, А. Н. Коптев, Т. М. Гареев

Самарский государственный аэрокосмический университет имени академика С. П. Королёва (национальный исследовательский университет)

Рассмотрен метод моделирования современных гидравлических комплексов на основе применения теории динамических аналогий. Описан алгоритм моделирования состояния гидравлической системы по параметрам рабочей жидкости.

Гидравлическая система, рабочая жидкость, диагностика, упреждающее обслуживание, чистота рабочей жидкости, математическое описание.

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

Своевременно выполненный анализ позволяет определить начало интенсивного

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

Рис. 1. Электрогидравлический рулевой привод с механической упругой позиционной обратной связью: 1 - входное звено от механической системы управления; 2 - отсечной клапан; 3 - электрогидравлический сервозолотник; 4 - поршневой преобразующий гидроцилиндр; 5 - центрирование управляющего поршня; 6 и 10 - датчик линейных перемещений контура обратной связи; 7 - пружинный стопор автопилота; 8 - поршень отключения автопилота; 9 - главный управляющий золотник

ния состояний, в которых может находиться подобная система в период эксплуатации.

Цель данной работы - получение математического описания сложных автоматизированных технологических систем с учётом общеизвестных законов, справедливых для данной области физических явлений [1]. Эти законы отражают опыт, накопленный ранее для других аналогичных систем. Наиболее часто используются законы Кирхгофа (в электрических системах), принцип Даламбера (в механике), уравнения кинетики химических реакций, уравнения энергетического (в частности теплового) баланса, уравнения материального баланса, вытекающие из общих законов сохранения энергии и массы, и т. д.

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

Подтверждением правильности такого выбора являются результаты, полученные на основе метода динамических аналогий [1]. Этот метод использует аналогию сходных величин различных физических систем. Такие аналогии для величин электрической и гидравлической систем представлены в табл. 1.

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

сывающего исследуемую систему. Для этого мысленно расчленяя физическую систему в направлении передачи энергии на ряд элементов, которые соответствуют участкам системы, обладающим инерционностью, сопротивлением, емкостью и т. д., заменяем каждый элемент аналогами в соответствии с табл. 2. Множители р и 1/р в табл. 2 означают операторы дифференцирования и интегрирования, соответственно. Затем эти члены объединяются в уравнение в соответствии со структурной схемой анализируемого объекта.

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

Однако построить эквивалентную схему реальных ГС в терминах электрических

Таблица1

Электр ические системы Механические системы Гидравлические (пневматические) системы

Поступательное движение Вращательное движение

Напряжение, и Сила, Крутящий момент, М Давление, Р

Ток, 1 Линейная скорость, V Угловая скорость, N Расход, Ч

Заряд, Ч Линейное перемещение, X Угловое перемещение, ф Объем, V

Активное сопротивление, К Коэффициент вязкого трения, Б Сопротивление вращению, Кв Коэффициент гидравлического сопротивления, Кг

Индуктивность, ь Масса, т Момент инерции, I Инертность, J

Емкость, С Эластичность, 1/К Эластичность, 1/См Сжимаемость, сг

Таблица 2

Электрические системы Гидравлические системы

Сопротивление иш = Ш Шг Q

Индуктивность иь = рЫ р^О

ТТ 1 • Емкость ис =-1 с рС 1 О РСГ

1-й закон Кирхгофа £ /\ = 0 3 £ о, = о 3

2-й закон Кирхгофа £и , = и ] £ Рз = р 3

Потери энергии в сопротивлении Ш 2 Шг О2

И2 Запас энергии в индуктивности ^ ^2 2

си2 Запас энергии в емкости —— Сг Р2 2

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

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

В качестве примера рассмотрим участок ГС, схематически изображённый на рис. 2.

Рис. 2. Участок гидравлической системы

Вследствие трения в элементе 1 давление падает. По аналогии с электрическим сопротивлением (активным) падение давления АР1 = Яг0. Элемент 2 может аккумулировать энергию движущейся жидкости (аналог электрической емкости):

АЛ =-

2 C

Г" \ Qdt .

Соответственно для элемента 3 полу-

чим:

А Р,

I

ж

Для участка в целом модель изменения параметров рабочей жидкости будет иметь вид

Реых - Рвх = £ АР = IГ Ж-+С-1 Qdt + ШГ Q .

г=1 Ж С Г

Таким образом, сложный электрогидромеханический комплекс аналогично рассмотренному примеру можно описать Л-сетью и исследовать средствами тензорного анализа цепей [2]. При этом, если число элементов систем настолько велико, что расчёт её как целого производить слишком трудоёмко или невозможно, то необходимо применять метод расчёта по частям, который получил название диакоптика. Модель сис-

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

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

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

N =Ьрй! 612,

где N - мощность, кВт; Ар. - перепад давления, кгс/см2; Qi - расход, л/мин.

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

Расходными характеристиками называют зависимости между двумя гидравлическими параметрами Ар=реьх - рех=/^). Тогда импеданс /-го элемента гидравлической подсистемы может быть определён как производная от Ар по Q :

^ =

в литературе [4]. Например, для фильтров, использующих фильтроматериалы со сложной поровой структурой, используют расход -ную характеристику

ьр = т Q

(1)

Расходные характеристики элементов гидросистем достаточно подробно описаны

кБ „

где к - удельная пропускная способность единицы поверхности фильтроматериалов при перепаде давления в 0,1МПа; Бф - площадь фильтроэлемента; ц - динамический коэффициент вязкости. Соответственно, импеданс фильтра

2Ф = ц/(кБ ф ).

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

Ьр л = к а, Ьрт = кЛ2,

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

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

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

1. Закон узлов. Сумма всех «продольных переменных» в любом узле системы должна быть равна нулю.

2. Закон контуров. Сумма всех «поперечных переменных» вдоль любого замкнутого контура системы равна нулю.

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

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

п

Рн =ТАР], (2)

1=1

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

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

При любых преобразованиях структуры цепи количество переменных (размер-

ность пространства цепи) остаётся постоянным, равным количеству ветвей.

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

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

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

На основе контурного способа описания гидравлической подсистемы с учётом (1) уравнение, описывающее преобразование потока отдельной ветвью ГС, можно представить как Ар = г • Q .

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

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

Р = 1 • Q

(3)

где Р - компаунд-тензор значений давления источников питания, по контурам гидравлической подсистемы; Q - компаунд-тензор значений расхода жидкости на элементах

контуров; Z - компаунд-тензор импедансов элементов. Компаунд-тензор импедансов элементов гидравлической подсистемы будет составляться на основе матрицы M, описывающей объект - конкретную гидравлическую подсистему с учётом расходных характеристик каждого её компонента.

На основе тензорного выражения (3) можно составить эквивалентную модель ГС, в которой каждый элемент гидравлической подсистемы будет заменен эквивалентным элементом электрических цепей с импедансом, описываемым матрицей Z. Для расчёта этой эквивалентной модели можно использовать аппарат тензорного анализа электрических сетей (цепей) [2].

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

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

Данная статья написана по результатам проведения поисковой научно-исследовательской работы в рамках реализации ФЦП «Научные и научно-педагогические кадры инновационной России» на 2009 - 2013 годы.

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

1. Кулик, В. Т. Алгоритмизация объектов управления: справочник [Текст] / В. Т. Кулик. - Киев: Наукова думка, 1968. - 363 с.

2. Петров А. Е. Тензорная методология в теории систем [Текст]/ А. Е. Петров. -М.: Радио и связь,1985. - 152 с.

3. Гареев, А. М. Упреждающее обслуживание гидравлических систем летательных аппаратов [Текст]/ А. М. Гареев, С. Н. Тиц. -Самара: Изд-во СНЦ РАН, 2010. - 112 с.

4. Шумилов, И. С. Системы управления рулями самолётов [Текст]/ И. С. Шумилов. - М.: Изд-во МГТУ им. Н.Э. Баумана, 2009. - 469 с.

USING THE THEORY OF DYNAMIC ANALOGIES FOR MODELLING AUTOMATED MANUFACTURING SYSTEMS

© 2011 A. M. Gareyev, A. N. Koptev, T. M. Gareyev

Samara State Aerospace University named after academician S. P. Korolyov (National Research University)

The paper presents a method of modeling up-to-date hydraulic complexes on the basis of using the theory of dynamic analogies. An algorithm of modeling the state of a hydraulic system on the basis of working fluid parameters is described.

Hydraulic system, working fluid, diagnostics, proactive maintenance, fluid contamination, mathematical description.

Информация об авторах Гареев Альберт Минеасхатович, кандидат технических наук, доцент кафедры эксплуатации авиационной техники, Самарский государственный аэрокосмический университет имени академика С. П. Королёва (национальный исследовательский университет). E-mail: [email protected]. Область научных интересов: оперативный контроль технического состояния авиационной техники; неразрушающие методы контроля.

Коптев Анатолий Никитович, доктор технических наук, профессор, заведующий кафедрой эксплуатации авиационной техники, Самарский государственный аэрокосмический университет имени академика С. П. Королёва (национальный исследовательский университет). Область научных интересов: оперативный контроль технического состояния авиационной техники; неразрушающие методы контроля.

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

Gareev Albert Mineaskhatovitch, Ph.D., Associate Professor, department of aircraft system maintenance, Samara State Aerospace University named after academician S. P. Korolyov (National Research University). E-mail: [email protected]. Area of research: fast control of aircraft system technical condition, non-destructive methods of control.

Koptev Anatoliy Nikitovich, Ph.D., Professor, Head of aircraft system maintenance department, Samara State Aerospace University named after academician S. P. Korolyov (National Research University). Area of research: expeditious control of aircraft system technical condition, non-destructive methods of control.

Gareev Talgat Mineaskhatovitch, post-graduate student, department of aircraft system maintenance, Samara State Aerospace University named after academician S. P. Korolyov (National Research University). Area of research: fast control of aircraft system technical condition, nondestructive methods of control.