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

книги из ГПНТБ / Бесконтактный контроль потока жидких металлов

..pdf
Скачиваний:
10
Добавлен:
23.10.2023
Размер:
19.29 Mб
Скачать

часты уравнения (5.11). Тогда для нормальной компоненты вто­ ричного поля имеем

— Ь

Общее решение

уравнения (5.13)

с учетом

граничного

условия

(5.12) имеет следующий вид:

 

 

 

 

 

 

 

 

а

 

 

Ь

 

 

iiioob

/ s h a ( x - f o )

Г

Г

а

Г

ч ,

]

 

 

— а

 

 

—о

 

 

 

х

Ь

 

 

 

 

 

X s h a ( a - | ) d | - / [ t o - — /

и (g, (/) ch/ ]sh a ( x - g )

} .

Зная нормальную компоненту индукции на поверхности при­ емного индуктора у = Д, нетрудно определить величину э.д. с . наведенную в одном витке приемной катушки:

а

с

 

<g = - — J

J bv (х, z) dxdz.

 

Значение интеграла \'bu{x)dx

может быть получено

интегрпро-

—а

 

 

ванием (5.13) по х в пределах от —а до а.

 

Таким образом, после несложных преобразований

получим

a

'^ff

 

 

 

 

[ by (х) dx =

Вc m

[ со ( S h а а

-

ch аа )

•I

а2А sh 2aa

 

L \ aa

 

'

-a

 

 

 

 

 

a

b

 

 

 

 

4ab

f v(x'y) c^

oi,xdxdy — 4:abv0

ch aa j

~~h if

—a —о

 

 

 

 

-

j .

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

аЬ

c h a a - [ j ( f v(x, у) Ay} ch с ш і х ] [ 4a&o0 ]

Для турбулентного течения может

 

 

 

 

 

 

быть принят степенной закон распре­

 

 

 

 

 

 

деления скорости по сечению канала.

 

 

 

 

 

 

Так

как

в

выражение

(5.14)

 

входит

 

 

 

 

 

 

скорость, усредненная по координате у,

 

 

 

 

 

 

достаточно

задаться

только

зависимо­

 

 

 

 

 

 

стью скорости от координаты

х

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5.15)

 

 

 

 

 

 

По

формуле

(5.14)

была

вычислена

О

 

2

4

оса

 

 

 

 

 

 

. гг

 

со — avo

Рис.

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

отно­

относительная погрешность А V =

со

сительной

погрешности

от

при

 

 

 

 

 

 

 

 

аа

для канала конечной ши­

определении средней скорости ПО

 

 

 

 

рины.

формуле

(5.1)

для трех случаев: лами­

/

ламинарное течение;

2

турбулентное

течение,

т = 0 , 2 ;

нарного

течения (и = им(\—х2))

 

 

и тур­

3

турбулентное течение,

ш =

булентного

течения

при пг — ОА и т =

 

 

 

 

=0,1.

 

 

 

 

 

 

= 0,2. Как

видно

из

рис.

5.2, где

 

 

 

 

 

 

представлены результаты расчета, погрешность тем меньше, чем равномернее распределение скорости по сечению канала, причем погрешность измерений может быть снижена путем увеличения а/т, т. е. ширины канала. Как следует из (5.14), эта погрешность связана с симметричными отклонениями профиля скорости от равномерного распределения. Даже большие асимметричные ис­ кажения профиля скорости не будут вызывать в рассматривае­ мом случае погрешности измерений.

Таким образом, увеличение ширины канала либо применение хорошо проводящих боковых стенок канала позволяет значи­ тельно уменьшить погрешность расходомерного устройства, оп­ ределяющего расход по синхронной частоте: co=«t>o. Из формулы (5.14) следует, что погрешность в случае канала конечной ши-

рины обусловлена отличной

от нуля производной - ~

. Из­

вестно

[4], что -^-^1

обращается в нуль в том случае,

если

 

ОХ х=±<х

 

 

 

канал

ограничен идеальным

проводником. Все же применение

токозамыкающих шин может привести к появлению дополнитель­ ной погрешности. Так как проводимость шин не равна бесконеч­ ности и шины расположены в рабочем зазоре датчика, при а = а у 0 вторичное поле, даже при равномерном профиле скорости, в зоне движущейся проводящей среды не равно нулю [4]. Это связано с наличием индуцированных в шинах токов. Чем больше отноше­

ние проводимости шин к проводимости контролируемой

среды,

тем меньше будут проникать индуцированные в шинах

токи в

її*

 

область жидкого металла. Как уже отмечалось выше, в предель-

/СТщ

\

ном случае

vco наводимые в шинах токи не вызывают

\СТЖ

/

дополнительной погрешности. Следовательно, чтобы ослабить ин­

дуцированные в шинах

токи, шины целесообразно

размещать

вне зоны датчика.

 

 

Случай произвольной

величины относительного

зазора дат­

чика. Как следует из сказанного выше, при малой величине от­ носительного зазора датчика 2А/т влияние неравномерности про­ филя скорости по толщине канала незначительно. Задача об оп­ ределении этой погрешности при произвольном зазоре датчика решена аналитически в предположении малых магнитных чисел Рейнольдса и бесконечной ширины канала.

Полагая, что токи в проводящей жидкости распределены в виде плоских слоев, параллельных стенкам канала,'легко полу­

чить интегральное уравнение для векторного

потенциала

 

 

 

 

 

ь

 

 

 

А(я)=Аош{у)+

 

То

A

I

[ w - a « ( 6 ) ] c h o ( E + A)

X

 

a sn /аА ^ v

 

 

 

X c h a O / - А)Л (l)de-f- / [ ( o - a r j ( £ ) ] c h a ( | - A ) X

 

 

 

 

у

 

 

 

 

X cha(y

+ A)A(l)dl]

,

(5.16)

где

 

 

 

 

 

 

 

 

Аст{у)

=

-

Hojo

c h a ( t / - A )

(5.17)

 

a

sh 2aA

 

 

 

 

 

 

 

Как и в предыдущем случае, предполагается, что влияние вторич­ ного поля на распределение токов в движущейся проводящей среде пренебрежимо мало. Тогда в правую часть интегрального уравнения (5.16) необходимо вместо A (g) подставить Л с т ( £ ) . Следовательно, для вторичного потенциала на поверхности при­ емного индуктора имеем

ь

Л в т |

= Jsh^aA

/ ^ - a y ( s ) ] c h a a + A ) c h a ( s - A ) d | .

 

 

