книги из ГПНТБ / Бесконтактный контроль потока жидких металлов
..pdfчасты уравнения (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 |
|
|
|
|
|
|
|
|
|
|
|
|
|