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

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

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

Приведем условия, при которых имеют место те или иные осо­ бенности измерительного процесса. Так, при 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

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