(5.18)

Вторичная э.д. с. обратится в нуль на частоте

 

 

 

 

и

 

 

 

 

 

 

 

с п 2 а Д + [ J

v(y) ch2at/dt/j

[2bv0

 

 

 

Ю

-

а Щ

ch2aA + s h 2 « 6 ( 2 a 6 ) - '

'

( 5

Л 9 )

По формуле

(5.19) была определена относительная погрешность

AV = -——

 

для

ламинарного

и турбулентного

течения

 

при

т=

= 0,1 и m = 0,2. Как показали расчетные данные (рис. 5.3), зави­ симость погрешности от а имеет максимум, положение которого зависит от коэффициента заполнения рабочего зазора датчика движущейся проводящей жидкостью k и формы профиля скоро­ сти. Погрешность стремится к нулю при устремлении а к нулю (плоскопараллельное поле) либо к бесконечности. Резко умень­ шается также погрешность с уменьшением коэффициента запол­ нения к. Последнее связано с выравниванием нормальной компо­ ненты магнитного поля по толщине канала с уменьшением его относительной толщины.

3

т

f

 

• Л .V• •• •і' . '

 

\

 

 

1

 

 

 

ш

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

/— намагничивающий индук­

тор; 2 —

участок

трубопровода

с контролируемой

средой; 3, •/ —

приемные

катушки.

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

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

 

 

Л в т Іу=д = (С, — I) c h a ( A - ' - l ) .

(5.20)

Зная частоту,

при которой Л в т обращается

в нуль, нетрудно оп­

ределить искомую погрешность.

 

 

Результаты

расчета для турбулентного

и ламинарного тече­

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

на рис. 5.5. Из результатов

расчета следует, что

погрешность, связанная с неравномерностью профиля

скорости

по толщине

канала, противоположна по знаку погрешности, об­

условленной

неравномерностью профиля

скорости по

ширине

канала. При Rem ^0,75 для турбулентного течения и Rem^0,3 для ламинарного эта погрешность мало зависит от магнитного числа Реннольдса.

Погрешность при ламинарном течении не превышает 5% при коэффициенте заполнения, равном 0,75. Естественно, что с умень­ шением коэффициента заполнения погрешность уменьшается, при этом зависимость погрешности от a становится менее выражен­ ной. Для турбулентного течения при т = 0,2 погрешность не пре­ вышает 1,6%, а при пг = 0,1 — 1%.

Очевидно, выбором соответствующей геометрии канала и дат­ чика в режиме турбулентного течения можно получить погреш­ ность, не превышающую 1%. На рис. 5.6 приведены расчетные зависимости относительной погрешности от а для бесконечно ши­

рокого канала для турбулентного

течения при R e m = l и для ла­

минарного течения при Rem = 0,5.

Как видно из рисунков, с уве-

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

Действительно, при Rem = l на синхронной скорости

глубина

проникновения при

толщине канала

10 мм и полюсном

делении

датчика т = 3 см

составит около 7,8

мм. Дальнейшее

увеличение'

Rem (скорости)

приводит к увеличению частоты питания

датчика,

а тем самым и к уменьшению глубины проникновения

магнитного

поля. Следовательно, растет и погрешность измерения

средней

скорости (расхода)

электропроводящей

среды.

 

 

Случай канала

конечной ширины.

Выше была

рассмотрена

погрешность расходомеров вследствие неравномерности профиля скорости для бесконечно широкого канала и для канала конеч­ ной ширины при малом рабочем зазоре датчика 2Д/т. Погреш­ ности в указанных случаях оказались противоположными по знаку.

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

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

у=

у 0 е *(ш<-ок).

Магнитная проницаемость

сердечника принимается равной бес­

конечности. Значения магнитных чисел Рейнольдса Rem

>Cl.

Ищется

частота,

при кото­

рой э.д. с, наводимая вто­

ричным

полем в

приемных

катушках,

расположенных

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

приемного

индуктора,

' обращается

в

нуль.

 

 

 

 

Рис. 5.7. К постановке задачи об опреде­ лении погрешности индукционного расхо­ домера.

Случай произвольной ве­ личины относительного воз­ душного зазора рассмотрен в работах [5, 6], где принято, что магнитное поле возбуж-

дения неизменно по ширине канала и обращается в нуль на гра­

нице

канала

при х=±а.

Предполагается,

что вторичное поле

также обращается в нуль на границе канала

(при х= ± а ) .

В

работе

[6] для равномерного профиля

скорости показано,

что вторичное поле вне канала и //^-компонента поля всюду равны нулю. Подобный результат имеет место и для неравномер­ ного профиля скорости, если Rem<§Cl [7].

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

бесконечного ряда по магнитному числу Рейнольдса

 

« - 2

Rem *Hv .

(5.21)

 

 

v=0

Уравнения для компонент индукции магнитного поля имеют

следующий вид:

 

 

 

 

 

 

 

 

дчх

д2Вх

-

 

 

 

 

 

 

- ^ + - ^ - x 2 5 . - c = - / R e m A l / 5 i ; ;

(5.22)

^ +

^ ~ ^

у

=

~ 1 R

e m

A

^ ;

( 5 - 2 3 )

 

 

дВх

 

дВу

 

 

 

 

где

- № + — ^ - + — ^ = 0 ,

 

 

(5.24)

л

= | 1 ; ^

а [

ш

^ { а

~

а

Щ ) ] ;

(5.25)

 

 

 

VQ

 

 

 

 

 

Для областей I и I I I необходимо положить а = А У = 0 . Воспользу­ емся разложением (5.21) и дифференциальными уравнениями (5.22) — (5.24), чтобы получить систему уравнений для Bv .

Таким образом, приравнивая коэффициенты при одинаковых степенях Rem , получим 'Следующую систему уравнений:

для областей I и I I I

+ — з Н г - - » 2 5 * ^ 0 ;

( 5 - 2 7 )

ag2 +

drf

- a 2 5 y v i = 0 ;

(5.28)

 

 

 

 

•iaB„*+

dB™1

.

dBvvJ

=0;

(5.29)

£

+ '

^

для области I I

 

 

 

 

 

х 2 в ї » п = - і А г а І ( » и і ) 1 1 ;

 

(5.зо)

 

- « 2 5 w n =

-iAVBy^

;

(5.31)

- i a S ; v n +

= £ — +

^ — = 0 .

 

(5.32)

Решения для областей I и I I должны удовлетворять следующим граничным условиям:

5 1 = 5 П ,

,

В1 и

,

В1 П.

;

(5.33)

Bxv* = BxvlI

= Bsvl

= B t v u = 0

при

г) = ± 1 ,

vSsl;

(5.34)

5 x o I = B * o u = B z 0

I = B z o I I = 0

при -Ї-І==1;

 

(5.35)

Ba^B^jo

 

 

при

н = - 1 ;

(5.36)

5^ 1 =5 ^ 1 =5 ^ , 1 = 0

при

l-voo,

v ^ l ;

(5.37)

 

[rotBv n] .

=0 .

 

 

(5.38)

При такой постановке задачи решение в нулевом

приближении

не отличается от полученного в работе [6]. Таким образом, если

53-0=0, то из уравнений

(5.30) и

(5.32)

и соответствующих гра­

ничных условий следует, что 5 K V = 0 . Принимая во внимание, что

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

при­

ближении при RCn-Cl отсутствует 5к -компонента поля, из

(5.32)

и (5.38) следует

 

 

 

 

 

 

 

д ' В ^ 1

- a 2 5 y V n = 0

при

l=±k0,

vSsl;

(5.39)

с5ті2

 

 

 

 

 

 

 

дВу^

 

dByJ

= 0

при

л = ± 1 , vSsl.

(5.40)

С7ї)

дї]

 

 

 

 

 

