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

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

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

Из (1.27) следует,

что

 

 

^ ( p ) - V i ^ f m ^ -

( L 2 9 )

При измерениях

с Хи = О Zc (р) = ] / и

у (р) = ] А | - -

Если тепловая связь между боковой поверхностью первичного преобразователя и средой слабая, т. е. а х = 0, то у (р) = ] / ^ ~ -

и Zc (р) = (рЛср)

Функциональные зависимости 7 (р) и Zc (р) характеризуют динамичность процесса измерения. Произведение у (р) Zc (р) яв­ ляется постоянной величиной, численно равной Исходя из

этого можно определить одну из названных функциональных за­ висимостей, если известна другая.

При анализе и расчете температурных полей широко приме­

няются критерии Fo =•

и Bi =

Критерий Фурье является

внутренним критерием,

величина которого определяется исходя

из значений физико-геометрических

параметров

твердого

тела

(в нашем случае — первичного преобразователя).

Критерий

Био,

или критерий краевого подобия, характеризует интенсивность теплообмена между контролируемой средой и первичным преобра­ зователем. При Bi <С 1 имеет место малая интенсивность тепло­ обмена и большая теплопроводность, что позволяет определить состояние контролируемой среды по изменению температуры пер­ вичного преобразователя во времени. Если Bi > 1, то температура первичного преобразователя равна температуре среды.

Величина /, входящая в критерии Fo и Bi, является характер­ ным размером и определяется как отношение объема тела к пло­ щади его поверхности. Для простейших геометрических тел зна­

чение / определяется так: для шара и для цилиндра

