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

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

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

направлении оси X (рис. 3.2). Единственная составляющая векторного-потенциала Ау удовлетворяет в проводящей области (для остальной области а = 0 ) уравнению

д2Аь

d2A,

дАъ

дАу

(3.22)

dx2

dz2

0.пУ - дх - = а ц о -

dt

с начальным условием Ay\t=0

Рис. 3.2. К расчету распростране­ ния импульса электромагнитного поля вдоль движущегося электро­ проводящего полупространства.

=0 (3.23) и граничными условиями

УІ ^ У 1 ! ^ ;

 

 

дАу*

 

dz

z-h

dz

 

= ц0І(і)[8{х-а)

-8(х+а)];

 

 

 

(3.24)

A 1 1 1

уш\

2=0

'

a y 1

 

& У "

 

dz

z=0

dz

2 = 0 '

 

 

 

где

Ау\ Ау11,

Ау

решения уравнения (3.22) для областей

z^h,

O^z^h

и z<0,

соответственно. Кроме того, дополнительно

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

тегральные преобразования Лапласа и Фурье к уравнению

(3.22)

и граничным условиям (3.24) с учетом начального условия

(3.23)

и условия ограниченности совершенно аналогично, как и для пре­ дыдущей задачи, получим уравнение

•у2Ау = 0

(3.25)

дг2

 

(для областей I , I I и I I I у2 = К2) с граничными

условиями

У и = У Ч z =

/t'

 

 

 

 

 

 

a y

 

a y 1

 

 

 

• sin

a\\

 

z = 7i

 

dz

 

 

np

(3.26)

dz

 

 

 

 

 

 

 

У І І г = 0 = У І І І z=0 '

 

 

 

 

 

 

dz

 

дАуш

 

 

 

 

 

 

z=0

dz

z=0

 

 

 

 

 

 

 

 

 

 

 

 

 

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

(3.25) для всех трех

областей

ищутся в сле­

дующем виде:

 

 

 

 

 

 

 

 

 

 

 

V = C , e - ^ - f t > ;

 

 

 

(3.27)

 

I v n = C 2

c h

Xz + C3

sh Kz;

 

 

(3.28)

 

 

I y ra=C 4 ev^ .

 

 

 

 

(3.29)

Д л я определения постоянных Сі получаем следующую

систему

алгебраических уравнений:

 

 

 

 

 

 

 

C, = C 2 c h U - t - C 3 s h U ;

 

 

 

 

 

1С у + ХС2 sh Kh + С3Х

chXh=-

np

sin Ka;

(3.30)

 

ХС3 =уС4

 

 

 

 

 

 

Решение имеет вид

 

 

 

sin

la

Що I

 

 

А 1-А И— —

 

 

 

 

np

Q-Hz+h)

A+y

2npK X

 

 

 

 

 

 

 

X [е-Мг-лі _ e-X(z+h)]

sin xa\

 

 

(3.31)

 

 

A v

m = - J * L e - " + * ^ .

(3.32)

Оригиналы выражений

(3.31) и

(3.32) находятся по формуле

об­

ращения для преобразования Фурье и таблицам для

операцион­

ного исчисления:

 

 

 

 

 

 

 

 

 

 

 

t

оо

 

 

 

 

 

 

Ау=-

f

dt

f е-^г +'1 Ш(А,,т)

s i n M * - U T ) sinXadX-i-

 

 

0

 

о

+ Ay°T-Ay™-t

 

 

(3.33).

 

 

 

 

 

 

 

 

t

 

°°

 

 

 

 

 

 

A y m = ^ ^ L

f

i x

 

[b-UiNfax)

sinX{x-vt)sinJMtt.,

(3.34>

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

 

/ 1 /

А,2т \

 

 

M(%, T)

=

• -а./лтщ; erfcl

V

I ;

 

 

 

 

 

 

Уяацот

 

v

'

ацо'

 

 

N(Kx)

 

= —=-exp(

- —

 

I

-

 

 

o\x0

 

\ '

'

o\xo

 

 

Векторный потенциал «отраженных» токов

 

 

 

 

 

y

 

(z + h)2+(x-a)2

 

 

V

'

Векторный потенциал сторонних токов (потенциал бесконечнодлинной прямоугольной катушки с током (3.11))

А у = [n(z-h)2+(x-a)2 ' ( 3 ' 3 6 )

Напряженность электрического поля над поверхностью полу­

пространства

 

 

ГІД

 

 

 

 

 

 

-Я„от + £ у от (

 

(3.37),

 

Еу=

-

^

У*

 

где электрическое поле, наведенное индукционными токами,

 

 

оо

 

 

 

 

 

Е

2fAp7

Г &

- U

z + h ) M n

t) smXix-vt)

smla&%\

