книги из ГПНТБ / Коротков П.А. Динамические контактные измерения тепловых величин
.pdfПриведем условия, при которых имеют место те или иные осо бенности измерительного процесса. Так, при Bi «сС 1 измеритель ный 'процесс характеризуется особенностями, свойственными первой зоне зависимости Lo 3 (Fo). Если Bi > 1, то имеют место особенности ее третьей зоны. И, наконец, при Bi «=; 1 приемный преобразователь вносит максимум искажения в контролируемый процесс. В этом случае проявляются особенности второй зоны зависимости Lo 3 (Fo).
Применив к выражению (1.99) преобразование Фурье, получим
|
|
|
|
Bi + 1 |
Lo3 (/со) = |
|
|
|
|
|
т ехр / |
{arctg |
1 |
|
I2 |
/ |
1 + Bi |
|
||
|
|
|
' aw V Bi |
|||
1 |
ш |
I |
В |
|
|
|
|
|
|
|
|
|
(1.100)
|
|
Bi-W |
0 |
<1 |
i4S |
|
Рис. 18. Фазово-частотная характеристика измерительной среды по тепловому потоку при различных значениях параметра Bi
Условия, при которых имеют место те или иные проявления измерительного процесса, будут те же, что и относящиеся к вы ражению (1.99). Одноемкостная схема замещения в основном находит применение при интенсивном теплообмене между средой и первичным преобразователем. При относительно умеренном теплообмене может быть применена двухъемкостная модель тепло передачи.
Для двухъемкостной схемы замещения измерительной среды будем иметь
4-Lo3 (Bi, |
Fo) = |
|
ехр (- |
|
(1.101) |
где |
|
|
|
|
|
|
/ 2 |
/ |
1 - f Bi \ |
l\'i ( / 11++2BJX2Bix |
\ |
|
ax |
V |
Bi |
ai V Bij |
J' |
(aot) |
|
|
|
|
|
T3 |
= |
|
|
|
|
|
(aaj) |
|
(aa,) |
|
|
40
При т —> О |
|
|
|
|
|
А |
L o 3 ( B i , |
Fo) = |
^ _ |
+ - T - . |
( 1 Л 0 2 ) |
|
|
|
' 3 |
' 4 |
|
Если т = оо , то |
Lo3 (Bi, Fo) = |
0. |
|
||
При а —-> о о и |
т « |
1 0 1 |
; |
|
(1-103) |
^a, Bi,.
|
|
|
|
|
/2 |
|
|
|
|
а |
п |
|
-/f L o 3 ( B i , F o ) = l + — Р а |
|
|||
|
|
|
1 + |
2 |
|
|
|
|
|
a |
|
где B i ^ - |
^ - ; а 1 3 = — ; |
я = |
— . |
|
|
Положив |
в (1.102) а = аг, |
I = |
/ ь |
будем |
иметь, что |
- £ - Lo s ( B i , Fo)->2,5.
При а —> 0 с учетом условия (1.103)
a B i 2 2 \
£ L o , ( B i , Fo) ^ 1 + i +
Bii
где Bi, |
|
Из (1.104) при — в ^ 2 - — 1 найдем, что |
|
- £ - Lo 3 (Bi, Fo) |
|
Положив в (1.104) / = 0 или 1г = 0 в зависимости |
от условий |
опыта, получим расчетную формулу для определения |
- ^ - L o 3 (Bi, |
Fo) применительно к одноемкостной схеме замещения. |
|
Сравнивая результаты исследований зависимости |
Lo 3 (Bi, Fo) |
одноемкостной и двухъемкостной схем замещения, видим, что
значение Lo3 (Bi, Fo) |
одноемкостной |
схемы при Fo — const, |
||||
Bi = |
const |
будет всегда меньше, чем двухъемкостной. |
Это вы |
|||
звано |
влиянием второго звена на первое, которое вызвано |
шунти |
||||
рующим действием последующего звена предыдущим. |
В |
резуль |
||||
тате |
этого |
положение |
характерных |
зон зависимости |
|
Lo3 (Fo) |
двухъемкостной схемы замещения будет смещено в сторону более малых значений параметра Fo по отношению к одноемкостной схеме замещения.
41
Если
^ |
а |
» 1 и |
|
" |
С1РЛ + cpl
а- 2 1а
то процесс теплопередачи в зоне исследования подчиняется особен
ностям первой |
характерной |
области |
изменения Lo 3 (Fo). |
|
|
|
г- |
|
<? |
При |
U Bi « 1 и |
/2 |
~ Г В |
12«1 |
+ 2 |
/ 2 |
|||
|
|
а |
м |
|
|
|
|
|
а
в зоне исследования проявляются особенности третьей области изменения Lo 3 (Fo).
Частотные свойства двухъемкостной схемы замещения измери тельного процесса по тепловому потоку описываются выражением
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
1 + ©- |
|
|
IL |
|
|
|
Lo3(/co) |
= |
|
|
|
|
|
|
||
|
|
|
|
|
|
X |
|
|||
|
|
|
|
|
|
|
a21 |
a Bi |
|
|
|
|
|
|
|
|
|
1 / Г |
|
|
|
|
|
X e x p { - / a r c t g [ ^ ( - ^ |
+ |
^ - ) |
|
(1.105) |
||||
Очевидно, что |
при |
, а |
^> - у - -{- -тр- |
измерительный |
процесс |
|||||
|
|
|
|
|
/j |
1 |
|
|
|
|
подчиняется |
условиям, |
свойственным |
первой |
области |
Lo3 (Fo). |
|||||
Если при этом |
<С |
|
, то |
|
|
|
|
|
|
|
|
|
|
|
|
2 \ т - |
|
|
|
|
|
|
|
|
|
1 + |
- |
|
|
|
|
|
Я Bi |
Lo8 |
(/©) = • |
1 + |
Bj2q2 |
e x p ( - / a r c t g ^ ) . |
|
||||
|
|
|
|
0)2/2 |
|
|
|
|
|
|
В том случае, когда |
-^w~ <С -Щ- + |
~тг- > |
измерительный |
процесс |
||||||
подчиняется |
условиям, |
свойственным |
третьей |
области |
Lo 3 (Fo). |
|||||
При этом |
|
|
|
|
|
|
|
|
|
|
ш ь ° з |
ущ )=е х р Н a r c t |
g Ш"^+ a t ) }} • |
|
42
Из (1.105) видим, что свойства измерительного процесса за висят не только от значений физико-геометрических параметров среды и первичного преобразователя (измерительной среды), но также от частотного диапазона области существенных частот контролируемого теплового потока1 . Поэтому при выборе физикогеометрических параметров первичного преобразователя следует знать частотные свойства контролируемого сигнала. Эти свойства, как правило, устанавливаются исходя из внешнех условий опыта.
В |
заключение |
приведем пример использования рассмотрен |
ных |
характеристик. Пусть по условиям опыта требуется измерить |
|
тепловой поток, |
поступающий от газа к термоприемнику. Этот |
поток должен быть образован конвекцией. Таким образом, в ор ганизации величины теплового потока не должны принимать участие теплофизические свойства термоприемника. Очевидно, что значения физико-геометрических свойств термоприемника должны быть выбраны такими, при которых зависимость Lo 3 (Fo)
будет определять |
первую |
характерную |
зону. |
Известно |
также, |
|||
что Bi = |
Ю - 1 , т и |
= 1 0 - 1 с и |
что процесс теплопередачи в измери |
|||||
тельной схеме подчиняется |
уравнению (1.6) с условиями |
(1.10) — |
||||||
(1.12). Из |
рис. 16 |
заключаем, что Fo ^ |
1 и — |
^ Ю - |
1 . |
Если |
||
материал |
термоприемника — медь, то / = |
3. 10~3 |
м, так как а = |
|||||
= 0,36 м2 /ч. |
|
|
|
|
|
|
|
|
После |
того как определены необходимые |
значения |
|
физико- |
геометрических параметров термоприемника, при которых из меренная величина теплового потока будет соответствовать ис тинной, необходимо для подбора соответствующей аппаратуры оценить величину вторичного сигнала. Условием для выбора ре гистрирующей аппаратуры является обеспечение надежной реги страции вторичного сигнала. Это условие будет выполнено, если уровень полезного сигнала превысит уровень шума. Для опреде ления уровня вторичного сигнала с учетом условий измерений и выбранных значений физико-геометрических параметров термо приемника достаточно воспользоваться переходными или частот ными характеристиками измерительного тракта (см. п. 1).
П р и м е ч а н и е . Нестационарный теплообмен, как правило, |
исследуется |
с учетом уравнения Фурье, записанного в виде (1.5). При этом для |
нахождения |
взаимосвязи вторичного сигнала с истинным (1.6) берется с условиями (1.10)— (1.12). Как показано в [57], соотношение (1.94) предусматривает бесконечную фазовую скорость распространения тепла. Однако в действительности скорость распространения тепла имеет конечное значение [57 ]• Закон Фурье с учетом урав
нения (1.6) с условиями (1.10)—(1.12) |
имеет |
вид |
|
L o 3 ( B i , Fo) + |
т,а d Lo3 (Bi, Fo) |
|
|
7X" |
dFo |
|
|
= Bi exp (Bi2 Fo) erf (Bi l^Fo), |
(1.106) |
1 Контролируемый тепловой поток есть поток тепла, подлежащий определению при измерениях.
43
где ir — модуль сдвига. Решением (1.106) является
J_ |
, |
|
„ |
, |
/ |
Fo / |
2 |
|
X |
- L o 8 |
( B i . |
Fo) = |
e x p ( - ^ | ) X |
|
|||
|
( |
FoFo |
|
|
|
|
|
|
Х-J5Б1+B i |
JJ Bii exp(hiB-i-i~)+ i i ) X |
|
||||||
|
|
|
|
|
2 |
2 |
|
|
|
|
|
о |
|
|
|
|
|
|
|
|
|
Bi V'. |
|
. dx. |
(1.107) |
|
X |
1 — — • |
f |
exp (г2 ) dz |
|||||
|
|
V л |
J |
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
Исходя из начального условия — Lo 3 (Bi, 0) = Bi,
|
|
|
Л |
|
i - |
Lo 3 (Bi,Fo) = Bi exp ( - ^ - ) - f - |
|
||
|
|
Fo |
|
|
X |
|
2 |
Bi]fP° |
(1.108) |
1 |
-== |
exp (—22 ) dz |
||
|
|
Vn |
J |
|
Рассмотрим поведение функции (1.108). Для этого исследуем ее производные:
l-3L°3^ |
|
) = |
Bi2 exp |
(Bi 2 |
Fo)erf |
( B i l A F o ) ( l |
+ |
2 / F o |
/ 2 ' |
|
|
||||
F 0 |
т,а Bi |
, |
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
J |
i |
/ |
/2 |
Fo \ |
|
|
|
|
|
|
|
(1.109) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d |
2 L o , ( B i , |
Fo) |
_ . , |
, „ . „ _ . |
. / . . . / - p - v / , |
, |
2 l / F o / 2 |
\ |
) |
, |
+ |
||||
|
3 - — ; |
|
= B i 4 exp (Bi 2 |
Fo) erf (Bi V Fo) I 1 -f- |
_ _D . |
|
|
||||||||
|
3Fo 2 |
|
|
|
|
|
|
|
|
|
i,a Bi |
|
|
|
|
+ |
B i 2 exp ( B i 2 |
Fo) 2 J |
^ 1 e |
x p ( _ B i |
/ Г о ) ( |
|
+ |
1 ) X |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
axr Bi |
|
|
|
|
||
|
|
Bi2 exp(Bi2 Fo)erf |
(Bi VFO) |
x,a Bi / F o |
|
|
|
|
|
||||||
|
|
|
/ 2 |
/ |
f2 Fo |
\ |
|
|
|
|
|
|
(I.110) |
||
|
|
|
T,a' • exp |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
Численные расчеты и анализ зависимостей (1.109), (1.110) показали, что функция (1.108) имеет экстремумы (максимум и минимум функции) при значениях Fo, весьма близких функции LoH (Fo). Поэтому можно положить, что термоприем ник потребляет максимальное значение тепловой мощности от среды и таким об разом значительно искажает поле температур в зоне измерения, если его физикогеометрические параметры подобраны такими, что при них проявляются особен ности второй характерной зоны зависимости Lo3 (Fo). Аналогичная картина будет иметь место при отличных от (1.10)—(1.12) условиях, но с учетом второй зоны зависимости Lo3 (Fo).
44
Рассмотрим зависимость следующего вида [53]:
^д?(р, х)
|
~Ка(Р)ф |
Т-^—, |
Х = 0. |
(1.Ш) |
|
|
|
|
t(P) |
|
|
Она может быть |
получена |
из |
(1.93), |
если положить |
Bi = оо. |
Таким образом, |
Kq (Fo) = |
Lo3 |
(Bi, Fo), Bi = oo. |
|
|
Функциональную зависимость Kq (Fo) иначе называют вход |
|||||
ной проводимостью первичного |
преобразователя [52]. Приведем |
лЧ
Я1 |
I |
I |
IX |
|
I |
I |
I |
10'' |
10° |
10' |
!0} |
Ю3 |
Ю" |
Ю Fo |
|
|
Рис. |
19. Зависимость -j- |
|
Kq |
от Fo: |
|
|
/ — с учетом условий |
(1.16)—(1.18); |
2 — с учетом |
условий |
||||
|
|
|
(1.13)-(1.15) |
|
|
|
ее выражения применительно к рассмотренным выше условиям теплопередачи.
Значение входной проводимости модели первичного преобра зователя, приведенной на рис. 1, при t = const изменяется во времени следующим образом:
|
|
' |
л п= 0 |
|
Используя лишь первый член ряда (1.112), получаем выраже |
||||
ние для входной проводимости модели первичного |
преобразова |
|||
теля, |
изображенной на рис. 2. |
|
|
|
На |
практике обычно используют не более первых трех чле |
|||
нов ряда (1.112). Это вызвано |
тем, что с увеличением п значения |
|||
членов данного ряда резко |
убывают. |
|
||
На |
рис. 19 приведены |
зависимости вида -у- Kq |
(Fo), полу |
ченные с помощью выражения (1.112). Анализ кривых приводит
45
к выводу, что при Fosg: 10"1 можно воспользоваться лишь первым членом ряда. С увеличением значения параметра п влияние «отра женных тепловых волн» на формирование теплового потока воз растает. Таким образом, значение п должно выбираться в каждом конкретном случае отдельно.
Входная проводимость одноемкостной схемы замещения пер вичного преобразователя изменяется во времени согласно выра
жению |
|
Kq (Fo) = А ехр (—Fo). |
(1.113) |
Этим выражением обычно пользуются при интенсивном тепло обмене между средой и первичным преобразователем или при тех значениях Fo, при которых выполняются условия регуляр ного режима. При значениях параметра Fo, существенно меньших единицы, следует пользоваться выражениями входной проводи мости, учитывающими особенности «работы» первичного преобра зователя как элемента с распределенными параметрами с, р и К.
Входная проводимость двухъемкостной схемы замещения пер вичного преобразователя рассчитывается по формуле
К, (Fo) = |
f { d - d l ) " Р |
- |
exp ( - ^ |
), |
(1.114) |
где |
|
|
|
|
|
Применяя к |
выражениям |
(1.112) — (1-114) |
преобразование |
||
Фурье, получим |
|
|
|
|
|
|
со |
|
|
|
|
|
п = 0 |
|
|
|
(1.115) |
|
— ехр [ —• ( 2 п + 1 ) / ] / ^ ] ( ; |
|
|||
|
|С„(/со)| = УЩ-, |
п = 0; |
|
(1.116) |
|
Са (/'©) | = -J- [ ! - е х Р |
(-ш)] |
|
(1.117) |
||
1 Cq (/со) | = |
|
|
|
(1.118) |
|
|
У(1-ту)+ту |
|
|
|
Использование выражений (1.112)—-(1.114) или (1.115) — (1.118) определяется удобствами наблюдения контролируемых особенностей изучаемого процесса. Выражения (1.112) — (1-П4)
46
применяются главным образом при расчете величины теплового потока, так как Kq (Fo) входит в формулу т для расчета теплового потока
q(Fo) = Kq (0) /(т) +хJ Kq-(х- n)/(Ti)dx{.
0
Исходя из этого и учитывая, что степень сложности выраже ний (1.112) — (1.П4) различная, можно их применять также при упрощениях расчетных операций. Эта задача решается путем сравнения зависимостей Kq (Fo) с учетом заданной величины па раметра Fo и погрешности расчета.
Использование полученных соотношений для расчета величины теплового потока предусматривает знание физико-геометрических свойств измерительной среды. Их значения на практике известны, как правило, с некоторой погрешностью. Величина погрешности нахождения теплового потока рассчитывается по формуле
|
|
|
я |
|
|
|
(1.119) |
|
|
|
|
|
|
|
|
|
|
где |
Aq—абсолютная |
величина |
погрешности |
нахождения тепло |
||||
вого |
потока. |
|
|
|
|
|
|
|
На практике значение Atc |
не всегда известно. В таких случаях |
|||||||
вместо (1.119) |
применяется |
иное |
соотношение: |
|
||||
|
|
|
|
Я |
|
Kq ' |
(1.120) |
|
|
|
|
|
|
|
|
||
где |
AKq—абсолютная |
величина |
погрешности нахождения |
Кч. |
||||
Задаваясь |
значением AFo, |
исходя при |
этом, например, |
из |
возможностей методики измерений таких параметров первичного преобразователя, как а, I, К, можно определить Aq с помощью выражения (1.120). Так, если процесс теплопередачи в зоне из мерения подчиняется уравнению (1.6) с условиями (1.10) —(1.12), то погрешность нахождения величины теплового потока может
быть рассчитана |
по |
выражению |
|
|
|
|
Д Bi ^ |
|
|
|
|
-вГ) |
|
|
|
г |
2ДВ1 (Ш |
»2 ] / |
|
*1 |
х 1 |
+ -вГ+{ж) |
J ( 1 + |
|
я |
|
erf ( B i K F o ) exp(Bi2 Fo ) |
|
|
X erf |
[ ( B i ] / F o ) ( l + ^ ] / l 4 AFoFo — 1. |
(1.121) |
47
Если условия измерений удовлетворяют требованиям первой характерной зоны зависимости Lo 3 (Fo), то, как следует из (1.121),
доминирующим фактором, влияющим на величину - ~ , будет ABi.
Если по условиям опыта |
измерительный |
процесс удовлетворяет |
|||||||||
требованиям |
третьей |
характерной |
зоны |
зависимости |
Lo 3 (Fo), |
||||||
то основным источником погрешности будет AFo. |
|
||||||||||
Если |
измерительный |
процесс |
удовлетворяет |
требованиям |
|||||||
(1.7) - (1 . 9) , |
то |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Fo |
i |
"г Fo |
J |
|
|
|
|
|
|
— ехр |
( я + |
I ) 2 |
|
|
|
|
||
|
|
|
|
|
Д Р о \ |
|
|
|
|||
|
|
|
|
Fo |
|
|
|
|
|||
|
|
|
|
|
~ F o / |
|
|
(1.122) |
|||
|
|
|
|
|
|
|
|
|
|||
|
q |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(n + |
l) 2 |
|
|
|
|
|
2 e x p ( - l o - ) - e x p |
|
|
|
|
|||||
|
|
|
|
Fo |
|
|
|
||||
Если требуется измерить тепловой поток, вызванный конвек |
|||||||||||
цией, но при этом условия измерений удовлетворяют |
требованиям |
||||||||||
третьей |
зоны |
зависимости |
Lo 3 |
(Fo), |
то |
величина |
погрешности, |
||||
с которой находится |
тепловой |
поток, |
вычисляется |
по формуле. |
|||||||
В (1.121) и (1.122) через |
ABi обозначена абсолютная |
величина |
|||||||||
погрешности |
определения |
параметра |
Bi |
|
|
|
|
со
(1.123)
:=0
Стремление снизить инерционность контактного первичного преобразователя приводит к уменьшению его геометрической формы. Кроме того, если в качестве такого преобразователя при меняется голая термопара, то изменяются также и соотношения ее основных размеров. В этом случае термопары можно по форме уподобить усеченному конусу. Поскольку голые миниатюрные термопары широко используются, приведем для них основные соотношения, связывающие контролируемые температуру и теп ловой поток с величиной вторичного сигнала.
48
Итак, имеем усеченный конус с начальным сечением г = го при х = 0 и углом раствора р. Боковая поверхность теплоизоли-
рована, т. е. выполняется условие О |
< 5. 10"a . В цилин- |
к
дрической системе координат термоприемник ограничен поверх ностями
х = 0; г = r„ - j - х tg Р; 0 ==S ф ==S 2я. Таким образом, для всех точек, принадлежащих телу,
|
0«Sx=sSoo; 0 < |
г < r 0 + |
x t g P < о о . |
|
Упростим уравнение |
(1.3) |
с учетом |
рассматриваемой задачи. |
|
Учитывая |
изотропность |
тела |
и то, что термопара является телом |
|
вращения, |
получим |
0L — о- 2—•0•. |
||
|
|
|||
|
|
Эф |
' Эф |
|
В результате уравнение теплопроводности (1.3) можно значи
тельно |
упростить до |
вида |
|
|
|
|
|
|
|
(дЧ |
1 |
dt |
, |
дЧ |
|
(1.124) |
|
|
|
|
г |
дг |
' |
дх2) |
dt |
|
|
|
|
|
|||||
Введем |
|
|
|
|
|
|
|
(1.125) |
5 |
(х) |
= л (г |
S (ж) |
|
(1.126) |
|||
0 + |
х tg Р)2 . |
|||||||
осредненную |
по |
сечению |
х |
температуру |
|
где
После интегрирования уравнения (1.124) по сечению х получим следующее выражение:
|
|
|
dt \ |
п х _ n dt |
|
|
2 « ( r „ - f * t g P ) ( f | z - 2 t g p f | J - |
|
|||||
- 2 я |
tg2 p^, |
, |
a2 (so |
|
(1.127) |
|
|
dx2 |
|
64 ' |
|||
|
|
|
|
|
||
где s — боковая поверхность |
преобразователя. |
|
||||
Уравнение (1.127) после некоторых преобразований |
приводится |
|||||
к виду |
|
|
|
|
|
|
dt |
(дЧ^ |
|
2\i dt |
•) (Bi0 tgP)2 |
(1.128) |
|
dFo |
\ dW% |
|
dW |
|
|
|
где |
W = B i 0 + |
BitgP; |
|
|
||
|
|
|
||||
|
|
|
|
ax |
i — t g 2 |
|
|
|
|
|
; м- |
l — L t g » p |
4 П . А. Коротков |
49 |