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

книги из ГПНТБ / Мясников, Л. Л. Новые методы измерений в подводной акустике и радиотехнике

.pdf
Скачиваний:
14
Добавлен:
22.10.2023
Размер:
12.12 Mб
Скачать

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

на

поверхности хорошо проводящих тел

были предложены

М.

А. Леонтовичем и носят название

условий

Леонтовича.

 

В случае идеального проводника Ех

и Еу становятся исчезающе

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

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

{рис. 3.1).

Напряженность электрического по­ ля падающей волны можно предста­ вить в виде

Е г = Е йе~1

где г о— орт направления на точку наблюдения; s — орт волнового вектора; /*0. s = ах + ру + yz\ а, р, у —• напра­ вляющие косинусы. Соответствующие выражения могут быть написаны для напряженностей в отраженном и пре­ ломленном лучах.

Условия на

границе (г — 0) будут

иметь вид

 

 

 

_

 

Е■'-г1

е 2

~ ЕХ2,

*z2

Ч

 

 

Рис. 3.1. Падение плоскополяризованной монохроматической электромагнитной волны на гра­ ницу раздела двух сред.

Отсюда вытекают известные законы отражения и преломления:

sin# = sin i,

т. е. # =

i,

и

 

 

^2 _ -1 / Ч _

П2

* Г ~

V

 

где п[ и П2 — показатели преломления первой и второй среды соот­ ветственно.

В случае комплексного значения ег значение sin i также стано­ вится комплексным. Это означает, что преломленная волна пере­ стает быть плоской и плоскости равных амплитуд и равных фаз не будут совпадать.

Для низких, звуковых частот комплексная диэлектрическая про­

ницаемость

 

е' = е -f- i о

(3.1)

.

со

имеет сильно выраженную мнимую часть, причем k 0 =

— .

4

51

Для морской воды можно считать е = 80, а = 1 Сим/м и, таким

образом, для частоты порядка

1000 Гц

е' = 80

+ i - 18 • 106.

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

По определению, коэффициент отражения есть отношение ампли­

туды напряженности

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

поля отраженной волны Е 01

к амплитуде напряженности падающей волны

Е 0:

Коэффициент преломления равен

 

 

 

v —

S

 

где £ 02— амплитуда

напряженности

электрического поля прелом­

ленной волны.

 

 

 

Коэффициенты v

и £ могут быть

записаны

в виде

2kx cos #

k2cos # — kxcos r

М cos г -{- k 2cos О ’

kxcos r + k%cos i

Это — известные

формулы Френеля.

Связь между

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

напряженности магнитного

будет аналогичной, потому что

 

 

Нъ — У £2

Нг — у Гг2Ег.

(3.2)

поля

(3.3)

В случае поляризации в плоскости падения волны (т. е. когданапряженность электрического поля перпендикулярна плоскости падения, а напряженность магнитного поля лежит в ней), получаются аналогичные соотношения, только в этом случае удобно сначала рассматривать компоненты напряженности магнитного поля. Условия на границе при 2 = 0 будут иметь вид

Hxi = Нх2 при 2 = 0;

НА = На при 2 = 0.

Поскольку обе среды считаются немагнитными,

И-i = И-2 = 1 и Y

— 1 •

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

52

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

4

cos

у~е'2 ■

sin

4

cos

+ у g'

(3.4)

sin2 О

а коэффициент преломления для магнитных компонент

cos O' •

У 4 — sin2 0

cos 0

W

k2

vm=

 

sin2 ft

 

 

ъ2

cos О 4

 

 

 

: - Y e2 ■

 

cos 0

 

Ro

 

 

 

 

k2

 

 

 

 

 

sin 2 ft

(3.5)

sin2 0

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

1 оси 2:

 

 

дЕ

»■

дЕ,

 

 

 

дЕ2

чч

дЕ,

 

 

дг

дх

 

 

дг

*

ду

 

 

дН2

 

 

дН2

 

дН2

 

дН2

 

 

дг

 

 

дх

 

 

дг

 

^

ду

Из

уравнений Масквелла для обоих сред и из граничного условия

 

 

1НЬ1 = Ну2 при z — О,

считая в силу сказанного

дЕХ1

 

 

 

дх

малыми, получаем

 

дг

 

 

 

 

 

 

дЕг1

 

 

 

дН,

 

 

 

 

 

дх

 

 

 

 

 

У1

 

 

 

 

 

 

 

 

 

 

dt

 

(3.6)

 

 

дЕХ2 __

■ц0ц

< ? Я ,У2

 

 

дг

 

 

 

 

dt

 

 

При 2 = 0 правые части в

(3.6)

 

 

 

 

дЕ

 

численно^ равны, [значит д *‘

и

dEXi

также численно

 

равны.

Отсюда,

принимая в расчет выра­

дг

 

жение для плоской волны и граничное условие

имеем

Ех1 =

Ех2 при

z =

0,

Е,

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-*1

 

 

 

 

V

 

1

 

 

2яе0/

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5$

При

о

 

 

2 Я 8 08 2/ » 1

 

 

 

е г9 - -

 

‘• - К 1

 

 

 

 

V

7 2яе0/

(3.7)

 

Ex2

I EX1

 

 

 

1 EZI

 

тде / — частота.

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

tg 0 =

«он очень мал.

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

Перейдем теперь к рассмотрению электромагнитного поля источ­ ника, в качестве которого предлагается или электрический диполь (короткий по сравнению с длиной волны прямой проводник с пере­ менным током) или магнитный диполь (рамка с переменным током).

В теории распространения радиоволн, развитой Зоммерфельдом, вводится вектор Герца Z, удовлетворяющий волновому уравнению

 

d2Z

dZ г,

(3.8)

V2Z- ■£08fX0fX

— ‘ О-

Напряженности полей при этом выражаются

так:

Е =

rot rotZ;

|

 

Н =

е0е rot

. |

' 3-9)

В случае монохроматической волны, когда зависимость от времени «будет иметь вид eriat, уравнение для Z приобретает вид уравнения Гельмгольца:

 

AZ -j- &2Z =

0.

(ЗЛО)

Тогда

 

 

 

 

Е = grad div Z +

k2Z\

(З.П)

/

H — Uое0е rot Z.

 

54

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

k2 = k l ^ + i ^ y

(3.12>

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

Z = |0, О, Zz),

причем z — это вертикальная ось; Zz — функция координат г и г. В дальнейшем будем ее обозначать просто через Z.

Решение уравнения (ЗЛО) должно удовлетворять следующему

условию.

Если ввести радиус-вектор /?, имеющий длину R =

= | / r 2 +

z2, то в точке О решение должно обращаться в бесконеч­

ность определенным образом (как HR), а на бесконечности должно обращаться в нуль. Величины, относящиеся к воздушному полу­ пространству (т. е. при z > 0), будем отмечать индексом нуль. Таким образом, k 0 будет волновым числом для воздуха (оно веще­ ственно и равно ю/с). Тогда решение уравнения (3.10) будет иметь вид

2k2

pik0R

при

z >

0;

— ,2

Zo

r + K

 

 

 

(3.13)

2

eikR

 

 

2k\

 

 

 

0

Z'

при

z <

0.

Z ~ k 2 + kl

R

 

 

 

Функции Zo и Z' должны удовлетворять уравнению Гельмгольца для амплитуд. Будем искать выражения для этих функций в виде интеграла Зоммерфельда

Zo = J /о (Я) J0 (Яг) е ~ ^ ~ гdX]

°.

(314

Z' =f { XJ) J o ( X r ) e V v - v

о

где Jo (Яг) — функция Бесселя первого рода нулевого порядка. Данные выражения удовлетворяют дифференциальному уравне­

нию (3.10) и требуемым условиям,

если f 0(X) и f(X) убывают по край­

ней мере как

ИХ2 при

X —>оо.

 

 

Желательно выразить целиком Z0 и Z через интегралы. Это вполне

возможно,

так

как

из

теории

бесселевых

функций следует:

 

 

 

03

Л ---- уо (Яг)

1 Г 2

2

 

J

Г

v х ~feo g dX при г > 0;

R

О

V x 2- k l

 

 

(3.15)

 

 

 

 

 

 

 

 

 

со

 

 

 

 

 

 

 

 

 

 

 

 

 

4 r e ik* =

Г

,

Я

JolXrt eVv^ *dX при г < 0 .

R

 

J

Y x 2 — k2 ov

 

 

о

55

Если

2 = 0,

то

R = г.

 

 

 

 

В качестве предельных ставятся следующие условия:

1) Zokl =

Zkl

при 2 =

0;

 

 

 

2)

производная функции Z по z должна быть непрерывной.

В

результате

получается

 

 

 

 

 

 

?

X

(lr)e

- V

x 2 - k 2 Z

 

 

Z0= 2k2J

 

1

k° z dx-

(ЗЛ6)

Z = 2kl j- \ j , { X r ) e V v - k * z d K

о

где N = k2 ] / V — ko k20j / V — k2.

Выполняя операции над вектором Герца по формулам (3.11), можно найти компоненты напряженности магнитного поля Н и напряженности электрического поля Е.

Приближенные формулы для интегралов приводят к так на­ зываемой формуле Вейля — Ван-дер-Поля [37].

Вектор Герца для воздушной среды равен

 

 

 

 

 

 

 

 

 

Jk„R

 

 

 

 

 

 

 

 

 

Z 0 =

A _ _ s ,

(3.17)

где 5 — коэффициент

ослабления,

определяемый выражением

 

 

 

 

 

 

 

 

 

 

х + т

 

 

 

 

5 = 2 —4хе~<и+т>2 |

da,

(3.18)

в

котором

 

 

 

 

 

 

 

 

—/со

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

------т 711

 

 

 

 

 

 

 

_

R

е

4

 

 

] / '

к2 sin2 О ;

 

 

 

 

 

4

 

 

 

 

У Ч

 

 

sin

ft

(3.19)

 

 

 

 

 

 

 

 

 

3 .

 

 

 

 

 

 

 

 

 

 

Г Т Ш

 

 

 

 

 

 

 

 

 

 

^ —5- cosft.

 

 

 

 

 

 

 

 

 

 

sin й

 

 

 

Обычно

формулу

(3.18)

пишут

в

упрощенном виде,

полагая

 

 

 

 

 

 

 

 

 

 

^2

 

 

•sin 0 = 1 и

пренебрегая

 

величиной ~

под радикалом.

 

 

Обозначив

 

 

 

 

 

 

 

 

 

 

 

 

 

 

_

 

v 2

__

k 0R

/ k 0 \2

 

 

 

 

Р

 

 

 

 

 

 

 

и

имея в виду,

что т ^ О

(так как

c o sft^ O ), получим

 

 

 

 

 

 

2 — 4 У ре -р

Ур

 

 

 

 

5 =

j

ea2da.

(3.20)

— /со

Величина р носит название «числового расстояния».

56

На рис. 3.2 показана зависимость коэффициента ослабления 5 от «числового расстояния» р. Начиная с некоторых достаточна больших значений р коэффициент ослабления выражается в логариф­

мическом масштабе прямой линией.

s' 2

 

S/2J5

«Числовое

расстояние»

 

р

с

учетом

 

проводимости

воды,

когда

можно поло-

 

 

 

жить ' 2 < < i A

r

 

можно

 

 

приравнять

 

 

 

 

 

 

 

 

 

 

 

 

 

вещественному числу

 

 

 

 

 

 

 

 

 

я R

 

 

 

 

 

 

 

 

(3.21)

от -

 

 

Р =---------Т

 

 

 

 

 

 

V

*

+ (

2яе0/

 

 

 

 

0,0010,01 0,1

1,0 10 too

woo р

Если полагать

морскую

 

воду идеаль­

 

Рис. 3.2. Зависимость коэф­

ным проводником, то коэффициент 5 бу­

дет постоянным

и

S -----

2.

 

 

 

 

фициента ослабления от «чис­

Если

к +

т велико

по

сравнению с

лового расстояния»

для воз­

 

духа.

 

единицей,

то на том

участке,

где

подын­

 

 

 

тегральная функция заметно отличается от нуля, формула (3.18) переходит в выражение

cos г1) —

k

l

----sm“

 

____

cos {1

A

l /

i _ J ^

k

V

k2 sin2 d

Здесь — коэффициент отражения Френеля для магнитных ком­ понент.

Случай скользящего падения волны (# близко к л/2) требует применения формулы Вейля — Ван-дер-Поля.

§ 3.2. ИССЛЕДОВАНИЕ РАСПРОСТРАНЕНИЯ СВЕРХДЛИННЫХ РАДИОВОЛН ПОД ВОДОЙ

Рассмотрим теперь распространение электромагнитной волны, излучаемой диполем на границе раздела воздух—вода, в морской воде. В этой области г < 0 , и, для того чтобы с абсолютным возра­ станием z поле не увеличивалось, выбраны те знаки, которые при­ ведены во второй из формул (3.16). Преобразуя, имеем

Z = 2kl

и 0 (\r) е^§**-** zdX__

(3.22)

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

Z =

i k R

(3.23)

S.

5Т

Коэффициент ослабления S также выражается формулой Вейля— Ван-дер-Поля

S = 2 — 4xe-(*+T)‘ и+Jт ea2da,

— * со

но параметры х и т будут другими: з .

------ --------------Г - Я 1

 

___________________

kR

е 4

k

Г .

k2

sin2 ф

'

sin d

' k0 ] /

~~

*2

S1

 

 

3

.

(3.24)

 

 

ти

 

I /

4

 

 

г ‘

 

 

V

2 '

sin -O' COS Ф.

 

Коэффициент ослабления

5

для

морской

воды гораздо меньше

(т. е. ослабление значительно больше), чем для воздушной среды. Кроме того, вследствие большой мнимой части волнового числа k существенную роль играет экспоненциальное поглощение (благодаря

множителю eikR). Ослабление

амплитуды в е раз

достигается при

толщине скин-слоя б, определяемой формулой

 

 

_

 

S =

2 у,

(3.25)

где у — поправочный коэффициент (близкий к единице), зависящий от расходимости волн и коэффициента ослабления S.

При описании распространения радиоволн над морем необходимо учитывать сферичность Земли, в особенности для сверхдлинных радиоволн. Пусть имеется диполь, ориентированный по направлению земного радиуса. В данном случае единственной компонентой вектора Герца будет Zr. При этом сферическая система координат г, ф, 0 имеет начало в центре земного шара с полярной осью, проходящей через указанный вертикальный диполь. Вектор Герца уже не будет удовлетворять уравнению Гельмгольца (3.10). Но если ввести функ­ цию Герца U такую, что

 

 

Zr = Ur,

(3.26)

то

U удовлетворяет уравнению

 

 

 

A U + k2U = 0.

 

Решение этого

уравнения записывается в виде

 

 

 

ikR

(3.27)

 

 

U ^ ^ - S ,

где

R — У г2 +

а2— 2arcos'&, причем а есть радиус

Земли.

 

Когда величина ka велика и велико абсолютное значение диэлек­

трической проницаемости, функция ослабления будет медленно изменяться на протяжении одной длины волны. Это следует из того, что в обширной области формула Вейля — Ван-дер-Поля

58

переходит в выражение S = 1 + £, где £ — коэффициент Френеля. Однако для очень длинных волн таким приближением ограничиться нельзя.

В. А. Фоком задача была решена с учетом сферичности Земли. В этом случае функция Герца равна

 

 

 

 

p ika■0'

У-

9)-

(3.28)

 

 

 

 

 

где

 

 

 

 

 

 

 

 

V(x,

 

.

Я

 

elx1w (t — у) dt.

 

 

у, q) = e

 

 

 

 

 

 

 

 

w ' (t ) — qw (t )

 

Здесь интегрирование производится

по

замкнутому

контуру С

в комплексной плоскости t,

охватывающему корни знаменателя

в подынтегральном выражении; х и

у — безразмерные

величины,

через которые г

и О выражаются следующим образом:

 

 

 

 

r= = a{ l + T v ) ’‘

* =

(3-29>

,

 

1

 

 

 

 

( ka\ 3

; Я— параметр, зависящий от кривизны земной

причем А =

(-тг)

 

поверхности в данной точке.

Функция w (/), введенная В. А. Фоком, выражается через функ­

цию Ганкеля первого рода порядка 1/3 по формуле

 

w(t) = e 3

(3.30)

Далее, теория дает, что в области тени коэффициент ослабления быстро убывает, а в области полутени выражается через функции Эйри [37]. Для сверхдлинных волн область тени отодвигается по мере увеличения длины волны и может быть сведена на нет (обегание радиоволны вокруг земного шара).

Ограничим задачу тем кругом условий, которые имеют отношение к измерениям в области низкочастотных электромагнитных колеба­ ний. Для стандартной звуковой частоты 1000 Гц длина электро­ магнитной волны в воздухе равна 300 км, и, таким образом, можноограничиться рассмотрением поля на расстоянии, малом по сравне­ нию с длиной электромагнитной волны в воздухе или сравнимом с ней, т. е. ограничиться ближним полем. Специальные вопросы дальнего приема к указанной области измерений не относятся. Следовательно, имеются основания пренебречь сферичностью Земли и рассматривать распространение радиоволн над поверхностью моря как над плоскостью. При этом, учитывая скользящее движение волны, нужно исходить из теории Зоммерфельда.

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

59

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