УСТАНОВКА ДЛЯ ОПРЕДЕЛЕНИЯ ДИНАМИЧЕСКИХ ХАРАКТЕРИСТИК СЕНСОРОВ ТЕПЛОВОГО ПОТОКА
ТЕПЛОВЫЕ РЕЖИМЫ ПРИБОРОВ И СИСТЕМ
УДК 536.629.7
Н. В. ПИЛИПЕНКО, Г. В. ПОЛЬЩИКОВ, И. А. СИВАКОВ
УСТАНОВКА ДЛЯ ОПРЕДЕЛЕНИЯ ДИНАМИЧЕСКИХ ХАРАКТЕРИСТИК СЕНСОРОВ ТЕПЛОВОГО ПОТОКА
Представлена конструкция стенда для задания нестационарного теплового потока и измерения динамических характеристик датчиков. Приведены результаты эксперимента, показывающие возможность применения батарейного приемника теплового потока (датчика Геращенко) в нестационарной теплометрии.
Ключевые слова: нестационарная теплометрия, датчик Геращенко, обратная задача теплопроводности, параметрическая идентификация.
Решение задач прикладной теплометрии для таких отраслей науки и техники, как теплоэнергетика, металлургия, авиа-, судо- и ракетостроение, мощное приборостроение и пр., определяет дальнейший потенциал их развития. В некоторых практически важных случаях интерес представляет исследование переходных процессов в изучаемых объектах, когда измерения постоянных или переменных во времени плотностей тепловых потоков выполняются в нестационарном режиме работы приемников теплового потока (ПТП).
Один из интенсивно развивающихся в последние 30—40 лет подходов к этой проблеме — использование методов решения обратных задач теплопроводности (ОЗТ) и различных алгоритмов оптимизации. Особенностями данного подхода являются расчетное определение (восстановление) теплового потока по измеряемым температурам в отдельных точках ПТП и возможность определения температур, которые по каким-либо причинам нельзя непосредственно измерить, уточнение глубин заделки термопар, начальных распределений температур и др.
Существуют различные методы решения обратных задач. Выбор наиболее универсального, помехозащищенного и вычислительно эффективного метода решения для определенного класса задач является достаточно сложной проблемой.
В работах Д. Ф. Симбирского, А. В. Олейника, Н. В. Пилипенко, Дж. Бека и других авторов предложено использовать рекуррентные (последовательные) методы параметрической идентификации, базирующиеся, в частности, на модифицированном алгоритме цифрового фильтра Калмана. Изначально разработанные для использования в измерительных системах реального времени, данные методы доказали свою эффективность при решении ряда граничных задач для различных типов ПТП [1, 2]. Оценка возможностей подхода на основе решения обратных задач теплопроводности показала [3—5], что его применение позволяет получить приемлемые для практики результаты.
В настоящей статье рассматривается метод восстановления нестационарного теплового потока на основе параметрической идентификации модели теплопереноса в ПТП, представлена конструкция разработанной установки для задания нестационарных тепловых потоков и определения динамических характеристик ПТП, а также приведены результаты эксперимента по восстановлению плотности теплового потока с использованием тепломера Геращенко.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 5
72 Н. В. Пилипенко, Г. В. Польщиков, И. А. Сиваков
Рассмотрим кратко способ получения оптимальных оценок теплового потока q(τ) . Ди-
намика тепловых процессов в общем случае описывается одним или несколькими параболи-
ческими уравнениями в частных производных с соответствующими граничными условиями и
представляется в виде температурного поля в некоторой многосвязной области. Процесс теп-
лопереноса в ПТП может быть описан с использованием дифференциально-разностной моде-
ли (ДРМ), которая в общем случае имеет следующий вид [1, 2]:
d dτ
T(
τ)
=
FT
(
τ
)
+
GU
(
τ)
,
где F — n × n -матрица обратных связей, G — n × 2 -матрица управления; вектор управления
U (τ) , в частности для условий 2-го рода на торцах ПТП, определяется как
U (τ) = ⎣⎡q1 (τ) q2 (τ)⎤⎦T .
Тепловая схема, подробный вывод и топология ДРМ батарейного ПТП, а также структура матриц управления и обратных связей приведены в работе [6].
В приемниках теплового потока измерению подлежат либо температуры в отдельных точках, либо их разности, либо среднеобъемные температуры, что отражено в m × n -матрице измерений D универсальной модели измерений:
Yk = DTk + εk ,
где Yk — m ×1 -вектор измерений в k-й момент времени, εk — m ×1 -вектор случайных погрешностей измерений температур.
Искомый тепловой поток задается уравнением [4]
r
∑q(τ) = q jϕ j (τ) , j=1
где ϕ j (τ) — система базисных функций времени; q j — априори неизвестные коэффициенты,
которые объединяются в r ×1-вектор искомых параметров Qz = q1 q2 " qr T .
В качестве базисной функции используются B-сплайны 1-го порядка.
Такую аппроксимацию q (τ) называют параметризацией ОЗТ. Тогда задача восстанов-
ления q (τ) сводится к последовательному получению оптимальных оценок Qˆ z,l вектора ис-
комых параметров Qz путем минимизации по Qz квадратичной функции невязки:
( ) ( )∑l
Φ(Qz ) =
Yk − Yˆ k (Qz ) Т R−1 Yk − Yˆ k (Qz ) ,
k =1
где Yˆ k (Qz ) — модельный вектор измерений, являющийся аналогом вектора Yk и рассчиты-
ваемый по ДРМ теплопереноса в ПТП для различных значений вектора Qz ; R — ковариа-
ционная m × m -матрица вектора εk случайных погрешностей.
Для получения оптимальных оценок Qk+1 вектора искомых параметров Qz в (k +1) -й
момент времени используется рекуррентный алгоритм фильтра Калмана по искомым параметрам [3—5]:
Qk+1 = Qˆ k + Kk+1 ⎡⎣Yk+1 − Yˆ k+1(Qˆ k )⎤⎦ ,
( )Kk+1
=
Pk
H
Т k
H
k
Pk
H
Т k
+R
−1 ,
Pk+1 = Pk − Kk+1Hk Pk ,
где Pk , Pk+1 — ковариационные матрицы погрешностей оценок параметров для моментов времени τk = k∆τ и τk+1 = (k +1)∆τ ; Hk — матрица коэффициентов чувствительности изме-
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 5
Установка для определения характеристик сенсоров теплового потока
73
ряемой температуры ПТП к изменению искомых параметров в момент времени τk+1 ; Kk —
весовая матрица. Теоретическая возможность использования высокоинерционных датчиков при исследо-
вании быстропротекающих процессов показана в работе [6], однако полученные в ней результаты требуют экспериментальной проверки.
В качестве объекта исследования был выбран батарейный ПТП (тепломер Геращенко) как наиболее распространенный и обладающий рядом положительных характеристик. Для оценки возможности применения батарейных ПТП в нестационарной теплометрии была создана установка для задания нестационарного потока, структурная схема которой приведена на рис. 1.
КГЛ ОС М
Д
К СКТ
БП
220 В 50 Гц
СУ1
СУ2
220 В 50 Гц
PCS-500
Персональный Щ14-13 компьютер
Рис. 1
В качестве излучателя использовалась галогенная лампа типа КГЛ-100-12V с оптиче-
ской системой (ОС) в виде эллиптического отражателя. Использование такого отражателя по-
зволяет сфокусировать на объекте значительно большую часть излучаемой лампой мощности
по сравнению с линзовыми конденсорами. Питание лампы осуществлялось через трансфор-
маторный блок (БП). Между оптической системой и исследуемым датчиком (Д) располагался
модулятор (М), используемый при проведении испытаний с модулированным потоком излуче-
ния. Согласующий усилитель (СУ1) обеспечивает работу датчика, сигнал которого регистри-
руется компьютерным цифровым осциллографом типа PCS500 фирмы “Velleman” (Бельгия).
Так как спектр излучения лампы находится в широком спектральном диапазоне (от ультрафиолетового
12
до дальнего инфракрасного), а пучок излучения, фор-
3
мируемый эллиптическим отражателем, имеет сложную пространственную структуру, то измерить тепло-
4
вой поток лампы в плоскости исследуемого датчика стандартными средствами не представляется возмож-
5
ным. Поэтому был разработан и изготовлен калори-
метр (К) с квазиобъемным поглотителем и апертурной
Рис. 2
зеркальной диафрагмой. Конструктивная схема калориметра приведена на рис. 2. Измери-
тельное ядро 1 калориметра, представляющее собой алюминиевую пластинку, на одну сторо-
ну которой нанесено поглощающее покрытие (остальные грани ядра свободны), закреплено с
помощью двух кернов 4 с большим тепловым сопротивлением относительно боковых стоек 3,
жестко установленных на основании 5. Коэффициент отражения внутренних поверхностей
всех деталей составляет не менее 0,95. Зеркальная диафрагма 2 формирует входной зрачок
калориметра 10×10 мм2. Для уменьшения тепловых потерь за счет конвективных потоков
расстояние между конструктивными элементами не превышает 2—3 мм.
Для обеспечения эффективного поглощения излучения во всем спектральном диапазоне
использовалось комбинированное поглощающее покрытие. Первоначально наносилось глу-
бокое аноксовое покрытие, а затем наносился лак Паркинсона с медленной сушкой в малом
свободном объеме. Эффективность такого покрытия характеризуется коэффициентом погло-
щения 0,98—0,99.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 5
74 Н. В. Пилипенко, Г. В. Польщиков, И. А. Сиваков
Система контроля температуры (СКТ) — см. рис. 1 — содержит медьконстантановые термопары. Согласующий усилитель (СУ2) позволяет использовать для регистрации сигнала термопар любой вольтметр постоянного тока; в данном случае выбран вольтметр Щ14-13.
Галогенная лампа, эллиптический отражатель и модулятор потока излучения конструктивно объединены в блок излучателя, общий вид которого показан на рис. 3. Основанием блока служит алюминиевая плита 1, обеспечивающая жесткость всей конструкции. На оси 3 расположен двухлопастной диск модулятора 2, который приводится во вращение синхронным двигателем типа Г32. На оси двигателя установлен шкив, обеспечивающий различные скорости вращения диска модулятора. Частота модуляции (от 0,1 до 10 Гц) зависит от количества оборотов двигателя, количества лопастей модулятора и диаметра шкива двигателя.
1
2
3
Рис. 3
При разработке электрической схемы установки учитывалось влияние помех, шумов и температурного дрейфа всех ее элементов. При этом в первом каскаде использовались операционный усилитель с автоматической балансировкой начального смещения и малым уровнем шума, система гальванической развязки для снижения уровня помех, а также дополнительный каскад усиления.
После сборки и юстировки установки с помощью калориметра определялась плотность теплового потока, выполнялась стабилизация галогенной лампы в течение длительного времени, после чего проводилась 20-секундная экспозиция входного окна калориметра и записывалась динамика изменения температуры ядра. На основе полученных данных установлено, что плотность теплового потока, генерируемого галогенной лампой, составляет порядка q0=6·104 Вт/м2.
Так как форма модуляции близка к синусоидальному виду, можно считать, что при включенном модуляторе действующий тепловой поток определяется соотношением
q(τ) = 0,5q0 (1+ sin ωτ) , где ω — частота вращения диска модулятора.
В ходе испытаний плоскость установки приемника теплового потока совмещалась с плоскостью, в которой находилась диафрагма калориметра, проводилась юстировка ПТП относительно светового пятна излучателя и определялась постоянная времени датчика — в рассматриваемом случае равная примерно 5 с. В ходе эксперимента фиксировалось изменение перегрева (∆Т) по толщине батарейного ПТП во времени. Затем датчик подвергался воздействию модулированного лучистого потока того же источника при различных частотах модуляции. Между экспериментами выдерживалась пауза, достаточная для остывания датчика и конструктивных элементов установки.
Полученные экспериментальные результаты сохранялись в компьютере, после чего производилась их обработка с помощью программного комплекса “Heat Flow”, разработанного на кафедре компьютерной теплофизики и энергофизического мониторинга СанктПетербургского национального исследовательского университета информационных технологий, механики и оптики. На рис. 4, а, б соответственно в виде графиков зависимостей ∆Т(τ) и q(τ) приведены результаты эксперимента по восстановлению плотности теплового потока,
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 5
Установка для определения характеристик сенсоров теплового потока
75
изменяющегося с частотой 2 Гц, с использованием батарейного приемника теплового потока с постоянной времени порядка 5 с.
а) ∆Т, К 8 6
4
2
б) q, Вт/м2 6000
4000
2000
00
0 0,5 1,0 1,5 2,0 2,5 3,0 3,5 τ, с
0 0,5 1,0 1,5 2,0 2,5 3,0 3,5 τ, с
Рис. 4
Проведенные экспериментальные исследования батарейного датчика показали возмож-
ность его использования в нестационарной теплометрии. Согласно полученным результатам
предложенный метод параметрической идентификации при восстановлении плотности теп-
лового потока позволяет существенно расширить возможности некоторых видов датчиков.
СПИСОК ЛИТЕРАТУРЫ
1. Пилипенко Н. В. Методы параметрической идентификации в нестационарной теплометрии. Часть 1 // Изв. вузов. Приборостроение. 2003. Т. 46, № 8. С. 50—54.
2. Пилипенко Н. В. Методы параметрической идентификации в нестационарной теплометрии. Часть 2 // Там же. 2003. Т. 46, № 10. С. 67—71.
3. Кириллов К. В., Пилипенко Н. В. Алгоритмы программ для решения прямых и обратных задач теплопроводности при использовании дифференциально-разностных моделей // Науч.-техн. вестн. СПбГУ ИТМО. 2010. № 5. С. 106—110.
4. Пилипенко Н. В., Сиваков И. А. Метод определения нестационарного теплового потока и теплопроводности путем параметрической идентификации // Измерительная техника. 2011. № 3. С. 48—51.
5. Pilipenko N. V., Sivakov I. A. A method of determining nonstationary heat flux and heat conduction using parametric identification // Measurement Techniques. 2011. Vol. 54, N 3. P. 318—323.
6. Пилипенко Н. В., Лазуренко Н. В., Лебедев П. В. Параметрическая идентификация нестационарных тепловых потоков с помощью тепломера типа „вспомогательная стенка“ // Изв. вузов. Приборостроение. 2005. Т. 48, № 9. С. 47—51.
Сведения об авторах
Николай Васильевич Пилипенко
— д-р техн. наук, профессор; Санкт-Петербургский национальный ис-
следовательский университет информационных технологий, механи-
ки и оптики, кафедра компьютерной теплофизики и энергофизиче-
ского мониторинга; E-mail: Pilipenko38@mail.ru
Георгий Владимирович Польщиков — канд. техн. наук, доцент; Санкт-Петербургский национальный иссле-
довательский университет информационных технологий, механики и
оптики, кафедра компьютерной теплофизики и энергофизического
мониторинга; E-mail: polschikov@list.ru
Иван Анатольевич Сиваков
— аспирант; Санкт-Петербургский национальный исследовательский
университет информационных технологий, механики и оптики, ка-
федра компьютерной теплофизики и энергофизического мониторин-
га; E-mail: I.A.Sivakov@gmail.com; SivakovIA@mail.ru
Рекомендована кафедрой компьютерной теплофизики и энергофизического мониторинга
Поступила в редакцию 04.05.12 г.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 5
УДК 536.629.7
Н. В. ПИЛИПЕНКО, Г. В. ПОЛЬЩИКОВ, И. А. СИВАКОВ
УСТАНОВКА ДЛЯ ОПРЕДЕЛЕНИЯ ДИНАМИЧЕСКИХ ХАРАКТЕРИСТИК СЕНСОРОВ ТЕПЛОВОГО ПОТОКА
Представлена конструкция стенда для задания нестационарного теплового потока и измерения динамических характеристик датчиков. Приведены результаты эксперимента, показывающие возможность применения батарейного приемника теплового потока (датчика Геращенко) в нестационарной теплометрии.
Ключевые слова: нестационарная теплометрия, датчик Геращенко, обратная задача теплопроводности, параметрическая идентификация.
Решение задач прикладной теплометрии для таких отраслей науки и техники, как теплоэнергетика, металлургия, авиа-, судо- и ракетостроение, мощное приборостроение и пр., определяет дальнейший потенциал их развития. В некоторых практически важных случаях интерес представляет исследование переходных процессов в изучаемых объектах, когда измерения постоянных или переменных во времени плотностей тепловых потоков выполняются в нестационарном режиме работы приемников теплового потока (ПТП).
Один из интенсивно развивающихся в последние 30—40 лет подходов к этой проблеме — использование методов решения обратных задач теплопроводности (ОЗТ) и различных алгоритмов оптимизации. Особенностями данного подхода являются расчетное определение (восстановление) теплового потока по измеряемым температурам в отдельных точках ПТП и возможность определения температур, которые по каким-либо причинам нельзя непосредственно измерить, уточнение глубин заделки термопар, начальных распределений температур и др.
Существуют различные методы решения обратных задач. Выбор наиболее универсального, помехозащищенного и вычислительно эффективного метода решения для определенного класса задач является достаточно сложной проблемой.
В работах Д. Ф. Симбирского, А. В. Олейника, Н. В. Пилипенко, Дж. Бека и других авторов предложено использовать рекуррентные (последовательные) методы параметрической идентификации, базирующиеся, в частности, на модифицированном алгоритме цифрового фильтра Калмана. Изначально разработанные для использования в измерительных системах реального времени, данные методы доказали свою эффективность при решении ряда граничных задач для различных типов ПТП [1, 2]. Оценка возможностей подхода на основе решения обратных задач теплопроводности показала [3—5], что его применение позволяет получить приемлемые для практики результаты.
В настоящей статье рассматривается метод восстановления нестационарного теплового потока на основе параметрической идентификации модели теплопереноса в ПТП, представлена конструкция разработанной установки для задания нестационарных тепловых потоков и определения динамических характеристик ПТП, а также приведены результаты эксперимента по восстановлению плотности теплового потока с использованием тепломера Геращенко.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 5
72 Н. В. Пилипенко, Г. В. Польщиков, И. А. Сиваков
Рассмотрим кратко способ получения оптимальных оценок теплового потока q(τ) . Ди-
намика тепловых процессов в общем случае описывается одним или несколькими параболи-
ческими уравнениями в частных производных с соответствующими граничными условиями и
представляется в виде температурного поля в некоторой многосвязной области. Процесс теп-
лопереноса в ПТП может быть описан с использованием дифференциально-разностной моде-
ли (ДРМ), которая в общем случае имеет следующий вид [1, 2]:
d dτ
T(
τ)
=
FT
(
τ
)
+
GU
(
τ)
,
где F — n × n -матрица обратных связей, G — n × 2 -матрица управления; вектор управления
U (τ) , в частности для условий 2-го рода на торцах ПТП, определяется как
U (τ) = ⎣⎡q1 (τ) q2 (τ)⎤⎦T .
Тепловая схема, подробный вывод и топология ДРМ батарейного ПТП, а также структура матриц управления и обратных связей приведены в работе [6].
В приемниках теплового потока измерению подлежат либо температуры в отдельных точках, либо их разности, либо среднеобъемные температуры, что отражено в m × n -матрице измерений D универсальной модели измерений:
Yk = DTk + εk ,
где Yk — m ×1 -вектор измерений в k-й момент времени, εk — m ×1 -вектор случайных погрешностей измерений температур.
Искомый тепловой поток задается уравнением [4]
r
∑q(τ) = q jϕ j (τ) , j=1
где ϕ j (τ) — система базисных функций времени; q j — априори неизвестные коэффициенты,
которые объединяются в r ×1-вектор искомых параметров Qz = q1 q2 " qr T .
В качестве базисной функции используются B-сплайны 1-го порядка.
Такую аппроксимацию q (τ) называют параметризацией ОЗТ. Тогда задача восстанов-
ления q (τ) сводится к последовательному получению оптимальных оценок Qˆ z,l вектора ис-
комых параметров Qz путем минимизации по Qz квадратичной функции невязки:
( ) ( )∑l
Φ(Qz ) =
Yk − Yˆ k (Qz ) Т R−1 Yk − Yˆ k (Qz ) ,
k =1
где Yˆ k (Qz ) — модельный вектор измерений, являющийся аналогом вектора Yk и рассчиты-
ваемый по ДРМ теплопереноса в ПТП для различных значений вектора Qz ; R — ковариа-
ционная m × m -матрица вектора εk случайных погрешностей.
Для получения оптимальных оценок Qk+1 вектора искомых параметров Qz в (k +1) -й
момент времени используется рекуррентный алгоритм фильтра Калмана по искомым параметрам [3—5]:
Qk+1 = Qˆ k + Kk+1 ⎡⎣Yk+1 − Yˆ k+1(Qˆ k )⎤⎦ ,
( )Kk+1
=
Pk
H
Т k
H
k
Pk
H
Т k
+R
−1 ,
Pk+1 = Pk − Kk+1Hk Pk ,
где Pk , Pk+1 — ковариационные матрицы погрешностей оценок параметров для моментов времени τk = k∆τ и τk+1 = (k +1)∆τ ; Hk — матрица коэффициентов чувствительности изме-
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 5
Установка для определения характеристик сенсоров теплового потока
73
ряемой температуры ПТП к изменению искомых параметров в момент времени τk+1 ; Kk —
весовая матрица. Теоретическая возможность использования высокоинерционных датчиков при исследо-
вании быстропротекающих процессов показана в работе [6], однако полученные в ней результаты требуют экспериментальной проверки.
В качестве объекта исследования был выбран батарейный ПТП (тепломер Геращенко) как наиболее распространенный и обладающий рядом положительных характеристик. Для оценки возможности применения батарейных ПТП в нестационарной теплометрии была создана установка для задания нестационарного потока, структурная схема которой приведена на рис. 1.
КГЛ ОС М
Д
К СКТ
БП
220 В 50 Гц
СУ1
СУ2
220 В 50 Гц
PCS-500
Персональный Щ14-13 компьютер
Рис. 1
В качестве излучателя использовалась галогенная лампа типа КГЛ-100-12V с оптиче-
ской системой (ОС) в виде эллиптического отражателя. Использование такого отражателя по-
зволяет сфокусировать на объекте значительно большую часть излучаемой лампой мощности
по сравнению с линзовыми конденсорами. Питание лампы осуществлялось через трансфор-
маторный блок (БП). Между оптической системой и исследуемым датчиком (Д) располагался
модулятор (М), используемый при проведении испытаний с модулированным потоком излуче-
ния. Согласующий усилитель (СУ1) обеспечивает работу датчика, сигнал которого регистри-
руется компьютерным цифровым осциллографом типа PCS500 фирмы “Velleman” (Бельгия).
Так как спектр излучения лампы находится в широком спектральном диапазоне (от ультрафиолетового
12
до дальнего инфракрасного), а пучок излучения, фор-
3
мируемый эллиптическим отражателем, имеет сложную пространственную структуру, то измерить тепло-
4
вой поток лампы в плоскости исследуемого датчика стандартными средствами не представляется возмож-
5
ным. Поэтому был разработан и изготовлен калори-
метр (К) с квазиобъемным поглотителем и апертурной
Рис. 2
зеркальной диафрагмой. Конструктивная схема калориметра приведена на рис. 2. Измери-
тельное ядро 1 калориметра, представляющее собой алюминиевую пластинку, на одну сторо-
ну которой нанесено поглощающее покрытие (остальные грани ядра свободны), закреплено с
помощью двух кернов 4 с большим тепловым сопротивлением относительно боковых стоек 3,
жестко установленных на основании 5. Коэффициент отражения внутренних поверхностей
всех деталей составляет не менее 0,95. Зеркальная диафрагма 2 формирует входной зрачок
калориметра 10×10 мм2. Для уменьшения тепловых потерь за счет конвективных потоков
расстояние между конструктивными элементами не превышает 2—3 мм.
Для обеспечения эффективного поглощения излучения во всем спектральном диапазоне
использовалось комбинированное поглощающее покрытие. Первоначально наносилось глу-
бокое аноксовое покрытие, а затем наносился лак Паркинсона с медленной сушкой в малом
свободном объеме. Эффективность такого покрытия характеризуется коэффициентом погло-
щения 0,98—0,99.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 5
74 Н. В. Пилипенко, Г. В. Польщиков, И. А. Сиваков
Система контроля температуры (СКТ) — см. рис. 1 — содержит медьконстантановые термопары. Согласующий усилитель (СУ2) позволяет использовать для регистрации сигнала термопар любой вольтметр постоянного тока; в данном случае выбран вольтметр Щ14-13.
Галогенная лампа, эллиптический отражатель и модулятор потока излучения конструктивно объединены в блок излучателя, общий вид которого показан на рис. 3. Основанием блока служит алюминиевая плита 1, обеспечивающая жесткость всей конструкции. На оси 3 расположен двухлопастной диск модулятора 2, который приводится во вращение синхронным двигателем типа Г32. На оси двигателя установлен шкив, обеспечивающий различные скорости вращения диска модулятора. Частота модуляции (от 0,1 до 10 Гц) зависит от количества оборотов двигателя, количества лопастей модулятора и диаметра шкива двигателя.
1
2
3
Рис. 3
При разработке электрической схемы установки учитывалось влияние помех, шумов и температурного дрейфа всех ее элементов. При этом в первом каскаде использовались операционный усилитель с автоматической балансировкой начального смещения и малым уровнем шума, система гальванической развязки для снижения уровня помех, а также дополнительный каскад усиления.
После сборки и юстировки установки с помощью калориметра определялась плотность теплового потока, выполнялась стабилизация галогенной лампы в течение длительного времени, после чего проводилась 20-секундная экспозиция входного окна калориметра и записывалась динамика изменения температуры ядра. На основе полученных данных установлено, что плотность теплового потока, генерируемого галогенной лампой, составляет порядка q0=6·104 Вт/м2.
Так как форма модуляции близка к синусоидальному виду, можно считать, что при включенном модуляторе действующий тепловой поток определяется соотношением
q(τ) = 0,5q0 (1+ sin ωτ) , где ω — частота вращения диска модулятора.
В ходе испытаний плоскость установки приемника теплового потока совмещалась с плоскостью, в которой находилась диафрагма калориметра, проводилась юстировка ПТП относительно светового пятна излучателя и определялась постоянная времени датчика — в рассматриваемом случае равная примерно 5 с. В ходе эксперимента фиксировалось изменение перегрева (∆Т) по толщине батарейного ПТП во времени. Затем датчик подвергался воздействию модулированного лучистого потока того же источника при различных частотах модуляции. Между экспериментами выдерживалась пауза, достаточная для остывания датчика и конструктивных элементов установки.
Полученные экспериментальные результаты сохранялись в компьютере, после чего производилась их обработка с помощью программного комплекса “Heat Flow”, разработанного на кафедре компьютерной теплофизики и энергофизического мониторинга СанктПетербургского национального исследовательского университета информационных технологий, механики и оптики. На рис. 4, а, б соответственно в виде графиков зависимостей ∆Т(τ) и q(τ) приведены результаты эксперимента по восстановлению плотности теплового потока,
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 5
Установка для определения характеристик сенсоров теплового потока
75
изменяющегося с частотой 2 Гц, с использованием батарейного приемника теплового потока с постоянной времени порядка 5 с.
а) ∆Т, К 8 6
4
2
б) q, Вт/м2 6000
4000
2000
00
0 0,5 1,0 1,5 2,0 2,5 3,0 3,5 τ, с
0 0,5 1,0 1,5 2,0 2,5 3,0 3,5 τ, с
Рис. 4
Проведенные экспериментальные исследования батарейного датчика показали возмож-
ность его использования в нестационарной теплометрии. Согласно полученным результатам
предложенный метод параметрической идентификации при восстановлении плотности теп-
лового потока позволяет существенно расширить возможности некоторых видов датчиков.
СПИСОК ЛИТЕРАТУРЫ
1. Пилипенко Н. В. Методы параметрической идентификации в нестационарной теплометрии. Часть 1 // Изв. вузов. Приборостроение. 2003. Т. 46, № 8. С. 50—54.
2. Пилипенко Н. В. Методы параметрической идентификации в нестационарной теплометрии. Часть 2 // Там же. 2003. Т. 46, № 10. С. 67—71.
3. Кириллов К. В., Пилипенко Н. В. Алгоритмы программ для решения прямых и обратных задач теплопроводности при использовании дифференциально-разностных моделей // Науч.-техн. вестн. СПбГУ ИТМО. 2010. № 5. С. 106—110.
4. Пилипенко Н. В., Сиваков И. А. Метод определения нестационарного теплового потока и теплопроводности путем параметрической идентификации // Измерительная техника. 2011. № 3. С. 48—51.
5. Pilipenko N. V., Sivakov I. A. A method of determining nonstationary heat flux and heat conduction using parametric identification // Measurement Techniques. 2011. Vol. 54, N 3. P. 318—323.
6. Пилипенко Н. В., Лазуренко Н. В., Лебедев П. В. Параметрическая идентификация нестационарных тепловых потоков с помощью тепломера типа „вспомогательная стенка“ // Изв. вузов. Приборостроение. 2005. Т. 48, № 9. С. 47—51.
Сведения об авторах
Николай Васильевич Пилипенко
— д-р техн. наук, профессор; Санкт-Петербургский национальный ис-
следовательский университет информационных технологий, механи-
ки и оптики, кафедра компьютерной теплофизики и энергофизиче-
ского мониторинга; E-mail: Pilipenko38@mail.ru
Георгий Владимирович Польщиков — канд. техн. наук, доцент; Санкт-Петербургский национальный иссле-
довательский университет информационных технологий, механики и
оптики, кафедра компьютерной теплофизики и энергофизического
мониторинга; E-mail: polschikov@list.ru
Иван Анатольевич Сиваков
— аспирант; Санкт-Петербургский национальный исследовательский
университет информационных технологий, механики и оптики, ка-
федра компьютерной теплофизики и энергофизического мониторин-
га; E-mail: I.A.Sivakov@gmail.com; SivakovIA@mail.ru
Рекомендована кафедрой компьютерной теплофизики и энергофизического мониторинга
Поступила в редакцию 04.05.12 г.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 5