Двумерная модель радиационной электромагнитной наводки в кабельной линии
|
Рассмотрим участок кабельной линии (см. рис. 2).
Запишем второе уравнение Максвелла в интегральной форме [2]:
(1)
где S - поверхность, ограниченная
замкнутым контуром С;
Es
- касательная составляющая вектора электрического поля;
Bn - нормальная составляющая вектора магнитной индукции;
сэл - скорость распространения электромагнитных
возмущений.
Используем его для замкнутого
прямоугольного контура abcda. Циркуляция электрического поля
вдоль сторон замкнутого контура равна:
(2)
где Ezжил, Ezопл - касательные составляющие электрического поля на жиле и
оплётке соответственно;
U(z) -
напряжение между жилой и оплёткой кабеля в сечении z
Полагая, что Ezжил=RжилI и Ezопл=-RоплI, получаем:
(3)
где R=(Rжил+Rопл) - погонный
внутренний импеданс кабеля.
Поток вектора магнитной
индукции через контур abcda можно представить в виде:
(4)
где - магнитный поток
между жилой и оплёткой на единицу длины кабеля.
Таким образом, уравнение (1)
для контура abcda принимает вид:
(5)
Запишем уравнение
непрерывности в интегральной форме [2]:
(6)
где r - плотность
заряда в объёме, ограниченном замкнутой поверхностью S;
jn -
нормальная составляющая вектора плотности тока на поверхности S.
Применим
его к отрезку z, z+dz жилы (см.
рис. 3). Интеграл определяет заряд на отрезке жилы и поэтому равен произведению
погонной плотности заряда qжил на элемент
длины dz. Поверхностный интеграл разложим на сумму интегралов по
поперечным сечениям z, z+dz и боковой
поверхности провода S0 на участке
(z, z+dz):
(7)
где s0, e0 - проводимость и
диэлектрическая постоянная изоляции между жилой и оплёткой;
Dn - нормальная составляющая вектора электрической индукции.
Интеграл с Dn можно преобразовать, используя теорему Гаусса:
(8)
где sжил, eжил - проводимость и диэлектрическая постоянная жилы
кабеля.
Таким образом, получаем:
(9)
Учитывая, что s0 << sжил (проводимость диэлектрической изоляции гораздо меньше
проводимости металлической жилы), получаем:
(10)
Проведя аналогичные
рассуждения для оплётки, получим:
(11)
До этого момента для
построения уравнений использовались только уравнения Максвелла, поэтому все
полученные соотношения носят фундаментальный характер. Теперь же сделаем
предположение о локальном характере электромагнитных полей внутри кабеля, т.е.
-
магнитное поле в поперечном сечении z=const кабеля определяется исключительно
током I, текущем по жиле и оплётке в этом
сечении;
-
электрическое поле в поперечном сечении z=const кабеля определяется исключительно
погонной плотностью зарядов на жиле и оплётке кабеля (qжил, qопл), а также
распределением электрического заряда внутри диэлектрической изоляции d(t,r,z).
Предположение о локальности
магнитного поля позволяет легко связать магнитный поток между жилой и оплёткой
на единицу длины кабеля y с током
текущем в кабеле I:
(12)
где - погонная
индуктивность кабеля;
cэл = 2.9979´1010
см/c – скорость распространения
электромагнитных возмущений в вакууме;
m0 –
магнитная проницаемость диэлектрической изоляции.
Используя эту связь,
получаем:
(13)
Предположение о локальности
электрического поля используем для получения связи между параметрами погонной
плотности зарядов на жиле и оплётке кабеля (qжил, qопл) и напряжением
U (вполне очевидно, что это соотношение
будет зависеть и от распределения электрического заряда внутри диэлектрической
изоляции d(t,r,z)).
Пусть в момент времени t в поперечном сечении z=const кабеля сложилось следующее распределение электрических
зарядов (см. рис. 4):
-
qжил - погонный
заряд жилы кабеля;
-
qопл - погонный
заряд оплётки кабеля;
-
d(t,r,z) - распределения электрического заряда внутри диэлектрической
изоляции.
Рассматривая
кабель как бесконечный коаксиальный цилиндрический конденсатор, обкладками
которого являются жила (r1 = rжил) и оплётка
(r2 = rопл), с
объёмно распределённым электрическим зарядом в диэлектрической изоляции получим
выражение для напряжения между обкладками. Учитывая, что потенциал бесконечно
длинной тонкой цилиндрической поверхности радиуса a с погонной плотностью электрического заряда q равен:
(14)
получаем:
(15)
В выражении (15) в фигурных
скобках представлены алгебраические суммы потенциалов всех зарядов системы (в
жиле, оплётке и диэлектрической изоляции) соответственно на жиле и оплётке.
Считая d не
зависящей от r, выражение
(15) после интегрирования приведётся к виду:
(16)
Выразив из (16) qжил и подставив в (10), окончательно получим:
(17)
где - погонная ёмкость
кабельной линии;
- погонный коэффициент
утечки;
;
.
Таким образом, проводя
совместное интегрирование уравнений (13) и (17) с использованием граничных
условий:
,
(18)
где L – длина
кабельной линии;
R0(t) – нагрузка на левом конце
кабельной линии;
RL(t) – нагрузка на правом конце
кабельной линии,
мы получаем параметры
паразитной наводки, возбуждаемой в кабельной линии при импульсном радиационном
воздействии.
Для совместного
интегрирования уравнений (13), (17) с граничными условиями (18) и начальными
условиями:
(19)
используем
метод характеристик [3], полагая, что Lкаб = const, R = const и C = const.
Определим уравнений
характеристик как
(20)
где a = 45°;
v – скорость распространения электромагнитной волны.
Система уравнений (13), (17)
примет вид:
(21)
где .
Учитывая, что (r - волновое
сопротивление линии) и , преобразуем систему (21) к виду:
(22)
Введём новые функции f и c:
(23)
Тогда система (22) примет
вид:
(24)
где ,
.
Для решения системы (24)
применим неявную конечно-разностную схему Кранка-Николсона [4], шаблон которой
представлен на рис. 5.
Имеем:
(25)
где (Dz – шаг по пространству);
- шаг по времени;
- значения
соответствующих функций в момент времени в точке с координатой .
Систему
(25) можно представить в следующем виде:
(26)
где
.
Разрешая систему (26)
относительно , получаем:
(27)
Из соотношения (18) следует,
что на левом и правом концах линии в любой момент времени выполняются следующие
соотношения:
(28)
где ;
.
Используя соотношения (26) и
(28), получаем выражения для искомых функций на концах линии:
(29а)
(29б)
Таким образом, используя
соотношения (27), (29а), (29б), (23) и учитывая, что для всех n, можно определить значения U и I в любой
точке линии в любой момент времени. Условие устойчивости для данной
конечно-разностной схемы выглядит следующим образом:
,
(30)
где - шаг по пространству;
- шаг по времени;
v – скорость распространения электромагнитной волны в
линии .
3. Задание
исходных данных для расчета
Для расчета параметров
электромагнитной наводки в коаксиальной кабельной линии при воздействии
импульса ионизирующего излучения необходимо задание следующих исходных данных.
1.
Параметры и характеристики кабельной линии:
§
радиус жилы (см);
§
радиус оплётки (см);
§
длина (см);
§
погонное сопротивление жилы и оплётки (с/см);
§
электрическая проводимость материала изоляции (1/с);
§
относительная диэлектрическая проницаемость изоляции;
§
относительная магнитная проницаемость изоляции.
2.
Параметры и характеристики воздействия:
§
амплитудное значение погонного радиационно-стороннего
тока, натекающего на жилу кабеля ;
§
объёмная плотность заряда в диэлектрической изоляции
между оплёткой и жилой ;
§
временная зависимость погонного радиационно-стороннего
тока, натекающего на жилу кабеля.
3.
Требуемые параметры расчета:
§
количество ячеек по длине кабеля;
§
временной шаг записи результатов расчета на диск;
§
шаг записи результатов расчета по длине кабеля.
Тестирование методики
производилось на решении модельных задач электродинамики. Были проведены
расчёты задач зарядки и разрядки кабельных линий, результаты сопоставлялись с
аналитическими решениями.
Зарядка
кабельной линии неравномерным потоком ионизирующего излучения. Кабель
(длина
,
где e = 8.854´10-12 Ф/м;
eo =2.5;
rопл =
rжил =
;
Lкаб =
Результаты расчёта показали,
что течение процесса зарядки кабеля полностью соответствует классическим
представлениям [5], выход на стационарный режим завершает процесс перезарядки
кабеля, в результате которой напряжение между жилой и оплёткой становится
одинаковым по всей длине кабеля. Присутствуют колебания, обусловленные резонансом,
вызванным переотражением электромагнитной волны от разомкнутых концов кабеля
(полупериод колебаний соответствует времени прохождения волны по кабельной линии).
Разрядка
кабельной линии на согласованную нагрузку. Кабель (длина
1.
А.М. Волощенко, С. В. Гуков, В. В. Шаховский
«Исследование радиационного воздействия на коаксиальный антенный кабель»,
Сборник докладов Всероссийской конференции «Радиационная стойкость электронных
систем – СТОЙКОСТЬ - 2002», г. Лыткарино Московской обл.,
2.
Л. Д. Ландау, Е.М. Лившиц «Теоретическая физика. Теория
поля» - М.: «Наука»,
3.
А. Н. Тихонов, А. А. Самарский «Уравнения математической
физики» - М.: «Наука»,
4.
Д. Поттер «Вычислительные методы в физике» М.: «Мир»,
5.
В. И. Вольман, Ю. В. Пименов «Техническая
электродинамика» - М.: «Связь»,
6.
Г. В. Зевеке, П. А. Ионкин, А. В. Нетушил, С. В. Страхов
«Основы теории цепей» - М.: «Энергоатомиздат»,