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

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

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

практических целей обычно ограничиваются первыми тремя чле­ нами ряда (1.81).

В тех случаях, когда исследуемый процесс описывается уравне­ нием (1.6) с условиями (1.10) — (1-12) и при этом -^=- 5» 5. 10~2 , значение теплового потока рассчитывается по формуле

ехр (— тхх) X

 

X Vm

Ф (mtj) t' (х •— хх)

dxx.

(1.83)

Kq(la)

 

 

 

 

 

 

 

 

Чует

 

 

 

 

 

 

 

 

11

 

 

 

 

 

 

 

 

1 ,2

 

 

 

 

S

/

 

 

 

 

 

 

 

3

 

 

5

 

 

 

 

 

 

4

 

8

12 16 20

24

28

30 гК

 

Рис.

11. Зависимость

Ко

 

от хк:

 

 

Ко¥уСТ

 

 

 

 

 

 

 

 

кривые

/, 2,

3,

4,

5 построены при

различных

-значе­

ниях параметра

Bi , соответственно равных

1, 25, 50, 250,со

Положив в (1.83) т = 0, получаем простое выражение для расчета величины теплового потока

(1.84)

Выражение (1.84) может быть применено, когда тх > 2. Рассмотрим случай, когда процесс теплопередачи в зоне из­

мерения описывается уравнением (1.6) с условиями (1.7) (1.9), при этом -ri— :=г5 . 10" 2 и Л = 0. Величина теплового потока рассчитывается по формуле

30

где

 

 

 

об

 

 

 

 

 

 

 

K„{r)= S

Л е х р ( — п 2 Fo).

 

 

 

 

 

 

 

n=l

 

 

 

 

 

На практике значение п обычно не превышает трех.

 

4

 

Кривые

К

(х)

даны

на рис. 11, где тк =

/2

=

- р

5 - ^ -

 

—-^-1п-^-

 

А <?уст

 

 

а

п

6

 

 

 

/2

 

 

 

 

 

 

= 0,02915. \0~6с. Очевидно, что при ха > 30тк тепловой поток достигает стационарного значения. Кроме того, с увеличением значения параметра Вх скорость переходного процесса возрастает. Это вызвано шунтирующим действием тепловой емкости тепло­ вым сопротивлением Zt. При т > ха

