![](/user_photo/_userpic.png)
книги из ГПНТБ / Коротков П.А. Динамические контактные измерения тепловых величин
.pdfИз сравнения результатов численных расчетов, выполненных по выражениям (1.59) и (1.60), следует, в частности, что п р и / <
< 0 , 1 можно в уравнении (1.6) пренебречь значением Ъ, т . е . положить его равным нулю. Величину погрешности измерения температуры при b =j= 0 можно определить по кривым, приведен ным на рис. 3, где d = —f~-
Погрешность измерения температуры t из-за переходного процесса, имеющего место в первичном преобразователе, можно
d |
|
|
|
|
|
1,0 |
|
|
|
|
|
0,8 |
|
|
|
|
|
0,6 |
|
|
|
|
|
0,4 |
|
|
|
|
|
0,2 |
|
|
|
|
|
О |
1,0 |
2.0 |
3.0 |
4.0 |
5.011% |
Рис. |
3. Зависимость |
величины |
оттока тепла d от первичного |
||
|
преобразователя |
при т ^ > — .' |
|
||
с |
учетом условий |
(1.13)—(1.15); |
2 - е |
учетом условий |
(1.16)—(1.18) |
найти как расчетным путем, используя при этом соотношения (1.57) — (1-60), так и графики, приведенные на рис. 4 и 5. В по следнем случае для этого достаточно знать Fo и Л. Максимально возможная величина погрешности будет составлять К (Fo, А). Если Fo > 2,4, влиянием переходного процесса на формирование измерительного сигнала можно пренебречь.
Переходную характеристику первичного преобразователя можно определить как расчетным путем, используя для этого, например, приведенные выше соотношения, так и эксперименталь ным. В последнем случае для этого достаточно тем или иным спо собом поместить преобразователь в пламя газовой горелки и от метить при этом зависимость изменения регистрируемого сигнала во времени, которое как раз и будет графическим изображением переходного режима исследуемого термоприемника.
Аналитическая форма переходного процесса в этом случае может быть установлена, например, путем аппроксимации экс периментального графика Кэ (Fo) формулами (1.57) — (1.60). При этом следует иметь в виду, что графики Кэ (Fo) не могут полностью исключить из рассмотрения теоретический анализ решений уравнения теплопроводности, описывающих процессы
20
В измерительной среде. Каждый из методов должен лишь рацио* нально дополнять другой.
Снижение инерционности контактных первичных преобразо вателей приводит к уменьшению геометрических размеров их теплоприемной части. Это, в свою очередь, увеличивает влияние оттока тепла по токопроводам на измеренное значение темпера
туры. |
|
Переходный процесс в таких |
термоприемниках описывается |
уравнением (1.6) с условиями |
(1.13)—(1.15). Для простоты |
|
K(Fo) |
|
|
0,2 б 4 0,6 0,8 |
1,4 Го |
Рис. |
4. Переходная характеристика мо |
Рис. 5. Переходная |
характеристика |
дели |
термоприемника, удовлетворяющей |
модели термоприемника, удовлетво |
|
уравнению (1.6) с условиями (1.13)—(1.15) |
ряющей уравнению (1.6) с усло |
||
|
|
виями (1.16)—(1.18) |
исследований |
положим |
А = |
0 и ZL |
Ф о о . Решение |
этого |
урав |
|||||
нения, |
приведенное к виду, |
удобному для исследования |
пере |
||||||||
ходного |
режима, |
|
имеет вид |
|
|
|
|
|
|||
|
|
|
1 + |
в |
+ 2 |
V . ~т |
ехр (— ml Fo) sin (m„l) |
(1.61) |
|||
|
|
|
|
л = 1 |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
где |
An |
|
|
|
в |
В |
|
ZiX |
m„l |
|
|
sin 2mnl |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
||||
определяется |
2mnl |
|
|
|
|
|
= —Bm„/ или из формулы |
||||
из уравнения tg (tnj) |
|||||||||||
mnl |
= (n — 0,5 + |
|
P) я . В табл. |
1 приведены значения параметра р |
|||||||
при |
mnl |
= 1 -=-10 и Z/, равных 0,5; 1; 2. Используя |
лишь |
первый |
член ряда (1.61), получим К (Fo) с учетом условий (1.16) — (1.18).
Эти условия |
выполняются при |
Fo С 1- |
Если площадь воспринимающей тепло части первичного пре |
||
образователя |
S t соизмерима с |
площадью поперечного сечения |
21
токопроводов S2 , то переходный процесс в таких термопарах описывается выражением
|
K(Fo) = |
|
ехр х2 |
т |
sin %xdx |
A i ( l - A i ) |
|
|
|
|
|
|
|||
|
|
|
Ai(Ai + l) |
|
|
||
где |
|
|
Хер . |
|
|
I |
|
|
|
h,= |
|
1- |
|
||
|
|
ViPi ' |
|
V'a |
|
||
|
|
|
|
|
|
||
Индекс |
1 относится к токопроводам. |
|
|
||||
При |
I —> о о |
|
|
|
|
|
|
|
|
К (Fo) = erf |
[2 VFoj ' |
|
|||
|
|
|
|
|
|
||
Если h1 |
< 1, то в (1.62a) Fo = |
|
Если Лх > 1, то Fo |
(1.62)
(I.62a)
ax If
Значение параметра hlt входящего в выражение (1.62), зависит
Т а б л и ц а ! 0 Т т е п л о п Р о в ° Д н ы х свойств то-
Зависимость между параметрами Zt, Р и тп1
|
Значения |
р при L, , |
равных |
Ti |
0,5 |
1 |
2 |
|
|||
1 |
0,230 |
0,145 |
0,084 |
2 |
0,119 |
0,065 |
0,033 |
3 |
0,078 |
0,040 |
0,020 |
4 |
0,055 |
0,029 |
0,015 |
5 |
0,045 |
0,021 |
0,010 |
6 |
0,036 |
0,018 |
0,009 |
7 |
0,028 |
0,015 |
0,007 |
8 |
0,018 |
0,012 |
0,004 |
9 |
0,012 |
0,008 |
0,003 |
10 |
0,008 |
0,004 |
0,001 |
|
|
|
6 = |
Pi = - г ;
копроводов.
При технических расчетах наряду сточными выражениями, описывающими переходный про цесс, используются также при ближенные формулы:
к, ff„(6) = fli l - e x p ( - 6 ) £ J £ -
|
|
л = 1 |
|
|
|
|
|
(1.63) |
|
/UP ) = * i |
г- |
р2 |
X |
|
P , - P l |
||||
|
|
|
||
X ехр(— р\т) |
|
|||
+ Х ^ ? Г е х р ( " Р 2 Т ) ] ' |
( L 6 4 ) |
|||
где |
|
|
|
|
т; |
|
|
|
I
22
Данные выражения описывают переходный процесс & цепных схемах.
Если 0 ^ - ^ - ^ 5 - Ю - 2 , то |
|
А |
|
*„(б1) = 1 - е х р ( - 6 ) 2 - н г . |
(1..64а| |
п = 0 |
|
где б = — т.
с
Используя соотношения (1.63) |
— (1.64а), получим зависимости |
|
Кп Фд П Р И различных значениях |
параметра klt |
которые приве |
дены на рис. 6. |
(1.57) — (1.62) |
|
Замена точных выражений |
приближенными |
(1.63) и (1.64) означает замену распределенных параметров термо приемника сосредоточенными. Отсюда следует, что она прибли женно эквивалентна лишь при рассмотрении переходного про цесса и неприменима при изучении распределения тепловых волн вдоль преобразователя. Поэтому указанная замена справедлива лишь по отношению ко входу и выходу.
Выражения (1.63) и (1.64) не учитывают обратного влияния звеньев. Поэтому реализация цепных схем с учетом приведенных выражений возможна лишь при развязке их элементами с одно сторонней проводимостью или при таком подборе элементов, чтобы обратное влияние было менее допустимого. В последнем
случае выполняется условие п - 1 > — - .
Если по каким-либо причинам не представляется возможным развязать звенья между собой, переходный процесс в них следует рассчитывать с учетом обратного влияния. При этом выражения получаются весьма громоздкими. Так, например, при п = 2
К» = (бт) = 1 - |
+ - Ь . ) ехр ( - -Ji- т ) , |
(1.65) |
где |
|
|
т1 = |
-!Ь. |
|
{аах)
Т2
*а ' •*
Тепловое воздействие среды может быть определено также зада нием коэффициента теплообмена и температуры среды либо зада нием теплового потока,
23
Из (1.21) и (1.22) при Zz — о о удельный тепловой поток |
связан |
с контролируемой температурой операторным выражением |
|
t (р) = q (р) [Z, (р) ch [у (р) l]-\-Zc (р) sh [у (р) 1}}. |
(1.66) |
Если контролируемый процесс весьма кратковременный, когда
выполняется |
условие т„ « |
, то |
|
|
t(p)=~q(p)Zc{p). |
(1.67) |
|
Из (1.67) с |
учетом А = 0 к |
~ < 10"3 |
|
|
|
1 |
|
|
~t(p) = |
(Xcppf~-q(p). |
(1.68) |
Соотношение (1.68) обычно используется при применении пленоч ных первичных преобразователей. Эти соотношения охватывают большинство случаев использования контактных первичных пре образователей температуры с учетом уравнения (1.6).
Приведем примеры, иллюстрирующие применение соотношений (1.56) — (1.58).
П р и м е р |
1. |
Пусть |
первичным |
преобразователем выбрана |
термопара, |
||
выполненная в основном из меди |
(а = |
0,36 |
м2 /ч). Форма преобразователя темпе |
||||
ратуры — цилиндр. |
Геометрические |
параметры термоприемника |
следующие: |
||||
/ = Ю - 3 м, Sx = |
3-10"6 |
м, 5 2 = |
10~8 м. |
Термопара заделана в |
эбонитовую |
втулку. Длительность опыта составляет 10"2 с. Регистрирующий прибор имеет полезную длину шкалы 102 мм. Разрешающая способность прибора составляет
10"3 |
мм. Требуется выбрать формулу для расчета переходного процесса в первич |
|||||||||||
ном |
преобразователе. |
|
|
|
|
|
|
|
|
|||
|
|
Так как - ^ J - |
> 1, то в (1.21) можно положить Z/ = |
оо. |
|
|
|
|
||||
|
|
Значение параметра b можно положить равным нулю, так как |
к |
< |
1. Зна- |
|||||||
|
|
Л |
||||||||||
чение параметра |
Fo будет: |
10"3 |
|
|
|
|
|
|||||
|
|
|
|
|
|
|||||||
|
|
Из полученных данных следует, что переходный процесс может быть рассчи |
||||||||||
тан |
по |
формуле |
(1.23). |
|
|
|
|
|
|
|
|
|
|
|
П р и м е р 2. Требуется |
выбрать значение / таким, |
чтобы переходный про |
||||||||
цесс |
в первичном преобразователе подчинялся (1.60); |
при этом |
т а |
= |
10~3 |
с, |
||||||
а = |
0,36 |
м2 /ч. |
|
|
|
при Fo <g 1, Zi = оо, |
А = 0, |
|
||||
|
|
Так |
как выражение |
(1.60) справедливо |
то |
|||||||
I |
> |
Yaxa- |
Подставляя |
в данное выражение |
исходные |
данные, |
получаем |
/ = |
||||
= |
3 - Ю" 3 м. |
|
|
|
|
|
|
|
|
|
||
|
|
П р и м е р |
3. Переходный процесс в первичном |
преобразователе |
подчи |
няется выражению (1.60). Требуется определить значения п и рт, при которых
точную модель |
теплопередачи |
можно заменить приближенной, |
учитывая, что |
|||||
погрешность нахождения |
контролируемой температуры |
не должна |
превысить |
|||||
С = Ю°/о. |
|
|
|
|
|
|
|
|
Исходные данные: ch А = |
1, £ г = |
10"4: Fo = 0,8. |
Из рис. 4 находим, что |
|||||
К (Fo) = |
0,85. |
Разбиваем |
интервал Fo = |
0-^-0,8 на три равные части; |
при этом |
|||
К (0,8) = |
0,85; |
К (0,5) = |
0,45; |
К (0,2) = |
0,27. Из рис. 6 имеем, что с учетом за |
|||
данной погрешности п= |
1,0 |
и р т = 0 - ь 0 , 8 . |
|
|
|
|||
П р и м е р 4. Известно, что процесс теплопроводности в первичном преобра |
||||||||
зователе |
подчиняется точному |
выражению (1.58). При этом Fo = |
0,3; Вх = 5Х |
24
XlO 2 ~-'- = |
0,75. Требуется определить значения 6 и п, |
при которых точную |
А |
|
|
модель теплопередачи можно заменить приближенной Кп |
(бт), чтобы при этом |
|
погрешность |
определения t (т) не превысила, например, величины 10%. |
Определение значений параметров аппроксимирующего выражения Кп (бт) производится обычно по трем значениям параметра Fo, например при К (5-10"2) ±
± K ( 1 0 " 2 ) ; К ( Ю - 1 ) ± К ( Ю - 2 ) ; К (2-10~2) ± Л' (10~2).
Используя графические зависимости, приведенные на рис. 6, приходим к выводу, что при указанных величинах Fo точное вы ражение К (Fo) может быть заменено приближенным при п = 3. *
В тех случаях, когда значения теплофизических параметров термоприемника неизвестны, контролируемое значение темпера туры t может быть найдено
посредством |
трех |
измерении |
к„(6,) |
|
||||||
температуры |
соответственно |
|
|
|||||||
в моменты времени т, |
т + А т |
|
|
|||||||
и т -|-2Дт |
|
[118]. |
Если |
при |
|
|
||||
этом значение |
истинной |
тем |
|
|
||||||
пературы |
изменяется во вре |
|
|
|||||||
мени по |
экспоненциальному |
|
|
|||||||
закону, |
то |
|
|
|
|
|
|
|
|
|
t. |
2t2 |
- |
(h |
+ |
*,) |
(1.69) |
|
|
||
|
|
|
|
|
||||||
где tu |
t2 |
и |
|
t. |
• значения из- |
Рис. 6. Переходные характеристики мо |
||||
меренных |
|
температур |
|
соот |
дели термоприемников, удовлетворяющей |
|||||
ветственно |
в моменты |
време |
|
выражению (1.63) |
||||||
ни т, т + Ат и т |
2Ат. |
|
|
|||||||
В [118] показано, что погрешность, которую может иметь |
||||||||||
найденная |
|
контролируемая температура, |
вычисляется по формуле |
|||||||
|
Ми |
|
|
|
|
4 К т ( А ) 2 - 2 ( К т |
+ 1 ) т - |
|||
t |
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
А |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
|
|
|
(1.70) |
где Кх |
|
|
|
|
|
|
|
|
|
|
Как следует из [118], точность определения температуры глав |
||||||||||
ным образом зависит от значений —- и |
т. е. от длины участка |
|||||||||
кривой |
по времени между tx |
и t3, и от места расположения этого |
участка на экспоненте. Данный метод был опробован при аэроди намических испытаниях в газодинамической установке кратко временного действия. Погрешность результатов оценивается в ±20—40% . При анализе переходных процессов в измерительной среде можно использовать как переходные характеристики, так и ее частотные свойства. Выбор вида характеристик в основном определяется требованием простоты наблюдения контролируемых
25
•особенностей процесса. Временные свойства первичных преобразо вателей подробно изложены выше. Остановимся на их спектраль ных свойствах.
Для получения спектральных характеристик первичных пре образователей достаточно применить к временным характеристи кам преобразование Фурье (1.38)—(1.41). Так, применяя это
Рис. 7. |
Амплитудно-частотная характе- |
Рис. 8. Фазово-частотная характе |
ристика |
первичного преобразователя, |
ристика первичного преобразова- |
|
удовлетворяющего (1.71) |
теля, удовлетворяющего (1.71) |
преобразование к (1.58) и (1.60), получим следующие амплитудно- фазово-частотные характеристики измерительного тракта:
|
C(M |
= y~Bl + |
Bl exp [ - / a r c t g . ( - ^ _ ) ] |
; (1.71) |
|||
|
С(/ш) = |
( В ? + в И ~ е х р ( - / a r c t g | i - ) ; |
( I J 2 > |
||||
где |
S j «= COS D ch D; |
B2 |
= |
sin D sh D; D = |
I ~ \ [ • |
||
Для |
простоты |
положим, |
что - ^ д - < 5-10"2 |
и — ^ - < 5 - 1 0 " 2 . |
|||
На рис. 7—10 даны амплитудно-частотные и фазово-частотные |
|||||||
характеристики |
первичных |
преобразователей. На рис. 7 и 9 |
|||||
отмечены частоты среза сос. |
|
|
|
||||
По этим характеристикам можно достаточно просто |
установить |
степень искажения контролируемой температуры в первичном преобразователе. Для этого достаточно сравнить диапазон суще ственных частот контролируемого процесса с полосой прозрач ности первичного преобразователя. Если окажется, что диапазон существенных частот контролируемого процесса расположен в по
лосе |
прозрачности |
(пропускания) первичного преобразователя, |
то он |
не исказит |
исследуемый сигнал. |
-26
Если же область существенных частот контролируемого про цесса значительно перекрывает полосу прозрачности первичного преобразователя, то контролируемое и измеренное значения тем ператур окажутся различными. Для восстановления контроли руемого процесса следует воспользоваться выражениями (1.32) — (1.36).
Если контролируемое значение температуры изменяется во времени по синусоидальному закону и при этом процесс тепло
го-;
А
0,8
0,6
0,4
0,2
О—10,4 0,8 1,2 1,6 2.0
Рис. 9. Амплитудно-частотная |
харак |
Рис. |
10. Фазово-частотная характери |
||
теристика первичного преобразователя, |
стика |
первичного |
преобразователя, |
||
удовлетворяющего |
уравнению |
(1.6) |
удовлетворяющего |
уравнению (1.6) |
с условиями (1.16)—(1.18) |
с |
условиями |
(1.16)—(1.18) |
|
передачи в |
измерительной |
среде подчиняется |
уравнению (1.6) |
|
с условиями |
(1.13) — (1.15), то при qx |
= 0 |
|
Си (/©) = ехр (/со) -
Хехр (-
л3а |
п3 sin (лп) |
X |
|
я 4 и 2 а 2 |
|
|
16/г -f- |
|
я2 /г2 |
\ |
(1.73) |
|
|
Из (1.73) следует, что преобразователь температуры вызывает за паздывание. Величина запаздывания передачи тепловой энергии
определяется значением параметра |
При частотах контро |
||
лируемого |
процесса, для которых |
удовлетворяется |
условие |
4 1 2 ш > 1, |
можно приближенно считать, что tK (т) = |
t (т). |
Амплитудно-частотная и фазово-частотная характеристики цеп ных схем замещения описываются выражениями:
С (со) |
(1.74) |
(/If+1)3
27
ф(<о) = е х р ( — / a r c t g - J - ) , |
(1.75) |
где А г = - ^ - .
Формулы (1.74) и (1.75) справедливы при условии, что отдельные звенья разделены между собой элементами с односторонней про-
I2 Р
водимостью, и приближенно справедливы, если — ^ 10—^—5-
[22]. Из (1.75) следует, что для частот контролируемого процесса, для которых выполняется условие Ах > 1, первичный преобразо ватель можно считать прозрачным.
Комплексный коэффициент передачи двухъемкостной схемы за мещения первичного преобразователя с учетом обратного влияния второго звена на процессы в первом звене имеет вид:
|
|
|
|
i_ |
|
1С (/СО) 1 = ^ - 1 ) 2 |
+ ( |
7 |
» |
^ . |
( 1 7 6 ) |
[ ( 1 - 7 > 2 ) 2 |
+ ( 7 > ) 2 ] 2 |
|
|||
ф И = [ _ / a r c t g ( T |
^ - |
- |
^ |
) ] . |
(1.77) |
Рассмотрим примеры, иллюстрирующие использование ча стотных характеристик первичных преобразователей температуры.
П р и м е р |
1. Требуется |
определить, соответствует |
ли значение |
измеренной |
температуры контролируемой |
при следующих исходных |
данных: а = |
0,36 м2 /с; |
|
/ = 10~4 м; / в = |
35 1/с; / н = |
0. Амплитудно-частотная |
характеристика пер |
вичного преобразователя подчиняется зависимости (1.72). Длительность иссле дуемого процесса составляет 10"3 с, при этом длительность переднего фронта 10"4 с. При этих данных имеем, что верхняя граничная частота области сущест венных частот контролируемого процесса составляет 50 1/с. В качестве нижней граничной частоты можно принять нулевую частоту. Из рис. 9 следует, что Д/ =
= |
50 1/с. Таким образом, первичный преобразователь |
по |
отношению к контро |
|||
лируемому сигналу можно считать прозрачным и тогда |
/ и |
== г. |
|
|||
|
П р и м е р |
2. Определим значения |
параметров |
первичного |
преобразова |
|
теля 1м а, при которых с учетом выражения (1.72) ги = |
t. При этом <мв = 100 1/с |
|||||
и |
сон = 0. |
|
|
|
|
является медь |
|
Материалом для изготовления первичного преобразователя |
|||||
(а = 0,36 м2 /ч). |
Тогда |
|
|
|
|
|
|
|
/ ^ / о й ^ = |
6-10"3 м. |
|
|
|
4 . О П Р Е Д Е Л Е Н И Е Д И Н А М И Ч Е С К О Г О К О Н В Е К Т И В Н О - К О Н Д У К Т И В Н О Г О Т Е П Л О В О Г О П О Т О К А
Тепловой поток, в отличие от температуры, непосредственно измерен быть не может. Для его определения необходимо иметь информацию о физико-геометрических параметрах измеритель ного тракта и изменении температуры свободной поверхности термоприемника, когда х = 0 (см. рис. 1 и 2), во времени. Ди намический тепловой поток для любого момента времени и лю-
28
бой точки в соответствии с уравнением Фурье может быть пред ставлен выражением
д(х, х) = - Х д - ^ ^ . |
(1.78) |
||
Однако в основном представляет |
инетерес поток тепла, |
поступаю |
|
щий на свободную поверхность |
термоприемника: |
|
|
д(т) = -к*<£*}, |
х = 0, |
(1.79) |
Согласно закону непрерывности теплового потока, в линейном
приближении выполняется следующее |
граничное условие: |
|
а [ ( с ( т ) - / ( т ) ] = |
- |
^ , |
х = 0, т ^ О . |
(1.80) |
Отсюда следует, что тепловой поток находится как разность тем ператур tc — t, умноженная на величину коэффициента тепло передачи. Если значение t может быть определено посредством измерений и расчетов (см. п. 3), то величины tc и а в общем слу чае оказываются неизвестными. Поэтому вычисление теплового
потока |
производится главным |
образом |
по формулам |
(1.20) и |
|
(1.22) с учетом условия (1.79). |
|
|
|
|
|
Приведем расчетные формулы |
для |
вычисления |
теплового |
||
потока |
с учетом рассмотренных |
в |
п. 3 |
условий опытов. Выбор |
той или иной формулы определяется условиями измерений и зна чениями физико-геометрических параметров первичного преобра зователя.
Если процесс теплопередачи в измерительной среде подчи
няется уравнению |
(1.6) и условиям |
(1.7) — (1-9), то при |
< 1 |
||||
|
|
|
|
т |
|
|
|
q (т) = yj&R tg (A) t (т) + |
У |
>f- |
j |
exp ( - т т х ) X |
|
||
|
|
|
|
о |
|
|
|
Z4{ |
M l ( Л + 0,5)* Я» + |
Л * \ е |
Х р 1 |
F o 7 |
|
|
|
п=0 |
|
|
|
|
|
|
|
— ехр |
|
( т - |
|
|
|
||
|
|
Кт7 |
|
|
|||
|
|
|
|
|
|
||
Из (1.81) следует, |
что при т > - д ^ - |
|
|
|
|
|
|
|
q(x) = Vtepm |
th (A) t (т). |
|
|
(1.82) |
||
РЯД в выражении (1.81) быстросходящийся. |
Кроме |
того, с увели |
|||||
чением п значения членов ряда резко уменьшаются. |
Поэтому для |
29