/* С

для

куба

(L — длина ребра куба; г — радиус шара

и основа­

ние

цилиндра).

 

Стационарный и квазистационарный процессы измерения пол­ ностью определяются значениями параметров Fo и Bi . Динамич­ ность измерительного процесса оценивается величинами комби­

наций параметров Fo и BiFo.

 

 

Если Fo <С 1 и

BiFo < 1,

измерительный процесс1

подчи­

няется особенностям

граничных

условий третьего рода. При этом

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

имеет место переходный

режим.

1 Измерительный процесс — изменение величины и формы регистрируемого (вторичного) сигнала в зависимости от соотношения величин критериальных параметров, характеризующих исследуемый процесс,

10

Если Fo > 1 и BiFo < 1, измерительный процесс также подчи­ нен особенностям граничных условий третьего рода, но при этом

термоприемник можно

рассматривать

равномерно

прогретым.

Если BiFo > 1, в измерительной среде выполняются

особенности

граничных условий первого рода. Если

при этом Fo > 1, то

пер­

вичный преобразователь

равномерно прогрет; при

Fo < 1

как

в среде, так и в первичном преобразователе имеет место переход­ ный процесс. Если BiFo < 1, температура среды является гене­ ратором теплового потока для первичного преобразователя. При этом величина q в основном зависит от значений tc и а и практи­ чески не зависит от свойств первичного преобразователя. При BiFo > 1 температура первичного преобразователя приближенно равна температуре среды. В этом случае среду можно рассма­ тривать как генератор температуры для первичного преобразо­ вателя.

Наиболее неблагоприятные для измерений условия имеют место при B i Fo? « 1 и Fo ^ 1. При данных соотношениях кри­ териальных комбинаций параметров Bi и Fo в измерительной среде имеет место сугубо нестационарный процесс. Это объясняется тем, что при BiFo ^ 1 и Fo — 1 в организации теплового потока в равной мере участвуют свойства как среды, так и первичного преобразователя. Учет влияния первичного преобразователя на контролируемый процесс рассмотрен отдельно в п. 4.

2 . О С Н О В Н Ы Е Ф У Н К Ц И О Н А Л Ь Н Ы Е З А В И С И М О С Т И

ИП А Р А М Е Т Р Ы , Х А Р А К Т Е Р И З У Ю Щ И Е

Д И Н А М И Ч Е С К И Е С В О Й С Т В А И З М Е Р И Т Е Л Ь Н О Г О П Р О Ц Е С С А

Взаимосвязь регистрируемого (вторичного) сигнала и контро­ лируемого процесса устанавливается двумя основными методами; методом спектрального анализа (методом интеграла Фурье) и методом решения дифференциальных уравнений, описывающих нестационарную теплопроводность в измерительной среде. Вы­ бор метода определяется степенью обозримости контролируемых особенностей исследуемого процесса.

Второй метод, приспособленный для анализа переходного процесса, использует в большинстве случаев операторную форму1 .

Анализ переходных процессов этим методом сводится к ана­ лизу решений уравнений теплопроводности, записанных в виде 2

 

 

 

 

КА*) = -¥тт

=К,(Р)

 

 

(1-30)

 

 

 

 

 

 

(Р)

 

 

 

 

 

1

Обычно выражения, характеризующие динамические свойства исследуе­

мой среды,

записываются

при условии

tc =

0, т<=;0;

tc =

1,

т > 0 .

2

Если

Bi

оо или

BiFo >

1,

то

Кв

(Bi,

Fo) =

KJFo)

и"То3 (Bi, Fo) =

= Kq

(Fo),

где

Ka (Bi, Fo) =i Кв

(p)

и

Lo 3

(Bi,

Fo) =

Lo3

(p).

11

L O 3 ( T ) ^ < ^ - , x = 0, т ^ О .

(1.31)

tc (P)

 

Анализ процесса теплопередачи в измерительной среде методом интеграла Фурье использует возможность представления воздей­ ствия и свойств исследуемой среды в виде суммы простых гармо­ нических колебаний различной амплитуды и фазы, частоты ко­

торых

изменяются непрерывно от со =

0 до со = оо.

 

Частотные и фазовые свойства измерительной среды описы­

ваются

следующими функциональными

зависимостями:

 

и

С (/со) =

| С (/со) | ехр [/ф (со)]

(1.30а)

 

 

 

 

 

Lo3(/co) =

х - 0 ,

со^О,

(1.31а)

где | С(/ со) | и | Lo 3 (/со) | — амплитудно-фазово-частотные

характе­

ристики измерительного тракта соответственно по температуре и тепловому потоку. Характеристика Lo 3 (/со) представляет собой зависимость амплитуд и фаз гармонических составляющих сиг­ нала q (т, 0) при tc = 1 от частоты. По своему физическому со­ держанию Lo 3 (т) является входной проводимостью измеритель­ ного тракта. Характеристика С (/со) является комплексным коэф­

фициентом

передачи

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

тракта.

 

Связь между временными характеристиками контролируемого

процесса (измерительного тракта) и спектральными

свойствами

вторичного

сигнала

устанавливается

зависимостью1

 

 

 

т,

 

 

 

Я (/со) =

J е (х) ехр (— /сох) dx,

(1-32)

 

 

 

to

 

 

где е (х) — исследуемая

функция. Моменты времени

т 0 и хх со­

ответствуют началу и концу действия е (х). Переход оте (т) к Н (/со)

осуществляется

по

формуле

 

 

е(т) =

4 " 1 1#(/ю)|ехр[/(<от - ф)]Ж»

(1.33)

 

 

—со

 

 

или

 

 

 

 

 

 

+ о =

Т,

 

е (т ) ~

~2лГ

1 е х Р ( / ( О Т ) ^ ( 0

j е(х)ехр(—j(x>x)dx.

(1.34)

та

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

12

Если

 

 

 

Я (/со) = | Я (/со) | ехр [/ф (со)] =

Я х (со) + / Я 2 (со),

(1.35)

то

 

 

 

 

со

 

 

Я (/со) = Я , (0) +

А } - ^

L cos (сот) dco

(1.36)

 

о

 

 

или

 

 

 

 

со

 

 

Я(т) = 4

f - ^ - sln(a>T)dT .

(1.37)

 

о

 

 

Из (1.33) — (1-37) следует, что переходный процесс может быть полностью определен зависимостью от частоты одной веществен­ ной либо мнимой части стационарного параметра. Обратная за­ висимость в обобщенном виде подчиняется выражениям

 

 

со

 

 

Я х

(со) =

со J sin (сох) Я (х) dx

(1.38)

 

 

о

 

 

и

 

со

 

 

 

 

Я 2

(со) =

со j

cos (сох) Я (х) dx

(1.39)

 

 

о

 

 

или

 

 

 

 

 

 

 

со

 

Нх (со) = Я(о) +

| cos (сох) Я ' (х) dx

(1.40)

 

 

 

о

 

и

 

 

 

 

 

 

со

 

Я 2

(со) =

J sin (сох) Я ' (х) dx.

(1.41)

 

 

о

 

 

Формулы (1.38) — (1-41) позволяют определить спектральные свойства исследуемого процесса, если известны его временные ха­ рактеристики, и, наоборот, по временным характеристикам уста­ новить частотные свойства исследуемого процесса. Использование тех или иных характеристик дает возможность однозначно опре­ делить качество воспроизведения контролируемого процесса. Вы­ бор временного или спектрального метода исследования в каждом конкретном случае решается отдельно, в зависимости от нагляд­ ности наблюдения контролируемых особенностей изучаемого про­ цесса и сложности использования математического аппарата.

13

(*max).

Операторный метод представляет собой метод интегрирования дифференциальных уравнений. Вид их решения обычно представ­ ляется переходной функцией

 

 

Н(х)

=

Н{р);

 

(1.42)

 

 

СО

 

 

 

 

 

 

77 (р) = р J Я (т) ехр (— px)dx.

 

(1.43)

 

 

о

 

 

 

 

 

Переход (1.42) в общем случае осуществляется по формуле

Мелинна—Фурье

[90]

 

 

 

 

 

 

 

(7-f~/со

 

 

 

 

Я (

т ) = _ ] _ I

Я ( Р ) е х р ( р г ) ^

( М 4 )

 

 

а—/со

 

 

 

Связь между контролируемым и измеренным сигналами под­

чиняется любому из следующих

выражений:

 

 

 

 

 

 

т

 

 

 

Я и

(т) =

Я (0) А (т) +

J Я ' (т - T l ) А К )

(1.45)

 

 

 

 

о

 

 

 

 

 

 

 

т

 

 

 

Я и

(т) =

Я (0) А (т) +

J Я ' (та ) А (т -

т2 ) dTi ;

(1.46)

 

 

 

 

о

 

 

 

 

 

 

 

т

 

 

 

Я и

(т) = H(x)h (0) +

J Я (T l ) /г' (т -

тх ) dt i ;

(1.47)

 

 

 

 

о

 

 

 

 

 

 

 

т

 

 

 

Я и ( т ) =

Я(т)А(0) =

{ Я ( т — T ^ A ' f a ) ^ .

(1.48)

 

 

 

 

о

 

 

 

Если контролируемый процесс не может быть описан одним

аналитическим выражением в течение всего времени

его суще­

ствования, то

используется

кусочно-аналитическая

аппрокси­

мация с последующим применением метода элементарных вклю­ чений [90].

Элементарные включения (Э. В.) и начало их действия выби­ раются таким образом, чтобы их суммарное действие было экви­ валентно действию данного контролируемого процесса.

Рассмотрим основные определения и свойства импульсных процессов (сигналов). Контролируемый процесс при прохожде­ нии через измерительный тракт (среду) искажается. Для оценки степени искажения исследуемого процесса используются следую­ щие параметры.

1. Максимальное значение контролируемой температуры (те­ плового ПОТОКа) L

14

2. Практическая ширина импульса т ш , определяемая как ширина импульса 1С (т), измеренная на условной высоте (услов­

ном

уровне), т . е . d =

A 4 m a x -

Обычно

величина А равна

0,1

или

0,5.

 

 

 

 

3. Крутизна переднего и заднего фронтов зависимостей tc

(т),

взятых на его условном

уровне,

S± и

5 2 .

 

4.Время нарастания процесса 4, (h), т. е. время, в течение

которого tc (т) достигает максимального значения (или 0,9 макси­ мального значения).

5.Время спадания изучаемого процесса тн , т. е. время, в те­ чение которого исследуемый процесс спадает от своего максималь­ ного значения (или 0,9 максимального значения) до значения услов­

ной высоты, на которой определяется ширина процесса tc (т) (обычно до 0,1 максимального значения).

Перечислим коэффициенты

искажения вторичного

сигнала.

1. Коэффициент потери амплитуды (в обобщенном

виде)1

[Л,,

Ни

 

"max

 

"max

2. Коэффициент расширения вторичного сигнала по отношению

кконтролируемому процессу

3.Коэффициент потери крутизны

S h 2

4.Коэффициент замедления нарастания

5.Коэффициент замедления спадания

Следует отметить, что для прямоугольного импульса коэф­ фициенты ц не имеют физического смысла.

Связь формы спектра процесса | Я (/со) | с его изменением во времени устанавливается соотношениями:

 

+00

 

 

И { д

= -Ш-\

| Я (/ю)|ехр(/а>т)Ао;

(1.49)

 

— СЮ

+ 00

 

 

 

 

| Я (/со) | =

Я (со) -

} Я (х) ехр ( - /сох) dx,

(1.50)

* Индекс «и» относится ко вторичному сигналу (процессу).

15

где Я (со) — спектр обобщенного процесса, амплитуды которого увеличены в ]/2я раз.

Если Я (т) — четная функция, то

со

Я

(со) =

] /

~

— \

Н (x)cos(cox)cix;

(1.51)

 

 

 

я

о

 

 

 

 

 

 

со

 

 

Я

(г) =

]

/

j Я (со) cos (сох) cix.

(1.52)

 

 

 

я

о

 

 

Если Я (г) — нечетная функция, то

со

Я (со) =

] /

А

J

Я (х) sin (сох) dx-

(1.53)

 

 

 

о

 

 

 

 

 

со

 

 

Я (т) =

Т / [

я (со) sin (сох) dx.

(1.54)

 

г

я

о

 

 

Таким образом, форма графика Я (т) и форма его спектраль­

ной плотности

являются

взаимно-однозначно-обращаемыми.

 

Если

контролируемые

4 (т)

или

q (т)

длительностью

т Ш 1

воздействуют на измерительную

схему с

шириной

полосы

т ш ,

то при любых значениях

А со или т ш форма вторичного сигнала не

изменится

при

ти

Асо =

т ш , Асо! =

уи .

 

(1.55)

 

 

 

Из (1.55)

следует, что

 

 

 

 

 

 

 

 

где фн и ) — функция

передачи

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

 

3 . О П Р Е Д Е Л Е Н И Е П О В Е Р Х Н О С Т Н О Й Т Е М П Е Р А Т У Р Ы

 

Т В Е Р Д Ы Х Т Е Л П Р И Д И Н А М И Ч Е С К И Х И С С Л Е Д О В А Н И Я Х

 

Измерения

поверхностных

температур

твердых

тел произво­

дятся главным образом термопарами и термометрами сопротив­ ления.

Жесткие требования к инерционности данных первичных преобразователей трудно совместить с требованиями механиче­ ской прочности, в результате чего при измерениях в первичных преобразователях имеют место некоторое распределение темпера­ тур и отток тепла по проводам и с боковой поверхности. По этим

^л.р%чинам измеренное и контролируемое значения температур оказываются различными.

Для перехода от измеренного значения температуры к контро­ лируемому на практике используется уравнение (1.6). При этом

Ш

принимаются во внимание условия (1.13) — (1-15) или (1.16) — (1.18). Выбор условий определяется конкретными требованиями,

предъявляемыми к конструкции

первичного

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

и условиями

опыта.

 

 

 

Решениями

уравнения

(1.6) с

условиями

(1.13) — (1.15) или

(1.16) — (1.18)

являются

выражения (1.21) — (1.24). Они гро­

моздки и поэтому сложны для использования на практике. При определенных условиях можно пренебречь их некоторыми со­ ставляющими.

Если, например, значение комплекса - ^ - существенно меньше

единицы, то большая часть тепла, получаемого от контролируе­ мой среды, аккумулируется первичным преобразователем. При этом оттоком тепла по токопроводам можно пренебречь с доста­ точной для инженерных расчетов точностью. При условии, что

Osj ~т4- =s£ 5-10~2 , в выражениях (1.19) — (1.22) можно

поло­

ли/

 

 

 

жить ZT =

оо. Значение

параметра -щ- может быть найдено как

расчетным

путем, если

известны величины параметров /,

X, ZH

так и экспериментально. В последнем случае это достигается сравнением показаний регистрирующего прибора при тепловом контакте первичного преобразователя со средой с показаниями,

полученными при погружении преобразователя с токопроводами

в контролируемую

среду.

При

условии,

что ^- ^

5-10"2 ,

 

 

Ь в

 

 

 

к

 

значение

параметра

уравнении (1.6)

можно

приравнять

нулю.

 

 

 

 

 

 

 

Если

Fo <С 1, то при расчетах

температурных

полей

можно

использовать соотношение (1.23) вместо

(1.21), что также приводит

к упрощению расчетных соотношений.

Физически указанная за­

мена решений эквивалентна исключению из рассмотрения влияния отраженных волн на формирование измеренного сигнала. Дей­ ствительно, если входной сигнал весьма короткий, то за время его действия групповая волна не успевает дойти до конца пре­ образователя х = /, а следовательно, и отразится от него. Если входное воздействие изменяется скачком, то существенный ин­ тервал времени, входящий в выражение комплекса Fo, численно равен времени установления теплового равновесия в первичном преобразователе. Если же входное воздействие (контролируемая

температура)

изменяется во времени по произвольному закону,

то скорость тепловой волны

vr приближенно может быть рассчи­

тана по формуле

УТ = — ,

где та — время заметного изменения

контролируемого

процесса [90].

Термопары или термометры сопротивления выполняются в виде

стерженька

или

пластины.

Крепление первичных преобразова­

телей осуществляется либо зачеканиванием в твердое тело, либо приклеиванием или напылением на теплоизоляционную основу.

2 П- Л- Короткое

" 1 ':

1?

Выбор типа крепления обусловлен условиями измерений и об­ ластью исследований. Так, например, при измерениях в двига­ телях внутреннего сгорания или котлах первичные преобразова­ тели температуры обычно зачеканиваются в стенки устройства; при внешних аэродинамических исследованиях на газодинамиче­ ских моделях они приклеиваются или наносятся напылением. Независимо от способа крепления контактного первичного пре­ образователя при измерениях всегда имеет место отток тепла

х

 

 

 

 

 

 

 

 

 

 

 

 

 

X

Рис. 1.

Модель

термоприемника

Рис.

2.

Модель

термоприемника

и

график распределения темпера­

и

график распределения

темпера­

 

 

туры

в

ней:

 

 

 

 

туры

в ней:

 

 

/

— эквивалентный

стержень;

2 —

/

— стержень;

2 — теплоизолирован­

теплоизолированная

прослойка

(тер­

ная

прослойка

 

(термический

зазор)

мический

зазор)

между стержнем и

между стержнем

и стенкой

(телом);

стенкой

(телом);

3 — стенка; 4,5 —

3 — стенка; 4,5

— изменение

темпера-

изменение температуры в стержне при

туры

в

стержне

при x Q

<

и

 

I*

 

I2

 

 

т .

<

и т„ >

а

соответственно

 

 

 

'*

 

 

 

 

"

а

о

 

 

 

 

т_ >

соответственно

 

 

 

 

 

 

 

 

а

а

 

 

 

 

по токопроводам и с боковой поверхности. Кроме того, в зависи­ мости от условий измерений в первичных преобразователях мо­ жет иметь место переходный процесс. Все перечисленные факторы составляют погрешность измерения, поэтому необходимо уметь их учитывать. Перейдем к исследованиям соотношений, характе­ ризующих переходный процесс в измерительной среде. Необхо­ димо при этом иметь в виду, что связь между t и /и устанавли­ вается интегральным уравнением Вольтерра, ядром которого является функция Грина или, как ее еще называют в технике, передаточная функция. Поэтому задачей исследований является изучение переходного режима в первичном преобразователе. Рассмотрим различные условия измерений,

18

 

Первичный

преобразователь,

у

которого

A =h 0

и

 

= О,

х

= I, схематизированно

показан на рис. 1. Здесь приведен

также

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I2

график

распределения температур

соответственно при

та

<С —

 

 

I2

 

 

 

 

 

 

ta

=h t из-за

 

 

 

 

 

и

та >

— .

В

первом

случае

переходного

процесса

в

первичном

 

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

Во

втором случае — вследствие

оттока тепла

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

в окружающую среду.

 

 

 

Переходный режим в первичном преобразователе с учетом

оттока

тепла

 

от

боковой

поверхности описывается

выражением

 

 

v

1-е

 

\

 

 

о

 

 

со

(га-f-0,5) л

 

 

 

 

 

 

1

 

/

 

\ V 4

 

 

 

 

 

К

< F ° ) =

T h X ~

2 е Х Р < -

ОТТ)

ZJ (П + 0,5)Я + Л» X

 

 

 

 

 

xsln[(n +

0,5)jt]exp{[—

( 2 п + 1 ) 2 ]

n a - J - } ,

 

 

(1.57)

где Л =

/

 

 

m

=

т ^ - ;

коэффициент

теплопроводности

окружающей

среды.

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/2

 

 

 

 

 

Из (1.57) следует, что при

та >-g^- t =

 

т. е. переходный

процесс

считается

завершенным.

 

 

 

 

 

 

 

 

При 0 ^

 

 

^

5 • 10"2

оттоком тепла от первичного преобразо­

вателя к телу, с которым он находится в тепловом контакте

вдоль

боковой поверхности, можно пренебречь. В этом случае из (1.57) имеем

/((Fo) = 1 — ^ sin [(п + 0,5) я] ехр [ — ( 2 л + 1 ) 2 " 2 - j L ] • (1-58)

За

 

 

 

Если та <^ -я—, то из (1.57) можно получить следующую формулу;

K(Fo).

ехр (A) erf (^у=

+ V™

 

e x p ( - ^ ) e r f ( — -

>/тт

(1.59)

которая удовлетворяет

уже условиям

(1.16) — (1.18).

Схемати­

зированно первичный преобразователь, переходный процесс в ко­ тором удовлетворяет выражению (1.59), показан на рис. 2.

Если А = 0 и 0 < /

< : 1, то выражение (1.59)

значительно

упростится:

 

 

к

^ = е г ! Ш -

( L 6 0 )

2*

 

19

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