^ ) = ^ ) ( в Т П Л -

(1.86)

Если В х > 1, то

Погрешность измерения q из-за того, что Zt ф оо, рассчиты­ вается по формуле

KlyCT

(1.88)

 

Zl = C°

Значение безразмерного комплекса Вх для сварных голых тер­ мопар колеблется в диапазоне 50—250. Отсюда следует, что по­ грешность определения теплового потока без учета теплоотвода может достигать 50%. При прочих равных условиях, погрешность измерения теплового потока без учета оттока тепла по державкам

т

 

за счет теплопроводности зависит также от величины

.

Если в выражении (1.85) ограничиться лишь первым членом

ряда, то получим выражение для определения теплового

потока

сучетом условий (1.10) — (1-12).

Втех случаях, когда процесс теплопроводности в измеритель­ ной среде моделируется одноемкостной схемой замещения, ве­

личина теплового потока рассчитывается по формуле

ч ( t )

= Ш

(

J I 1 - е х р [ - J

^

M

~

]

^

d

x

-

 

 

 

о

 

 

 

 

 

 

 

 

 

 

Я -4- Я

 

 

 

 

 

 

 

 

Если

т >

^ 2

" , то из (1.89) следует,

что

 

 

 

 

 

 

31

Если в (1.89) положить Хи --= 0, то

 

q (т) = J ехр (— Ро) ? (т — тх ) dxv

(1-91)

При моделировании измерительной среды двухъемкостной схемой замещения с учетом обратного влияния второго звена на первое величина теплового потока рассчитывается по формуле

<?(т)

К-Кг

 

ехр (— K&x) Г (т — тх ) d%1

 

1

 

 

 

 

J

^ ( - - ^ t ' i r - x ^ d r ,

(1.92)

где

К

7*1 =

;

т2

 

/2

 

 

1\

_

/ 2а

+

2 —

1

 

 

Тг =

аах '

 

 

'

ai

7\;

 

 

а

Г

 

All '

 

 

 

 

 

 

^12 = <aPl

Из (1.92) следует, что введение второго звена в измерительную цепь приводит к увеличению теплового потока.

Если измерительная цепь состоит из двух идентичных зве­ ньев, то

 

/тп

 

/ 4

гр

 

 

 

гр ( гр

 

 

*•

о

1

 

 

Ji~~~lW'

У 2 ~ " % ;

 

'з—'ъ

 

 

 

7\ = 2

а,

 

 

а

 

 

 

— =

2 — .

При хЯ и ^

„ _

2

я

>

5. Ю

- 2

значения коэффициентов,

>

5-10"

и —-

 

входящих в выражение (1.92), могут быть рассчитаны по форму лам:

7\>

=

/2

ср

, J i £ i P i _ .

 

яГ

Ях +

Я щ '

 

 

 

3 — 7 41

 

7 4 — 1я I+1ян I 1

я + яи

32

5. О П Р Е Д Е Л Е Н И Е З Н А Ч Е Н И И Ф И З И К О - Г Е О М Е Т Р И Ч Е С К И Х П А Р А М Е Т Р О В П Р И Е М Н Ы Х П Р Е О Б Р А З О В А Т Е Л Е Й

В З А В И С И М О С Т И О Т У С Л О В И Й И З М Е Р Е Н И Й

Измерения динамических температур и тепловых потоков кон­ тактными преобразователями сопровождаются появлением воз­ мущения в исследуемой зоне. Это вызвано сложностью взаимо­ действия приемного (первичного) преобразователя со средой. Ин­ тенсивность возмущений определяется соотношением значений физико-геометрических параметров первичного преобразователя и среды.

Вследствие наличия возмущений в зоне расположения первич­ ного преобразователя измеренные значения температуры и теп­ лового потока отличаются от контролируемых. Контролируемые значения температуры и теплового потока могут быть определены двумя способами: введением необходимых поправок в результаты измерений и соответствующим выбором физико-геометрических параметров первичного преобразователя в зависимости от условий опыта, при которых можно установить равенство между измерен­ ными и контролируемыми величинами. Предпочтение следует отдать второму способу. Теоретические исследования, направлен­ ные на изучение влияния первичного преобразователя на исследуе­ мый процесс, весьма малочисленны [52, 53, 114].

Влияние первичного преобразователя на процесс теплопередачи в зоне измерения было обнаружено в [116, 117].

Значения физико-геометрических параметров первичного пре­ образователя выбираются исходя из результатов исследований зависимости вида

- f L o 3 ( B i , F o ) =

,

(1.93)

 

(Р)

 

характеризующей скорость изменения потока тепла в измеритель­

ной среде

во времени.

 

 

 

 

 

Рассмотрим зависимость вида (1.93) на примере решения урав­

нения

(1.6) с условиями

(1.10) — (1.12). Для

простоты положим,

что Zl

= о о и А = 0. В

результате получим

следующее

выра­

жение:

 

 

 

 

 

 

 

 

 

-j- Lo3 (Bi, Fo) = Bi exp (Bi2 Fo) erf (Bi i

Fo).

(1.94)

Из (1.94) следует, что

при

Fo —> о4^-—> 1.

С увеличением

значения

Fo величина

Lo3

монотонно уменьшается. В

пределе

при Fo —» оо Lo 3 —* 0. Следовательно, при малых значениях Fo, когда - j - Lo3 (Bi, Fo) л* Bi, тепловой поток зависит в основном от величины конвективного сопротивления и температуры среды.

Если, в частности, tc (т) = const и — Lo 3 (Bi, Fo)

Bi, то

3 П. А. Коротков

33

величина теплового потока практически не зависит от времени.

При

-1 Lo 3 (Bi,

Fo)

Bi будут проявляться особенности

пер­

вой

характерной

зоны

зависимости Lo 3 (Bi, Fo) = ср (Fo).

Для

данной зоны изменения Lo 3 (Fo) характерно то, что тепловой по­ ток при заданной температуре среды зависит в основном от зна­

чения

конвективного сопротивления.

При

этом

t < tc.

 

С увеличением Fo значение Lo 3 (Fo) уменьшается и увеличи­

вается

роль физических

свойств

первичного

преобразователя

в организации теплового

потока.

При

Lo3 (Bi, Fo)

-

 

lc3

 

 

 

 

 

 

 

 

fa'

Bi-5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

!Dl

 

 

 

 

 

 

 

 

W"

 

i

 

 

 

 

 

 

fr­

 

Г

\

 

 

 

 

 

 

 

 

 

 

it)'3

 

\

 

 

 

 

W

 

 

 

 

 

 

 

 

 

 

to1

 

w3

 

Fo

 

 

 

 

 

 

 

 

Рис. 12. Изменение величины входной проводи­ мости измерительной схемы во времени при раз­ личных значениях параметра Bi

величина теплового потока в равной мере зависит от свойств среды и первичного преобразователя. Диапазон изменения Fo, при ко-

 

 

 

В i X

 

 

 

 

тором Lo 3

(Bi, Fo)

- g - - у . будем называть второй

характерной

областью изменения зависимости Lo 3 (Fo).

 

 

При

дальнейшем

увеличении

Fo

величина

Lo3

(Fo) устрем-

 

 

 

 

 

 

Bi У

ляется

к

нулю. Приближенно

при

Lo 3 (Fo)

=^ ~yo~~f" ^ ~ ^ и

величина теплового потока определяется в основном теплофизи-

ческнми

свойствами первичного

преобразователя. Диапазон из-

 

 

 

 

X

менения

Fo, при котором 0

Lo 3

(Fo)

10"1 Bi -у-, будем назы­

вать третьей характерной зоной зависимости Lo 3 (Fo). Этой зоне свойственны особенности граничных условий первого рода. Отсюда следует, что при тех значениях Fo, когда измерительный процесс подчиняется особенностям второй характерной зоны, происходит переход граничных условий третьего рода в граничные условия первого рода.

На рис. 12 дана зависимость Lo 3 (Bi, Fo) при различных зна­ чениях параметра Bi . Приведенные кривые позволяют определить необходимые значения а, I и Bi, при которых измерительный

34

процесс будет находиться в строгом соответствии с условиями опыта. Так, например, если требуется измерить тепловой поток, обусловленный конвекцией, то значения параметров а, I и Bi должны быть такими, при которых Lo 3 = (0,95н-1) Bi. Для со­ блюдения равенства t = tc необходимо, чтобы Lo 3 = Bi (0-f- --5-10-2 ).

Проследить особенности контролируемого процесса по гра­ фику Lo8 (Bi, Fo), Bi = const при малых значениях Fo и B i <t 1 затруднительно из-за монотонности его изменения во времени.

81-30

15 \

— ' —

2 ^ -

05

0.2

Tl dUj_

\dFo

Рис. 13. Скорость изменения величины Lo3 во времени при [значениях параметра Bi от 30 до 1

Определить наступление того или иного состояния контролируе­ мого процесса в этом случае удобнее посредством зависимости

 

 

 

xl д Lo3 (Bi, Fo)

, .

(1.95)

 

 

 

 

 

тт

1 о

1

x/3Loo(Bi, Fo)

На рис. 13 и 14 приведены графические зависимости

г ^ п —

 

 

 

 

 

Лого

при Bi =

const. При малых значениях

Fo эти зависимости разли­

чаются между

собой;

степень различия

зависит от значения кри­

терия Bi. С увеличением Fo рассматриваемые зависимости сбли­

жаются. Начиная с некоторого значения Fo

(в зависимости от

величины Bi)

 

 

т/_ _d_Lo3 (Bi, Fo) _ 1

(196)

dFo

V Fo "

 

 

 

Если величина д Loo зависит от Bi, то в измерительной среде про­

являются

особенности первой характерной зоны зависимости

3*

35

L o 3

(Bi, Fo). В этом

случае действуют

условия

(1.10) — (1.12).

п

тех значениях

с

Fo, когда

d Lo, (Bi,

Fo)

 

„ .

При

 

 

 

—' не зависит от Bi,

имеет место третья зона зависимости Lo 3 (Bi, Fo). В этом случае действуют условия (1.16) — (1-18) и t = tc.

Определить наступление того или иного состояния контроли­ руемой среды можно также путем «нормализации» функциональ­ ной зависимости Lo 3 (Fo), подобно тому, как это делается в тео­ рии электрических цепей с учетом особенностей «работы» систем

0

fO'

Ю2

I0J 104 I05 Ю6 f0' 10s W3 10ю Fo

to

Bi-1

/31

r'•2W'

 

 

•re-'

 

 

 

U 8Lot -e ~l dFo'u

Рис. 14. Скорость изменения величины Lo3 во времени при значе­ ниях параметра Bi от 4-Ю"1 до 10~3

с распределенными параметрами [30]. «Нормализация» здесь достигается изменением масштаба графика Lo 3 (Fo) на величину,

численно равную

отношению функциональных

зависимостей

Lo 3 (Fo) и Кп (Fo),

где Кп

(Fo) — коэффициент передачи

первич­

ного преобразователя по

тепловому потоку от

сечения

х = 0

д о х = - у . Выбранная форма коэффициента передачи обусловлена тем, что максимальное влияние термоприемника на процесс тепло­ передачи в измерительной среде имеет место при достижении ра­ венства q (0, т) = q(y~Y, т ) . С учетом условий (1.10) — (1.12) коэффициент передачи рассчитывается по выражению

ехр ( Bi 2 Fo + Bi I erf

B i K F o +

1

 

L o 3 ( B i ,

Fo)

4V~a

(1.97)

 

 

 

где X0 — длина свободного пробега электронов, участвующих в передаче тепла от сечения х = 0 до х = у .

36

Анализируя

приведенные

на

рис. 15

«нормализованные»

за­

висимости

LoH

(Bi,

Fo), Bi = const,

видим,

что при

малых

зна­

чениях

Fo, когда

LoH (Fo)

Bi, в измерительной среде выпол­

няются

условия,

характерные

для

первой

зоны

зависимости

Lo 3 (Fo).

При

дальнейшем

увеличении

Fo

значение

LoH

(Fo)

резко

возрастает,

достигая

максимальной

величины

L o H m a x -

При этом измерительный процесс подчиняется особенностям вто­

рой

характерной

зоны

зави­

 

 

 

 

 

 

симости Lo 3

(Fo).

После

 

до­

 

 

 

 

 

 

стижения максимального зна­

 

 

 

 

 

 

чения зависимость Lo3 (Fo)

 

 

 

 

/ 1

 

резко

устремляется к

нулю.

 

 

 

 

«

Крутизна

 

спада

зависит

от

 

 

 

 

/,5 Iff},

/

 

to'

 

 

 

1

величины

 

Bi.

При

LoH

>

 

1 — s

> Lo„ m a x

 

происходит

пере­

10'

 

 

 

10-

/ L

ход

граничных

условий

 

че­

 

 

 

твертого

рода

в

граничные

to'

 

 

 

 

 

условия первого рода. Таким

 

 

 

 

 

 

образом,

при

 

L o H ^ 5 - 1 0 ~ 2

 

 

 

 

 

 

в измерительной

среде

вы­

 

 

 

 

 

 

полняются

условия,

харак­

 

 

 

 

 

 

терные

для третьей зоны

 

за­

10'

 

 

 

 

 

висимости

Lo 3

(Fo).

 

 

 

 

 

 

 

 

 

 

 

10'

 

 

 

\

 

 

Решением уравнения

(1.6)

 

 

 

 

с

условиями

(1.7) — (1-9),

to'-

 

 

 

 

приведенным

к

виду

(1.93),

 

 

 

 

 

 

 

 

 

будет

выражение

 

 

 

 

10"

 

 

 

\\

\Fo

 

X Bi

L o 3 ( B i ,

Fo)

=

 

 

to'-

15.

Изменение

 

 

 

 

 

 

 

 

 

 

 

Рис.

«нормализованного»

 

 

 

 

 

 

 

 

 

 

 

значения

входной

проводимости изме­

 

£

 

H„sin(C„Bi0 )

 

 

 

рительной среды

во времени при различ­

 

n = l

 

 

 

 

 

 

 

 

 

ных значениях

параметра Bi

 

 

 

 

 

 

Вп cos (C„ Bi0 )] ехр [ -

(Bi0 Cnf

Fo],

(I.97a)

где

B i 0 =

 

- ~ - ,

r0—радиус

 

стержня

(первичного

преобразова­

теля).

 

 

ла

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Значения параметров Ап,

Вп и Сп

в выражении (1.97) рассчи­

тываются

 

по

формулам:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

tg(С„Bi0 )tg(CnW,)

=

\ + ^

{

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n-l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

к

 

 

 

 

 

 

 

 

Bn

=

-An]tg(CnBi0).

 

 

 

37

Эти формулы получены путем подстановки v =-

^- в (1.131),

(1.132)

и (1.103).

 

Из (1.131) с учетом приведенных выражений для коэффициен­

тов

Ап

и Вп

получаем:

 

 

X В

я = 1

 

 

 

 

 

 

 

2CnW1

+ 1 c o s ( c n J ^ ) - c o s ( c „ ^ ± ^ ) ] ^ x

 

 

 

со

я—1

х

е х р [ - (С„Bi0 )>Fo] + 22J

A~ ^ Z t c X ' " X

xexp [ - ( C „ B i 0 ) 2 Fo].

На рис. 16 приведены графические зависимости

 

Lo 3 (Bi, Fo),

Bi = const, полученные с помощью выражения

 

(1.97а). Используя

 

 

 

 

 

 

их,

можно

описанным

вы­

 

 

 

 

 

 

ше образом

выделить харак­

 

 

 

 

 

 

терные области

и

границы

 

 

 

 

 

 

измерительного

 

процесса.

0.8

 

 

 

 

 

Сравнивая

между

собой за­

 

 

\BHO'

 

 

висимости,

приведенные на

0,8

 

 

 

 

 

рис.

12

и 16, видим, что во

 

 

 

 

 

 

втором

случае

переходный

OA

 

 

 

 

 

процесс

протекает

быстрее.

0,2

 

 

 

 

 

Вследствие

этого

 

характер­

 

 

 

 

 

ные

зоны

здесь

смещаются

 

 

 

 

 

50 Го

в сторону

меньших

значений

f0'!210*5 Ю210' 0.5 10

5,0 10,0

Fo и наступление того или

Рис. 16. Изменение

величины входной

иного состояния

измеритель­

ной

среды

начинается

при

проводимости измерительной

среды во

времени

при различных

значениях

пара­

меньших

 

значениях

Fo.

 

 

метра

Bi

 

 

 

Итак,

 

первой

зоне

зави­

ряют

граничные условия

третьего

симости

Lo3

(Fo)

удовлетво­

рода,

а третьей зоне данной за­

висимости— граничные условия первого рода. Отсюда вторая

зона

является переходной

между

этими

условиями.

 

 

 

 

Исследуемые особенности измерительной среды можно уста­ новить также посредством ее амплитудно-частотных или фазовочастотных характеристик. Этими характеристиками обычно поль­ зуются в тех случаях, когда известна область существенных частот контролируемого процесса или когда в качестве регистрирующей аппаратуры применяется прибор типа анализатора спектра.

38

Так,

для случая (1.94)

имеем,

что

 

 

 

~ LO,(/<D) = Fo-o.5 {[1 + (Bi V2Ai)-l}>

+ ( B i / M ^ j - e . s

x

X exp / { - f -

arctg (Bi / 2 T J - 1

[ 1 + ( B i уГЩ)-^

.

(1.98)

По

выражению

(1.98)

построены

графические

зависимости

I

(/со) | и ф (со), которые даны

на

рис. 17 и 18. Характер­

У | Lo 3

ные зоны и их границы выделяются

следующим образом. По от-

111-

3-Ю.

W

81-10

5

3

1

!

з-ю' w •' з-ю7

10-Г 1—~—

Рис. 17. Частотная характеристика измерительной среды по тепловому потоку при различных значениях парамет­ ра Bi

ношению к частотам контролируемого процесса, для которых ве­

личина

[ Lo3 (/со) (

соизмерима с Bi,

выполняются

особенности

первой

характерной

зоны зависимости

(1.93). Если

4- Lo 3

(/to)!

Bi

 

А

1

 

 

 

— т о в измерительной среде проявляются осооенности, свой­ ственные второй зоне зависимости Lo 3 (Fo). По отношению к ча­

стотам контролЕфуемого

процесса, для которых

величина

I Lo 3 (/со) | близка

к нулю,

выполняются

особенности,

свойствен­

ные третьей зоне

зависимости (1.93).

 

 

Фазовые соотношения

измерительного

процесса также позво­

ляют определить степень нестационарности исследуемого про­ цесса.

Входные свойства одноемкостной схемы замещения опреде­

ляются выражением

 

 

 

 

Lo a (Bi, Fo) = exp

J 2

/

J J - J B i \

(1.99)

 

 

 

a

\

Bi /

 

39

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