(3.38>

 

л о

 

 

 

 

 

 

электрическое поле отраженных токов

£ T _ H £ m i n l £ ± ^ ! ± i f ± ^

( 3 39>

где

^ = 0; *>0;

электрическое поле сторонних токов

ц р / 6

( 0

{z-h)>+(x

+ a)*

 

 

 

П ( z - A ) 2 + ( * - a ) 2 '

( d

'

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

 

2

 

Г - |

_ М ! _ e r f (І/J?.) ] W

( 3 . 4 i )

Найдем значение интеграла, соответствующее первому члену правой части соотношения (3.41):

J ^ - J ^ f е х р Г - + Ш>1(1Л,= - — L , X

x«p[-*£]/«p[-(*V^+ 'V^>-

= _ _ _ L e x P [ - ^ l / e - d « .

V 0"Ho

ь ' —

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

-ствнтельной оси от b І/ до 0 и по мнимой — от 0 до °°:

At

 

 

 

 

О

г • со

 

 

-

 

2 У

^

=

_ L ^ (

6

y ^ ) ,

 

2t

\

*

At

'

 

*

v

'

At

тде функция W | й ~ | /

)

— интеграл

вероятности от комплекс-

лого аргумента — является табулированной функцией [7].

Второй член правой части соотношения (3.41) вычисляется элементарно с учетом, что Reft<0:

Третий член правой

части

(3.41) сводится

к двукратному ннте-

тралу

 

 

 

 

 

 

/ 3 =

[

\

[

e-"'da } e i M > U K

 

<W>V"o

1

о

}

 

д л я вычисления

которого предварительно

рассмотрим интеграл

 

 

 

со

 

ЧіІа»°

 

 

 

21

 

 

 

 

 

 

о"Ц.оУя0

 

 

0

 

Интегрируя по частям, получим

/4=

<2—=-'У-^—

 

\ ехр (iXb)

[

e - " ;

d « T

+

f L ^ y _

L

Гexp( - - i

- ^ +

u b

\ d L

Первый член для Re6<0 обращается в нуль, второй содержит уже известный интеграл / ь следовательно,

t

Для получения значения интеграла / 3 остается продифферен­ цировать / 4 по параметру Ь:

J 3 = -

d / 4

L-w(bV^)--Lw(b~\/^)

 

 

 

 

 

 

т

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

db

a\.i0b2

 

v

'

4t

і

2t

 

*

'

• 4t

'

bino\x0t

 

 

 

nf]2

 

\

1

At

I

9t

 

\

У

• At

/

 

 

Таким образом,

окончательное

значение интеграла

(3.41) по­

лучается из суммы уже вычисленных интегралов:

 

 

 

 

 

 

 

 

J = J\—/2

+ ^3

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

*

 

 

1

 

W

 

At

 

 

 

 

a\x0b2

 

bfnoiiot

 

 

оц0Ь2

 

\

1

 

 

Вводя безразмерные переменные

 

 

 

 

 

 

 

 

 

-

і

 

-

,

 

-

4t

 

 

_

 

 

 

 

 

х=х/а;

 

z = z/a;

 

t=

a\i0a2

-;

h—h\a\

 

 

 

 

 

1

 

 

 

-

 

„ „

 

 

 

 

 

 

v=—

a[x0av;

 

2лсга2

;

 

 

 

 

 

 

 

EyH=—-—£„H

 

 

 

 

 

si —

1 — (ic

— Є?ї);

 

b=l+(x-vt);

 

 

l = z + h

 

и учитывая (3.38) и (3.42), запишем выражение для безразмер­ ной функции, определяющей электрическое поле наведенных то­ ков при z2s0:

v

i t

і

V

j t

'

y

i t

'

v

- - Щ - *

i t

1

( 3 .43>

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

v _ , [ ^ ( - * ) - - £ . p ( - j i ) M £ - £ l

(3.44>

JffiZ-

V"°.25

Щ

~ V=0.50

 

V= 1.50

V=(,50

 

.V=i,00

 

V=0,50

I I I / \ /

\ s

 

\ V=0,2.5

4 11 л V 1 \ /

v=0,K V=0,25 V=0

4

V=l,50 V=(,00

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

Остальные

члены

(3.39)

и

(3.40),

содержа­

щиеся

в

выражении для

электрического

поля

(3.37),

существенной ро­

ли не играют. «Отражен­

ные»

токи,

создающие

отраженное

 

поле

Еуот,

как

видно

 

из формулы

(3.39), зеркальны

по от­

ношению

к

первичным,

если

зеркалом

служит

поверхность

полупрост­

ранства.

«Отраженное»

поле

появляется

только

в начальный

момент и не

зависит от скорости. Им­ пульсная функция, содер­ жащаяся в выражении (3.37), вносит некоторую неопределенность при t= = 0, так как при решении

Рис. 3.3.

Зависимость

функции Evu

от t и v при

а 1=0; б — £ = 0 , 5 ; в —

С=2 .

задачи

предполагается идеальный скачок векторного потен­

циала

при t = 0. Практически эта неопределенность исчезает при

учете индуктивности контуров. Экспериментально легко добиться,, чтобы время существования наведенных токов, которое опреде­ ляется критерием

T=o[i0a2,

было гораздо больше продолжительности прямого взаимодейст­ вия между катушками. На рис. 3.3 приведены кривые зависимо­

сти Ё,/ от

t с

различными значениями v для х — 2 и

соот­

ветственно

£ = 0;

0,5 и 2,0. Расчетные кривые показывают,

что-

поле существенно затухает уже при незначительном удалении от

поверхности полупространства. Формула

(3.43) имеет одно важ ­

ное свойство

 

E v ' \ x . v t = °-

(3-45)'

Кроме того, нулевая точка, определяемая формулой (3.45), как' это непосредственно видно из зависимостей, приведенных на рис. 3.3, находится между двумя существенно выраженными мак­ симумами. Последнее обстоятельство позволяет легко определить, положение нулевой точки экспериментально. Измерив время

г = — и зная величину х, можно определить скорость движения: среды V.

§ 2. И М П У Л Ь С Н Ы Е И З М Е Р И Т Е Л И С К О Р О С Т И Э Л Е К Т Р О П Р О В О Д Я Щ И Х С Р Е Д

И И Х О С Н О В Н Ы Е Х А Р А К Т Е Р И С Т И К И ;

Расчетные кривые зависимости наведенной э.д. с. (см. рис. 3.1), которые подтвердились экспериментальными измерени­ ями, показывают, что одновременно с изменением скорости в из­ мерительной катушке изменяется амплитуда сигнала и происхо­ дит перемещение основных максимумов. Естественно, что стольотчетливые изменения могут быть использованы для определения скорости движения среды. Учет амплитудных изменений не при­ водит к сколько-нибудь эффективным результатам, так как ха­ рактер этих изменений сложен, и практически трудно добиться устойчивых зависимостей. Уменьшение ширины измеряемого им­ пульса при увеличении 'скоростей является более отчетливой характеристикой. Так, например, как показывают расчетные кри­ вые зависимости э. д. с. от времени, в случае передачи импульса между двумя круговыми контурами (см. рис. 3.1) при увеличеншг скорости среды происходит уменьшение времени появления мак -

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

Изложенное выше позволяет предложить импульсный способ измерения расхода по индикации максимума э.д. с. в измери­ тельной катушке [8]. Принцип работы устройства следующий (рис. 3.4). На трубопроводе с контролируемой жидкостью через некоторый интервал друг от друга расположены катушка воз­ буждения (передающая) и измерительная (приемная) катушка. Передающая катушка питается от генератора прямоугольными •импульсами тока. Импульсы, наведенные в измерительной ка­ тушке, не отличаются от рассчитанных по приближенной модели (см. рис. 3.1), когда все окружающее пространство заполнено движущейся проводящей средой. Сигнал с приемной катушки предварительно усиливается и подается в устройство, формирую­ щее прямоугольные импульсы, передний фронт которых совпа-

а; Передающая

^_ Р _^

Приемная

6)

катушка

котишка

 

Генератор

Л

Уси

импульсов

 

 

Л

Формирующее

Формирующее

устройство

устройство

1Измерительвременных интерЙолоВ

Рис. 3.4. Блок-схема измерителя расхода электропроводящей жид­ кости по способу индикации сдвига максимума наведенной -э. д. с.

визмерительной катушке.

аблок-схема; б — эпюры тока и э. д . с. в передающей н приемной катушках.

Рис. 3.5. Блок-схема импульсного метода измере­ ния скорости электропроводящей жидкости.

дает с максимумом импульсов, наведенных в измерительной ка ­ тушке. Устройство для измерения времени определяет интервалы т между моментами становления тока в передающей катушке и появления переднего фронта формируемого импульса, интервал времени обратно пропорционален скорости среды. На рис. 3.4,6" изображены эпюры тока и э. д. с. в передающей и приемной ка­ тушках. Условие максимума функции, определяемой формулой (3.21), довольно сложно. Приравнивая к нулю производную' функции, нетрудно убедиться, что при малых значениях і„ и боль­ ших у* условие максимума определяется соотношением

L=—;

t=—.

(3.46).

V*

V

 

Из расчетных кривых, приведенных на рис. 3.1, непосредственновидно, что соотношение (3.46) начинает выполняться с достаточ­ ной степенью точности с У» = 8, т. е. при сравнительно большой величине магнитного числа Рейнольдса.

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

Значительным преимуществом обладают импульсные методы

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