Методы и средства передачи информации (Лекция №9)
.pdf2 A −ε a µ a |
∂ 2 A |
−grad div A −ε a µ a |
∂ |
grad ϕ= −µ a J , |
∂t 2 |
|
|||
|
|
∂t |
а сгруппировав члены с учетом независимости производных по времени и по координатам, получим:
|
2 |
|
∂ 2 A |
|
∂ϕ |
|
|
||
|
|
A −ε a µ a |
|
|
−grad div A +ε a µ a |
|
|
= −µ a J. |
(9.25) |
|
∂t |
2 |
∂t |
||||||
|
|
|
|
|
|
|
Учитывая, что вектор А недоопределен, можно положить, что
∂ϕ |
|
div A = −ε a µ a ∂t |
(9.26) |
так как rot rot А при этом не меняется. Это действие называется калибровкой векторного потенциала. Этот вариант калибровки называется калибровкой Лоренца и удобен при расчете переменных полей. С учетом (2.26) уравнение (2.25) примет вид:
2 A −ε a µ a |
∂ 2 A |
= −µ a J. |
(9.27) |
|
∂t 2 |
|
|
Вообще говоря, полученного уравнения достаточно для достижения поставленной цели, а именно, для установления связи между векторным потенциалом и источником поля − сторонним током проводимости. Однако более простой вид связи между вектором А и током J следует из анализа уравнения (9.27). Для этого проведем дополнительные действия.
Одновременно, с учетом представления div grad ϕ = 2 ϕ и принятого вы-
ражения калибровки (9.26) из уравнения (9.24), поменяв местами порядок дифференцирования, для скалярного потенциала получим соотношение:
∂ |
−ε a µ a |
∂ϕ |
+ |
2 |
ϕ = − |
ρ |
, |
||
|
|
|
|
|
|
||||
|
∂t |
|
ε a |
||||||
∂t |
|
|
|
|
|
|
из которого следует уравнение для скалярного потенциала:
2 ϕ−ε a µ a |
∂ 2 ϕ |
= − |
ρ |
. |
(9.28) |
∂t 2 |
|
||||
|
|
ε a |
|
11
Сравнивая уравнения (9.27) |
и (9.28) видим, что |
скалярная функция |
ϕ( x , y , z ,t )и проекции векторной |
функции A( x , y , z ,t ) |
на координатные оси |
х, у, z удовлетворяют уравнению одного и того же типа − уравнению Даламбера:
2 Ф−ε a µ a |
∂ 2 Ф |
= −η. |
(9.29) |
|
∂t 2 |
|
|
Уравнение (9.29) относится к уравнению гиперболического типа и называется волновым. Решением уравнения (9.29) являются функции бегущей волны. Это проще всего можно показать для однородного волнового уравнения:
2 Ф−ε a µ a |
∂ 2 Ф |
= 0 . |
(9.30) |
|
∂t 2 |
|
|
Это уравнение является частным случаем уравнения Даламбера для области пространства свободной от сторонних источников поля.
Ограничимся рассмотрением на примере функции прямоугольной декартовой системы координат Ф( x , y , z ,t )= Ф( z ,t ). При этом уравнение (9.30) при-
мет вид одномерного однородного уравнения в частных производных
∂ 2 Ф |
−ε a µ a |
∂ 2 Ф |
= 0 , |
∂z 2 |
|
∂t 2 |
|
решением которого (что легко проверить подстановкой этого решения в уравнение) имеет вид:
|
|
z |
|
|
z |
|
|
||
Φ = Φ1 |
t − |
|
|
+Φ 2 |
t + |
|
|
, |
(9.31) |
|
|
||||||||
|
|
v |
|
|
v |
|
|
где Φ1 и Φ 2 – функции бегущих плоских (т.е., с плоским фронтом) волн, у кото-
рых v = |
1 |
– фазовая скорость. Причем волна Φ1 распространяется c фа- |
ε a µ a |
зовой скоростью v вдоль оси z, а Φ 2 – в обратном направлении. Заметим, что аналогично записываются прямые и отраженные (обратные) волны в длинной линии.
В отсутствии направляющей структуры в свободном пространстве в сферической системе координат для полей точечного источника (в пространстве вне
12
источника) волновое уравнение (9.30) с учетом коэффициентов криволинейности осей имеет вид:
∂ 2 Φ |
+ |
2 ∂Φ |
−εµ∂ 2 Φ = 0, |
(9.32) |
|||||
∂r 2 |
|
r ∂r |
|
|
∂t 2 |
|
|
||
14243 |
|
|
|
|
|
|
|||
|
2 Φ |
|
|
|
|
|
|
||
а с применением подстановки |
ζ = Фr запишется в виде: |
||||||||
|
∂ 2 ζ |
−ε a µ a |
∂ |
2 ζ |
= 0 . |
|
|
||
|
∂r 2 |
∂t 2 |
|
|
|||||
|
|
|
|
|
|
|
|||
Решение волнового уравнения для новой переменной ζ записывается как |
|||||||||
|
ζ = ζ1 |
|
r |
|
|
r |
, |
||
|
t − |
|
|
+ ζ 2 t + |
|
||||
|
|
|
|
v |
|
|
v |
|
откуда с обратной заменой переменных получим решение уравнения (9.32) в ви-
|
|
|
|
|
Φ1 |
|
r |
|
r |
|
|
|||
|
|
|
|
|
t − |
Φ 2 |
t + |
|
|
|
||||
де: |
|
|
|
Φ = |
|
|
v |
+ |
|
|
v |
|
|
|
|
|
r |
|
|
r |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
где |
Φ1 |
и |
Φ 2 |
– сферические волны (соответственно прямая и обратная или |
||||||||||
r |
|
|||||||||||||
|
|
r |
|
|
|
|
|
|
|
|
|
|||
отраженная). |
|
|
|
|
|
|
|
|
|
|||||
|
Вид функции Φ1 определяется источником, а обратной функции |
Φ 2 – при- |
||||||||||||
чиной её возникновения (отражения). |
|
|
|
|
|
|
|
|
||||||
|
Для точечного источника вид Φ1 можно установить, рассматривая решение |
|||||||||||||
задачи для точечного электрического заряда q, находящегося в точке r= 0. |
||||||||||||||
|
Уравнение Даламбера (9.32) при постоянном заряде q (т.е. при |
∂ |
→0 ) в |
|||||||||||
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
∂t |
окружающем заряд пространстве переходит в уравнение Лапласа 2 ϕ = ∆ϕ = 0
( ∆ – символ скалярного лапласиана) для потенциала ϕ( r ) (в силу центральной
симметрии зависящего только от расстояния r до точки наблюдения) электростатического поля Е, которое в сферической системе координат имеет вид:
1 |
|
∂ |
|
2 |
∂ϕ |
|
|
|
|
|
|
r |
|
|
|
= 0. |
(9.33) |
r 2 |
|
|
|
|||||
|
∂r |
|
∂r |
|
|
|
13
Решение уравнения (9.33), полагая радиус r величиной конечной, сводится
кинтегрированию уравнения
∂r 2 ∂ϕ = 0.
∂r ∂r
Откуда получим: r 2 |
∂ϕ |
=C 1 , где C 1 |
– постоянная интегрирования, и затем: |
||||||
|
|||||||||
|
|
|
|
|
|
|
∂r |
|
|
|
∂ϕ |
= |
C 1 |
и ϕ = − |
C 1 |
+ C 2 . |
|
||
|
∂r |
r 2 |
|
|
|||||
|
|
|
r |
|
|||||
|
|
|
Постоянные интегрирования C 1 |
и C 2 дифференциального уравнения по- |
лучают из известных значений искомой функции. Такие данные в теории поля называют граничными условиями. Это могут быть данные о значениях поля (или потенциала) в точках или на поверхностях анализируемой области пространства, например в бесконечно удаленных точках. Так как в нашем случае постоянных интегрирования две, то для составления двух уравнений необходимо иметь данные как минимум в двух точках анализируемого пространства, либо данные о функции и её производной в одной точке. В нашем случае можно ука-
зать только одно известное значение ϕ( r ) = 0 при r → ∞ (стремящемся к беско-
нечности), полагая физическое поле бесконечно удаленного источника убывающим. При этом C 2 = 0. Для определения постоянной C 1 прямых данных нет.
Поэтому воспользуемся косвенными данными.
Рассмотрим решение задачи исходя из теоремы Гаусса. Как мы знаем (см. лекцию №8), электростатическое поле точечного источника в свободном про-
странстве |
ϕ( r ) = |
|
q |
, поэтому, формально заменив постоянный точечный |
|||
|
πε a r |
||||||
|
4 |
|
|
|
|
||
заряд q на переменный q (t), запишем ϕ( r ,t ) = |
q ( t ) |
. |
(9.34) |
||||
|
|||||||
|
|
|
|
|
4 πε a r |
|
Теперь, зная, что решение волнового уравнения выражается в функции бегущих волн, если предположить, что при фазовой скорости v ≠ ∞,
14
|
|
r |
|
|
|
|
q t − |
v |
|
|
|
ϕ( r ,t ) = |
|
|
, |
(9.35) |
|
4 πε a r |
|
||||
|
|
|
|
то при v →∞, выражение (9.35) перейдет в (9.34), то есть решение (9.34) явля-
ется предельным переходом для более общего решения (9.35).
Заметим, что мы лишь показали, какой вид функции логично ожидать.
Из найденного выражения (9.35) для потенциала точечного источника с применением метода суперпозиции следует решение для потенциала в точке наблюдения М от распределенного в объеме V заряда ρ( x , y , z ,t ) :
ϕ ( r ,t ) = |
1 |
|
∫ |
ρ(t − r |
v |
) |
|
|
|
|
|
dv. |
(9.36) |
||||
4 πε |
a |
r |
|
|
||||
|
|
|
V |
|
|
|
|
|
Из выражения (9.36) следует, что в точке наблюдения M в момент времени t суммируются парциальные потенциалы от значений ρk(t,r), взятых в разные моменты времени t'= t – r'/v, что иллюстрирует рис. 9.3. Этот факт характеризует понятие «запаздывающие потенциалы».
r' |
М |
ρ1 |
|
|
r" |
ρ2
Рисунок 9.3 − Иллюстрация к понятию «запаздывающие потенциалы Аналогично, рассмотрев выражение для векторного потенциала А элемента
постоянного тока Il (закон Био-Савара – см. лекцию №8, c. 17) и формально перейдя от уравнения Пуассона для векторного потенциала постоянного элемента тока Il
2 A = −µ a J |
(9.37) |
к уравнению Даламбера, можно записать для векторного потенциала выражение
15
A (r ,t )= |
µ a |
∫ |
J (t −r v ) |
dV . |
(9.38) |
4 π |
r |
||||
|
|
V |
|
||
Поэтому, в отличие от скалярного ϕ и векторного А потенциалов электро- |
статического ( ∂/ ∂t →0 и J = 0) и стационарного ( ∂/ ∂t →0 и J ≠ 0) полей, потен-
циалы электромагнитного поля называют «запаздывающими» – т. е. учитывающими время распространения волны.
С точки зрения вопроса определения структуры магнитного поля, который мы анализируем, понятия «запаздывающие потенциалы», возникшее, что называется, попутно не имеет прямого отношения. Однако с мировоззренческой точки зрения, полученный результат весьма важен, так как подтверждает уже известное нам понятие эффекта «близкодействия» с которым мы познакомились при рассмотрении распространения волн в длинной линии.
Возвращаясь к рассматриваемому вопросу, запишем формулы (9.36) и (9.38) для линейных зарядов τ( r ,t ) и токов i(r, t), характеризующих решения за-
дач в двухпроводной линии. При этом формулы (9.36) и (9.38) принимают вид:
ϕ = |
1 |
|
|
∫ |
τ(t −r v ) |
dl , |
(9.39) |
|||
4 πε |
a |
r |
||||||||
|
|
|
|
|
|
l |
|
|||
A = |
µ a |
∫ |
i (t −r v ) |
dl . |
(9.40) |
|||||
4 π |
r |
|||||||||
|
|
|
|
l |
|
|
|
|
|
Заметим, что в отличие от скалярного выражения (9.39), выражение (9.40) соответствует векторной функции, поэтому распадается на три скалярных уравнения
Aξk |
= |
µ a |
|
i (t − r v ) |
dl |
ξk |
(9.41) |
|||
|
|
|
|
|
||||||
4 |
π ∫ |
r |
||||||||
|
|
|
|
|
||||||
|
|
|
|
l |
|
|
|
|
|
по координатным осям ξk. Но так как в двухпроводной линии ток совпадает пусть с осью z, то выражение для векторного потенциала ограничится только од-
ной координатной осью z, т.е. A = z 0 A z , где z 0 − орт продольной оси линии.
Выражение (9.39) окажется полезным при рассмотрении в дальнейшем изложении понятия перестановочной двойственности, применяемого для расчета
16
первичных параметров двусвязных систем. Выражение (9.41) является основой получения силовых линий магнитного поля в двусвязных структурах.
Объясняется это тем, что каждый из токов в проводниках двухпроводной линии обуславливает векторный потенциал, параллельный линиям тока и имеющий единственную продольную составляющую, например по оси z, не зависящую от координаты z. Очевидно, что результирующий векторный потенциал А=z0 Аz(х, у) будет равен алгебраической сумме векторных потенциалов отдельных токов, а результирующая магнитная индукция В =rotА определится дифференцированием:
|
|
|
|
x 0 y 0 |
z 0 |
|
|
x 0 y 0 z 0 |
|
|
∂A z |
|
∂A z |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
B = rot A = |
|
∂ |
|
|
∂ |
|
|
∂ |
|
= |
|
∂ ∂ ∂ |
|
= x 0 |
− y 0 |
. |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
∂x ∂y ∂z |
|
∂x ∂y ∂z |
∂y |
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
∂x |
||||||||||||||
|
|
|
|
A x |
A y |
A z |
|
|
|
0 0 A z |
|
|
|
|
|
|
|||||||||
Поэтому B x |
= |
∂A z |
; B y |
= − |
∂A z |
. |
|
|
|
|
|
|
|
|
|
|
|
(9.41) |
|||||||
∂y |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
∂x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Теперь запишем через векторный потенциал А уравнение линий магнитной индукции В. Для этого сначала запишем дифференциальное уравнение линий магнитной индукции, которое следует из определения вектора индукции поля как касательной к силовой линии (см. рис. 9.4).
В
r0
Рисунок 9.4 − К выводу уравнения силовой линии
Здесь вектор r 0 = x 0 dx + y 0 dy + z 0 dz , а вектор B = x 0 B x + y 0 B y + z 0 B z , а
условие колинеарности этих векторов выразится равенством векторного произ-
ведения [r 0 ,B ]= 0, которое в явном виде с учетом равенстваH z = 0 и плоскопа-
раллельной структуры поля запишется как
17
|
x 0 y 0 z 0 |
|
|
|
x 0 y 0 z 0 |
|
= z 0 B y dx − z 0 B x dy |
|
|
|
|
||||
0 = |
dx dy dz |
|
= |
|
dx dy 0 |
|
|
|
B x B y B z |
|
|
|
B x B y 0 |
|
|
или B y dx = B x dy . |
|
|
|
|
(9.42) |
Если теперь в дифференциальном уравнении (9.42) этих линий сделать подстановку из уравнения (9.41), то получим
Adxz dx + Adyz
Откуда Аz(х, у) = const и есть искомое уравнение силовой линии.
Теперь необходимо установить вид этой силовой линии, т.е. вид кривой Аz(х, у) = const в двухпроводной линии.
В общем случае на поверхности проводов Аz(х, у) ≠ const, так как магнит-
ные линии пересекают эту поверхность.
Ограничимся рассмотрением магнитных силовых линий вне проводов двухпроводной линии. Найдем вид функции Аz(х, у).
Исходим из того, что по закону полного тока на расстоянии r ≥ a , где а −
радиус провода, напряженность магнитного поля H α = 2 Iπr , а магнитная ин-
дукция B α |
= |
µ a I |
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
2 πr |
|
|
|
|
|||||
Векторный потенциал Аz(r) проще всего определить по формуле |
|
||||||||||
|
|
|
B α = rot α A = − |
dA z |
. |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
dr |
|
|
|
|
||
Таким образом, для внешнего магнитного поля одного провода векторный |
|||||||||||
потенциал |
|
|
|
|
|
|
|
|
|
|
|
|
|
A z = ∫− B α dr = − ∫ |
µ a I |
|
µ a I |
|
|||||
|
|
|
dr |
= − |
|
ln r + C . |
(9.43) |
||||
|
|
2 πr |
2 πr |
Для двух проводов векторные потенциал получим в виде суммы вкладов от каждого из проводов. При этом аналогично рассмотрению задачи для электростатического поля двухпроводной линии обозначим расстояния от «положитель-
18
ного» тока r + . Пусть это ток «уходящий за плоскость листа», соответствующий положительно заряженной оси (т.е. провод, к которому приложен положительный зажим источника напряжения). Расстояния от «отрицательного» тока (или обратного провода линии) r − . Тогда Суммарный векторный потенциал
|
|
+ |
|
− |
|
µ a I |
|
+ |
|
|
µ a I |
|
− |
|
|
+ |
|
− |
|
|||
A zΣ = |
A z |
+ A z |
= − |
|
|
ln r |
|
− |
− |
|
|
ln r |
|
|
+ C |
|
+ C |
|
= |
|||
2 |
π |
|
2 |
π |
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
= |
µ a I |
ln |
r − |
|
+ C Σ . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(9.44) |
2 πr |
r + |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Формула (9.44 повторяет по структуре выражение скалярного потенциала
электростатического поля двухпроводной линии. Также как и в случае скалярно-
го потенциала при равенстве r + = r − получаем |
ln1 = 0 , а значит приняв A zΣ=0 |
||||||||||
на оси у (т.е. при х = 0) получим константу C Σ = 0 |
и A zΣ |
= |
µ a I |
ln |
r − |
. |
(9.45) |
||||
|
|
|
|||||||||
|
|
|
|
|
2 π |
|
|
r + |
|
||
Итак, силовые линии магнитного поля A zΣ |
= |
µ a I |
ln |
r − |
=const повторяют |
||||||
2 π |
r + |
||||||||||
|
|
|
|
|
|
|
|
|
по записи эквипотенциальные кривые электростатического поля, имеющие вид
ϕ= |
|
τ |
ln |
r − |
. |
|
|
|
|||
2 |
πε a |
|
r + |
Этим доказано, что силовые линии магнитного поля двухпроводной линии
представляют собой окружности, центры которых определены аналогично центрам эквипотенциальных окружностей электростатического поля в плоскости поперечного сечения двухпроводной линии.
В итоге можно заключить, что линии электрического поля, если пренеб-
речь падением напряжения на проводах, в поперечном сечении линии также образуют окружности, нормальные к линиям магнитного поля. Картина поля пока-
зана на рис. 9.5. Заметим, что падение напряжения на проводах (за счет погон-
ного сопротивления R =1/σS ) приводит к изменению (уменьшению в направлении от генератора к нагрузке) напряжения в поперечных сечениях линии (это явление называют потерями напряжения). Это искажает напряженность электриче-
19
ского поля в линии за счет появления ненулевой тангенциальной его составляющей (поле теоретически не идеально плоскопараллельное). Однако тангенциальная составляющая электрического поля обычно ничтожно мала по
сравнению с нормальной составляющей. Так, если плотность тока в проводе J
= 1 А/мм², то падение напряжения на одном сантиметре медного провода
Et = |
J |
= |
|
1 |
≈1,76 10−4 B/м |
σ |
|
57 104 |
|||
|
0,01 |
|
|||
в то время как нормальная |
составляющая напряженности электрического поля |
на поверхности провода при таком значении плотности постоянного тока (определяемой приложенным к линии напряжении) составляет примерно 1 В/см и выше.
Рисунок 9.5 – Структура полей в двухпроводной линии
Исходя из плоскопараллельной модели поля в такой линии, найдем выра-
жение для мощности, передаваемой через её поперечное сечение. Для этого за-
пишем поток вектора Пойнтинга через дуговые четырехугольники, образован-
ные в плоскости поперечного сечения пересечением силовых линий электрического и магнитного полей (рис. 9.3). Углы при вершинах этих четырехугольни-
ков прямые (так как это углы между силовыми линиями электрического поля и
эквипотенциальными линиями. Площадь элементарного четырехугольника ds=da·db, поэтому поток плотности вектора Пойнтинга через него
20