УДК 681.586.35
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ И АЛГОРИТМ РАСЧЁТА ЧАСТОТНЫХ ХАРАКТЕРИСТИК ГАЗОВЫХ ИЗМЕРИТЕЛЬНЫХ ЦЕПЕЙ
ДИФФЕРЕНЦИАЛЬНОЙ СХЕМЫ
© 2012 А. А. Игонин, А. Г. Гимадиев
Самарский государственный аэрокосмический университет имени академика С. П. Королёва (национальный исследовательский университет)
Описаны алгоритм и программа расчёта частотных характеристик дифференциальной газовой измерительной цепи на основе разработанной её математической модели. Приведены основные формулы математической модели, дано описание алгоритма функционирования разработанной программы и основные правила работы с ней.
Газовая измерительная цепь дифференциальной схемы, тонкие трубки, ёмкость, мембрана, дросселирующий элемент, математическая модель, алгоритм и программа расчёта, частотная характеристика.
Важной проблемой, возникающей в процессе эксплуатации газовых измерительных цепей дифференциальной схемы, является обеспечение требуемой статической и динамической точности, которая зависит как от погрешности первичных преобразователей, так и от погрешности, возникающей при передаче давления от точки отбора в объекте контроля до преобразователя перепада давления. Необходимость применения газовых измерительных цепей в технологических и энергетических установках обусловлена многими причинами, основными из которых являются условия эксплуатации: высокий уровень температур в объекте контроля, ограниченность пространства под установку первичных преобразователей, жёсткие требования по точности измерения малых значений перепадов давлений. Например, при доводке газотурбинного двигателя возникает необходимость в измерении пульсаций воздушного скоростного потока на фоне изменяющегося статического давления [1]. Учитывая небольшие значения давления заторможенного потока воздуха относительно статического давления, можно рекомендовать для измерения параметров газового потока только измерительную цепь дифференциальной схемы. Однако в такой цепи могут возникать ре-
зонансные явления, вносящие существенную динамическую погрешность в результат измерения. Поэтому возникает необходимость во включении в измерительную цепь акустических корректирующих элементов, способствующих повышению точности измерения пульсаций скоростного напора.
Другим примером является то, что в магистральном газопроводе с измерительной системой на основе диафрагменных расходомеров необходим учёт влияния колебаний давления на точность измерения расхода газа [2]. Такой учёт возможен, если одновременно измерять стационарную составляющую расхода газа и пульсации перепада давления на мерной диафрагме. Применение датчиков динамического перепада давления на мерной диафрагме требует дополнительных материальных затрат и времени. Поэтому для компенсации пульсационной погрешности газовой измерительной цепи расходомера предлагается воспользоваться данными по колебаниям давления, измеренными датчиками давления и перепада давления собственно самой системы. Однако при передаче акустического сигнала - пульсаций давления - от точки отбора до преобразователя происходит его искажение, которое необходимо восстанавливать.
В статье рассматривается решение задачи по разработке алгоритма и программы для расчёта частотных характеристик (ЧХ) газовой измерительной цепи дифференциальной схемы, подключённой к двум точкам объекта контроля, в которых необходимо измерять как стационарные, так и пульсационные составляющие давления. На основе анализа результатов расчёта ЧХ можно рекомендовать параметры измерительной цепи, при которых удовлетворяются заданные требования к динамической точности измерительной цепи.
В математической модели газовой измерительной цепи приняты следующие допущения:
- пульсации давления и перепада давления представлены в виде полигар-монических сигналов;
- объект контроля со стороны обеих ветвей дифференциальной измерительной цепи описывается независимыми источниками колебаний давления, которые при необходимости анализа могут генерировать одинаковый спектр гармонических колебаний;
- каждый элемент цепи рассматривается как акустический четырёхполюсник;
- элементы измерительной цепи -линейные;
- движение газа в импульсных линиях рассчитывается по ламинарной модели с учётом нестационарности;
- разделительная мембрана преобразователя перепада давления рассматривается как упругость с учётом инерционных свойств и сил трения.
Основными элементами газовой измерительной цепи, для которых разработана её математическая модель с учётом принятых допущений, являются: трубка малого диаметра, ёмкость, мембрана, дросселирующий элемент [3,4]. В математической модели измерительная цепь представлена рядом соединённых каскадно четырёхполюсников, к концам которой подключены источники полигармониче-ских пульсаций давления [5]. Расчёт час-
тотных характеристик газовой измерительной цепи проводится импедансным методом.
Алгоритм расчёта частотных характеристик измерительной цепи построен на методе суперпозиции, основанном на раздельном расчёте колебаний давления и расхода газа в выделенных сечениях от каждого из источников с дальнейшим их суммированием. При этом последовательная измерительная цепь разбивается на ряд элементов с чётко обозначенными свойствами. Вначале расчёт ведётся при действии одного из источников, а другой источник при этом считается пассивным. Каждому элементу присваивается определённый номер, начиная от одного из источников колебаний давления. В алгоритме за первый элемент цепи принят тот, который присоединён к пассивному (второму) источнику.
Колебания давления в активном (первом) источнике колебаний давления описываются зависимостью:
Ссо=С+ірГє'4', (і)
к=1
где Р^0 - постоянная составляющая давления; (дк, р(1;к) - частота и комплексная
амплитуда к -й составляющей спектра колебаний давления активного (первого) источника; К - число учитываемых гармоник в спектре колебаний первого источника;] - мнимая единица; і - время.
После предварительного вычисления коэффициентов четырёхполюсников всех выделенных элементов цепи [3,4], начиная от первого элемента, подключённого к пассивному источнику с внутренним импедансом 7(2к), рассчитывают входные импедансы элементов для к -й составляющей спектра колебаний давления по формуле [3]:
л(\,к) 7(2,к) + п(1,к)
7(1,к) _ А ^ П1 (2)
£,(1,к) 7(2,к) + £)(1,к) ’ У '
где А; , п?’к) - коэффициен-
ты матрицы передачи первого относительно пассивного (второго) источника колебаний элемента цепи.
Далее рассчитывают последовательно входные импедансы каждого из І - ых выделенных элементов цепи вплоть до входа первого элемента цепи (І = I), подключённого к активному источнику колебаний давления по формуле:
л(\,к) 7(1,к) + о(1,к)
7(1,к) _ ^1 ^уй,(!-1) I
"Vк1 с (1,к) 7 (1,к) + П(1,к)
I 'Чк^-1) ^1
(3)
где ZV;k)(I-1) - входной импеданс предыдущего
дол) в('1’к) с^,к) о(~1,к)
р(1,к) = р(1,к) 7(1,к) /(7(1,к) + 7(1,к)); ^JI ,1 УЫ ук,1 \ Іві ук,1 /з
q(1.к) = р(1.к) /( 7(і.к) і 7 (1.к))
ЧI1 г ізі \ ізі ук1 )■
(4)
На входе в чувствительный элемент дифференциального датчика, соответствующего номеру г = у, комплексные амплитуды колебаний давления и расхода определяются по зависимостям:
р(1,к) _ р(1,к) А(\,к) (1,к) 0(1,к).
р],2 р]+1,1 Л]+1 Чі+1,1 В]+1 ;
_(1,к) _ -С(1,к) р(1,к) + п(1,к) п(1,к) Ч],2 +1 /О+и ^]+1Ч] +1,1-
(8)
Комплексная амплитуда колебаний давления на выходе чувствительного элемента датчика, соответствующего сечению (у -1,1), определяется с учётом ра-
венства q
X) _ q^12k) по формуле:
элемента цепи;
коэффициенты
матрицы передачи I - го, конечного элемента цепи для к -й составляющей спектра.
При известном входном импедансе цепи ) со стороны активного источника колебаний определяют комплексные амплитуды колебаний давления и расхода на входе в цепь по формулам:
(1,к) _ 7(1,к) Г1,к)
Л',1 ук, І-1Ч ],1 •
(9)
Аналогично описанной процедуре определяются комплексные амплитуды колебаний давления от первого источника на входе в чувствительный элемент дифференциального датчика и на его выходе для всех гармоник (к=1,..., К).
Расчёт комплексных амплитуд колебаний давления и расхода в измерительной цепи при активации второго источника осуществляется аналогичным образом. При этом источник колебаний описывается зависимостью:
Комплексные амплитуды колебаний давления и расхода по концам элементов измерительной цепи рассчитываются по зависимостям: для I - го элемента
р(1,к) _ р(1,к) А(\,к) (1,к) в
PI2 PI1 А qi1 В
(1,к)
(1,к) _ 2
с?-
к) р(1,к) FI1
+п
(1,к) q (1,к) , ІІ1 з
(5)
для стыка I - го и I -1 - го элементов
р(1,к) _ р(1,к) г (1,к)
iJI-1,1 iJI ,2 з У! -1,1
qI(12k);
для I -1 - го элемента
р(1,к) _ р(1,к) А(1,к) - q(1,к) о(1,к) 1JI-1,2 1JI-1,1 ^-1 ‘/I-1,1^7-1 з
^’к) _ ^ -1,2
С(1,к) р(1.к) + П(1,к) q(1,к) '“'I-1 УI-1,1 т LJI-1 ЧI-1 •
(6)
п= N
Р(2)(і) = Р^2) + V р(2’и)е]дп‘
1 ізі V) 1ізі,0~/^ У ізі е
(10)
П=1
где рЫ20 - постоянная составляющая давления; дп, Р(2'")
частота и комплексная
амплитуда и -й составляющей спектра колебаний давления активного (уже второго) источника; N - число учитываемых гармоник в спектре колебаний второго источника.
При сохранении нумерации элементов цепи входной импеданс I - го или первого относительно первого уже пассивного источника элемента определится по зависимости:
7(2.n)
■'vh.1
f(2.n) z(1’n)
V ^ist
B
(2.n )
C(2-n) z1-1’" ) + D(2.n)
((()
где - внутренний импеданс первого
источника колебаний;
А(12,пБ{2,пС]-2’”В{2,п} - коэффициенты матрицы передачи I - го элемента цепи для п -й гармоники спектра колебаний давления от второго источника.
С учётом (4)...(7) получим соответствующие зависимости при действии только второго источника колебаний:
p
(2.n) = j .1
q(2-n) =
4j .1
p(2.n) j(2,n) q(2.n) d(2,") . pj-(.і Aj-l - qj-(.і Bj-l ;
-C(2>n) p(2.n) + D(2>n)q(2.n) ^ j-і Fj-1.1 ^ ^j-і 4j-l.l з
((2)
а с учётом выражения ) = д(1,п), соответствующего сечению (; +1,1):
„(2,п) = -(2,”) ^(2,п) (13)
-Г;,2 ^уИ,]+1Ч],1 • Ч1*-*/
Если принять во внимание линейность входящих в состав измерительной цепи элементов, то в результате суммирования составляющих давления можно получить соответствующие величины на входах чувствительного элемента дифференциального датчика со стороны первого и второго источников колебаний давления:
к=K
n=N
P?f’(t) = P;0 + X p £) e* + x p;.2," ) e
к=1
к=K
n=1
(14)
P‘‘r’(t)=PS+X p j!;k) e* + X p j?) e
jwnt
к=1
((5)
Перепад давления на чувствительном элементе датчика будет равен
к=K
,(1.k)
Dp°-2)(t)=PS - ps+x(pj.■
™n k= (16)
- p « + X (p - p j2rV'v
Изложенная методика позволяет рассчитывать комплексные амплитуды колебаний давления на входе и выходе дифференциального датчика колебаний давления и. при необходимости. в произвольном сечении газовой измерительной цепи. По результатам расчёта можно определить. с каким усилением или ослаблением передаются гармоники спектра колебаний давления от точек измерения в объекте контроля к датчику перепада давления. Если задать только одну гармонику колебаний давления с постоянной амплитудой и переменной частотой от какого-либо источника. можно построить амплитудно-частотную (АЧХ) и фазовочастотную (ФЧХ) характеристики измерительной цепи и оценить её передающие свойства.
По изложенной методике разработана программа расчёта частотных характеристик газовых измерительных цепей с использованием языка программирования C++ и системы управления базами данных Firebird [5].
Разработанное программное обеспечение позволяет рассчитывать частотные характеристики газовых измерительных цепей с определённой пользователем конфигурацией и задаваемым набором параметров. Программа имеет возможность расчёта серии частотных характеристик для измерительной цепи с различными значениями одного или нескольких параметров. В программе можно задавать состав и конфигурацию газовой цепи (количество элементов может быть произвольным). значения параметров каждого элемента цепи. а также указывать на изменяющиеся параметры и задавать для них диапазон и шаг их изменения. С помощью запросов к системе управления базами можно осуществлять выборку данных для построения любой рассчитанной частотной характеристики или серии частотных характеристик.
Обобщённый алгоритм работы программного обеспечения можно представить в следующем виде:
1. Ввод пользователем конфигурации газовой дифференциальной измерительной цепи и параметров каждого её элемента, в том числе и сечений, в которых значения давления и расхода важны для пользователя.
2. Ввод пользователем параметров расчёта, внешних условий и рабочей среды.
3. Задание изменяющихся параметров, пределов и шага их изменения.
4. Расчёт и сохранение результатов в базу данных.
5. Формулирование запроса к системе управления базой данных и построение серии частотных характеристик.
После задания всех необходимых параметров происходит вызов подпрограммы, которая рассчитывает частотные характеристики газовой измерительной цепи по алгоритму, обобщённая блок-схема которого представлена на рис. 1. Входными данными подпрограммы являются:
- список изменяемых параметров, содержащий информацию о названии параметра, указателе на его расположение внутри модели газовой измерительной цепи, текущее значение параметра, минимальное и максимальные значения параметра, а также шаг его изменения;
- описание модели газовой цепи с заданными начальными значениями параметров и конфигурацией цепи.
Действие программы заключается в циклическом расчёте частотных характеристик заданного вида для всех комбинаций параметров и сохранении их в базе данных.
Программа расчёта характеристик измерительной цепи после запуска выдаёт на экран окно для задания параметров газовой цепи и расчёта, а после расчёта выдаёт результаты расчёта - графики в виде семейств ЧХ. Основной интерфейс программы состоит из трёх вкладок:
- «Г азовая цепь», в которой задаются
конфигурация и параметры газовой цепи, а также отображается её частотная характеристика для текущих значений заданных параметров;
- «Расчёт», с помощью которой формируется множество изменяемых параметров, задаётся имя и путь файла базы данных для хранения результатов расчётов и запускается расчёт;
- «Анализ», предназначена для считывания результатов расчёта из базы данных и отображения в виде графиков.
В программе можно задать конфигурацию цепи, перечислив элементы в списке в той последовательности, в которой они соединяются. Верх списка соответствует подключению к точке более высокого давления объекта измерения (например, вход мерной диафрагмы расходомера), нижняя строка - подключению к точке более низкого давления (например, выход мерной диафрагмы расходомера). Программа построена так, что в библиотеку элементов газовой цепи достаточно просто, дописывая текст программы, добавлять новые математические модели (например, глушители, гасители, фильтры и т. д.) в виде классов объектов, построенных по определённым правилам. Длина газовой цепи не является фиксированной, модель можно редактировать, добавляя в конец цепи, вставляя и удаляя элементы, а также изменяя порядок или тип элементов. Частотные характеристики цепи (АЧХ, ФЧХ) можно рассчитывать как по давлению, так и по перепаду давления.
Выходы газовой измерительной цепи как по давлению, так и по перепаду давления определяются заданными контрольными точками между теми элементами, в сечениях между которыми необходимо контролировать давление и расход. Количество контрольных точек, добавляемых в схему, не ограничено. Перед расчётом характеристик необходимо задать, из каких контрольных точек необходимо сохранить данные.
Рис. 1. Алгоритм расчёта частотных характеристик газовой измерительной цепи дифференциальной схемы
Изменяющиеся параметры для построения серии характеристик можно задать в закладке «Расчёт». В верхней части панели задаётся тип характеристики (АЧХ, ФЧХ) - величина, по которой строятся: характеристика (давление, перепад давления); частотный диапазон; количество точек характеристики и контрольные точки; значения давлений и расходов; частотные характеристики. На левой центральной панели закладки перечислены все параметры, которые можно сделать изменяющимися, на правой панели - те параметры, которые заданы уже как изменяющиеся. Если необходимо, чтобы нескольким параметрам цепи одновременно присваивалось одно и то же значение, следует создать «группу одновременно изменяющихся параметров» и добавить в неё необходимые параметры. Для каждого параметра задаётся начальное и конечное значения и приращение. Таким образом, определяется множество значений параметров, с которыми будут рассчитываться частотные характеристики. Если задано несколько параметров, то программа рассчитает частотную характеристику для каждого сочетания значений из множеств их значений.
Результаты расчёта сохраняются в базе данных, имя и путь - задаются в строке в нижней части содержимого вкладки. Каждый расчёт можно
комментировать одной строкой, далее результаты нужного расчёта будут выбираться по данной строке с комментарием. Кнопка «Расчёт» запускает расчёт характеристик.
Анализ рассчитанных данных предполагает построение одного графика или серии графиков частотных
характеристик. С помощью интерфейсных элементов на вкладке «Анализ» создаётся запрос базы данных, которая возвращает данные, необходимые для графического отображения частотных характеристик. Частотные характеристики можно хранить для последующей обработки, анализа и выбора параметров корректирующих элементов, обеспе-
чивающих требуемую динамическую точность измерения давления и перепада давления.
Библиографический список
1. Климнюк, Ю.И. Оптимизация конструкций компрессоров авиационных ГТД с учетом воздействия эксплуатационных неоднородностей воздушного потока [Текст] / Ю.И. Климнюк. - Самара: Изд-во Самар. гос. аэрокосм.ун-та, 2001. -248 с.
2. Кашапов, И.Д. Разработка методов и средств экспериментального исследования влияния пульсаций давления на погрешность газовых расходомеров. [Текст] / дис. . канд. техн. наук / И. Д. Кашапов. - Самара, 2001. - 190 с.
3. Гимадиев, А.Г. Методика расчёта пульсаций давления и расхода газа в измерительной цепи диафрагменного расходомера [Текст] / А.Г. Гимадиев, А. А. Иго-нин, В.Н. Илюхин // Вестн. Самар. гос. аэрокосм. ун-та. - 2006. - №2. - С. 270-275.
4. Шорин, В.П. Гидравлические и газовые цепи передачи информации [Текст] / В.П. Шорин, А.Г. Гимадиев, Н.Д. Быстров. - М.: Машиностроение, 2000. -328 с.
5. http://www.flrebirdsql .ога [Электронный ресурс].
SIMULATOR, ALGORITHM AND SOFTWARE FOR THE COMPUTATION OF DIFFERENTIAL GAS CIRCUITS FREQUENCY RESPONSE
© 2012 A. A. Igonin, A. G. Gimadiev
Samara State Aerospace University named after academician S.P. Korolyov (National Research University)
An algorithm and software for the computation of differential gas circuits frequency response series based on the mathematical model of electrical analogies is described .The basic formulas of the simulator, a description of the developed software functioning algorithm and the basic techniques of using the software are mentioned.
Simulator, frequency response differential gas circuits, thin pipes, volume, membrane, throttling element.
Информация об авторах
Игонин Андрей Александрович, инженер кафедры автоматических систем энергетических установок, Самарский государственный аэрокосмический университет имени академика С.П. Королёва (национальный исследовательский университет). Email: andy [email protected]. Область научных интересов: мехатроника и робтотехника.
Гимадиев Асгат Гатьятович, доктор технических наук, профессор кафедры автоматических систем энергетических установок, Самарский государственный аэрокосмический университет имени академика С.П. Королёва (национальный исследовательский университет). Е-mail: gimadiev [email protected]. Область научных интересов: динамика пневмогидравлических и топливных систем энергетических установок.
Igonin Andrey Aleksandrovich, engineer, the department of automatic systems of power plants, Samara State Aerospace University named after academician S.P. Korolyov (National Research University). E-mail: andy [email protected]. Area of research: mechatronics and robotics.
Gimadiev Asgat Gatyatovich, doctor of technical science, professor, the department of automatic systems of power plants, Samara State Aerospace University named after academician S.P. Korolyov (National Research University). Е-mail: gimadiev [email protected]. Area of research: dynamics of power plant pneumohydraulic and fuel systems.