С учетом условия

(5.40)

решение

(5.31)

можно искать в следую­

щем виде:

 

 

 

 

 

 

 

 

5 y v n

= 2jCn*(l)

cosX„T),

 

(5.41)

 

>yV

71=0

 

 

 

 

 

 

 

 

 

 

 

 

где Хп = пп.

Подставляя Byvn в (5.39), получим

 

B y v n = 0

при

l=±ko, v > l .

(5.42)

Решая соответствующие

дифференциальные

уравнения для

BlfV- и Szv-компонент поля при соответствующих

граничных усло­

виях, получим

 

 

 

B y v i = S I V = 0 ,

vSsl.

(5.43)

Таким образом, при неравномерном профиле скорости в ка- . нале конечной ширины, заполняющем весь рабочий зазор дат- ' чнка, поле вне канала оказывается невозмущенным. Кроме того, поперечная /^-компонента поля равна нулю и в области, занима­ емой каналом.

Следовательно, для нахождения нормальной компоненты ин­

дукции Ву

магнитного поля необходимо решить

дифференциаль­

ное уравнение

(5.31) при следующих граничных

условиях:

 

 

дВ „v

 

 

 

 

 

 

 

 

v

- 0

при

т і = ± 1 ,

v5=sl;

 

 

 

5„ v = 0

при

g = ± A 0 ,

vS==l.

 

Решение

уравнения

(5.31)

при

указанных

граничных

условиях

имеет вид

 

 

 

 

 

 

 

 

 

 

Bvv=2

Cnv{l)

cos А,„Ї],

 

 

 

 

 

 

•n=0

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

У An-+

У.'

о

 

 

 

 

 

 

 

 

ко

 

 

 

 

 

 

 

ch

Укпа+х%

 

 

 

 

 

 

 

 

 

і

 

 

 

 

 

 

 

/nv (6)

= ~ t

/ А ^ ( 1 . Л ) 5 „ С ^ 1 ) ( | , Л )

cosX„Tidti;

(5.44)

 

 

 

 

 

 

 

 

 

б а о = -

_ ^ р _ _ ^ Л _ +

\

c h y^a+^ag cos lnr\\

(5.45)

 

 

x

ch 2x

n = 0

 

 

 

 

 

 

 

 

 

 

 

 

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