Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Коротков П.А. Динамические контактные измерения тепловых величин

.pdf
Скачиваний:
8
Добавлен:
24.10.2023
Размер:
9.4 Mб
Скачать

Из сравнения результатов численных расчетов, выполненных по выражениям (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; Вх =

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 ) 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

Соседние файлы в папке книги из ГПНТБ