![](/user_photo/_userpic.png)
книги из ГПНТБ / Бесконтактный контроль потока жидких металлов
..pdfнаправлении оси 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, |
Ау1П — |
решения уравнения (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 |
|
\ ' |
4т |
' |
o\xo |
|
|
||
Векторный потенциал «отраженных» токов |
|
|
|
|
||||||
|
y |
|
4я |
(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)* |
|
|
4я |
|
П ( 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, т. е. при сравнительно большой величине магнитного числа Рейнольдса.
Полученные теоретические результаты справедливы только для очень больших скоростей, поэтому в большинстве случаев необходима градуировка расходомера.
Значительным преимуществом обладают импульсные методы