Министерство науки и высшего образования Российской Федерации
Санкт-Петербургский политехнический университет Петра Великого
Институт физики, нанотехнологий и телекоммуникаций
Работа допущена к защите
Директор ВШПФиКТ
_______________ Е.Н. Величко
«___» _________ 2020 г.
ВЫПУСКНАЯ КВАЛИФИКАЦИОННАЯ МАГИСТЕРСКАЯ
ДИССЕРТАЦИЯ
ИССЛЕДОВАНИЕ АМПЛИТУДЫ И ФОРМЫ ТОКОВ В
ЗАЗЕМЛИТЕЛЯХ ПРОМЫШЛЕННЫХ МОЛНИЕОТВОДОВ
по направлению подготовки 16.04.01 Техническая физика
Профиль 16.04.01_05 Радиофизика и электроника
Выполнил
Студент гр. 3441601/80501
В.В. Курикова
Научный руководитель
доцент, к.ф.-м.н.
К.В. Грешневиков
Научный консультант
старший преподаватель
Г.П. Жабко
Консультант
по нормоконтролю
Е.В. Савченко
Санкт-Петербург
2020
РЕФЕРАТ
На 54 с., 15 рисунков,
КЛЮЧЕВЫЕ
СЛОВА:
ЗАЗЕМЛИТЕЛЬ,
МОЛНИЕОТВОД,
НЕСИММЕТРИЧНЫЙ ВИБРАТОР, КОЛЕБАТЕЛЬНЫЙ ПРОЦЕСС, ТОК,
ИМПУЛЬ НАПРЯЖЕНИЯ, ПРОВОДИМОСТЬ
Тема выпускной квалификационной работы: «Исследование амплитуды и
формы токов в заземлителях промышленных молниеотводов».
Данная работа посвящена колебательному процессу. Главная задача
состояла в определении при каких условиях могут возникать колебательные
движения. Задачи, которые решались в ходе исследования:
1. Обзор современных методов расчета заземлителей.
2. Разработка
полупроводящую
цилиндрического
среду.
заземлителя,
Эквивалентная
погруженного
электрическая
схема
в
с
сосредоточенными параметрами.
3. Исследование процессов распространения волн тока в заземлителе как
в несимметричном вибраторе. Разработка электродинамической модели и
алгоритма расчета.
4. Переходные процессы при импульсном воздействии. Расчет уровня и
формы токов.
Были получены аналитически формулы тока, на основе представления
заземлителя как электрического вибратора, погруженного в среду.
Показано, что колебательные процессы можно представить как аналог
колебательного контура.
Все это полезно для тех, кто разрабатывает заземлители и занимается
созданием молниеотводов. Это дает возможность воспользоваться результатом
формул
для
исследования
заземлителей, земли и т.д.
тока
для
любых
импульсов,
параметров
ABSTRACT
54 pages, 15 pictures
KEYWORDS: GROUNDING CONDUCTOR, LIGHTNING COLLECTOR,
UNSYMMETRIC
VIBRATOR,
OSCILLATORY
PROCESS,
CURRENT,
VOLTAGE IMPEDANCE, CONDUCTIVITY
The subject of the graduate qualification work is “Study of amplitude and
shape of currents in earth electrodes of industrial lightning leads”.
The given work is devoted to oscillatory process. The main task was to
determine under what conditions oscillating movements could occur. Objectives of
the study:
1. Overview of modern methods of earth electrode calculation.
2. Development of a cylindrical earth electrode immersed in a semi-conductive
medium. Equivalent electrical circuit with concentrated parameters.
3. Study of current wave propagation processes in the earthing electrode as in
an asymmetric vibrator. Development of electrodynamic model and calculation
algorithm.
4. Transient processes under pulse action. Calculation of current level and
shape.
Analytically, current formulas were obtained based on the representation of the
earthing electrode as an electric vibrator immersed in the medium.
It is shown that oscillatory processes can be represented as analogue of
oscillatory circuit.
All this is useful for those who develop earth electrodes and create lightning
leads. This makes it possible to use the result of the formulas to examine the current
for any pulses, earth electrode parameters, earth, etc.
4
СОДЕРЖАНИЕ
Введение .................................................................................................................. 6
Глава 1. Обзор литературы..................................................................................... 8
1.1 Общие сведения ................................................................................................ 8
1.2 Молниезащита ................................................................................................... 8
1.3 Параметры и характеристики промышленных молниеотводов ................... 9
1.4 Электрические характеристики почв ............................................................ 11
1.5 Обзор современных методов расчета заземлителей .................................... 11
Глава 2. Ток в заземлителе ................................................................................... 19
2.1 Общие положения ........................................................................................... 19
2.1.1 Зависимость от времени .............................................................................. 19
2.1.2 Обозначения ................................................................................................. 20
2.2 Задача электростатики .................................................................................... 21
2.3 Задача постоянных токов ............................................................................... 22
2.4 Эквивалентная схема ...................................................................................... 25
2.5 Антенная задача .............................................................................................. 27
2.6 Импульсный ток .............................................................................................. 29
2.6.1 Импульсное напряжение ............................................................................. 29
2.6.2 Спектры сигналов. Восстановление функции i(t) ..................................... 33
2.7 Колебательные процессы ............................................................................... 33
Глава 3. Примеры расчетов .................................................................................. 36
3.1 Программа........................................................................................................ 36
3.2 Построение графиков ..................................................................................... 37
3.2.1 График импульса напряжения .................................................................... 37
3.2.2 Графики тока ................................................................................................ 38
3.3 О характере Z(w) ............................................................................................. 45
Заключение ............................................................................................................ 46
Список литературы ............................................................................................... 48
5
Список сокращений
СИ – система исчислений;
ЗУ – заземляющее устройство.
6
ВВЕДЕНИЕ
Основной темой будем считать заземлители.
Заземлители
представляют
из
себя
проводники,
непосредственно
соединенные с землей. Это может быть один или группа проводников.
Для защиты от удара молнии на здания и сооружения устанавливают
специальные устройства – молниеотводы.
Отличие данной работы в применении теории антенн к полупроводящей
среде для расчета токов и параметров заземлителей молниеотводов. Обычно
рассматривают в статике. Мы рассматриваем колебательный - более сложный
случай, как несимметричный вибратор.
Все это полезно для тех, кто разрабатывает заземлители и занимается
созданием
молниеотводов.
возникающими
Они
колебательными
обеспокоены
процессами.
большими
Это
дает
токами
и
возможность
воспользоваться результатом формул для исследования тока для любых
импульсов, параметров заземлителей, земли и т.д.
Наша главная задача состоит в том, что колебательный процесс
возникает. Мы рассматриваем вибратор, как длинную линию. Привычный
вариант – это когда волна тока доходит до конца, и пока дойдет весь ток
разойдется. А если вокруг цилиндра почва плохо проводящая, ток не успевает
выходить, доходит до конца и опять обратно. В итоге ток будет ходить «тудасюда» и возникают переколебания.
Мы аппарат антенн используем и исследуем при каких параметрах почвы,
размерах возникает колебательный процесс.
Объектом исследования в нашем случае – процессы тока в заземлителе.
Мы используем те подходы, формулы, которые используются в теории
антенн.
Отличие состоит в том, что размеры другие, частоты. Излучающая часть
антенны окружена воздухом без потерь. Нас интересует сколько она
электромагнитных волн излучит в пространстве. У нас нет никаких излучений,
7
она окружена проводящей глинной. Токи стекают в эту глины, почву или песок
и никакого излучения тут нет.
Наша задача состоит в том, чтобы токи уходили в землю. По сравнению с
проводимостью земли, наш металл, не зависимо от материала из которого он
изготовлен, будет считаться идеально проводящим. Мы будем говорить о
потерях окружающей среды.
Цель работы: определить при каких условиях могут возникать
колебательные движения.
Задачи работы:
1) Обзор современных методов расчета заземлителей.
2) Разработка
полупроводящую
цилиндрического
среду.
заземлителя,
Эквивалентная
погруженного
электрическая
схема
в
с
сосредоточенными параметрами.
3) Исследование процессов распространения волн тока в заземлителе как
в несимметричном вибраторе. Разработка электродинамической модели и
алгоритма расчета.
4) Переходные процессы при импульсном воздействии. Расчет уровня и
формы токов.
Новизна работы:
1. Впервые в отличие от обычного подхода получаем колебательный
процесс тока.
2. Получаем для любого заземлителя, для любых параметров расчет тока.
3. Получаем возможность исследовать ток для любых импульсов.
8
ГЛАВА 1. ОБЗОР ЛИТЕРАТУРЫ
1.1 Общие сведения
Молния – это капризное, случайное и непредсказуемое событие. Его
физические характеристики включают уровни тока, иногда превышающие
400 кА, температуры до 50000 градусов по Фаренгейту и скорости,
приближающиеся к одной трети скорости света.
Молния – это разряд электрического поля, которое образовалось между
облаком и землей и другими облаками. Это происходит с нарушением
атмосферы в локальной области.
Удары молнии могут вызвать большие разности потенциалов в
электрических системах, которые распределены вне одной защищенной
структуры или здании. Большие величины тока могут течь в земле и
металлическом материале в попытке выровнять потенциал. Часто в проводе
между двумя изолированными точками протекают большие импульсные токи.
Кроме того, удар создает большое электромагнитное поле, которое может
вызывать помехи в электронных и коммуникационных сигналах.
1.2 Молниезащита
Молния встречается на разных уровнях. Система должна быть
спроектирована так, чтобы справляться с типичным и большим ударом в это
области. Удары не могут быть остановлены, но энергия может быть отклонена
контролируемым образом. Попытка защитить от прямых ударов может повлечь
за собой чрезмерные расходы.
Молниезащита включает в себя обеспечение контролируемого пути с
низким импедансом к земле для энергии, а не позволяет ей выбирать
случайный путь с более высоким импедансом через структуру.
Во время удара наземная система будет насыщаться и иметь повышенный
потенциал относительно окружающей территории. Повышенный потенциал
9
будет сохраняться до тех пор, пока переходный процесс не распространится
через систему в землю.
Оборудование в области повышенных потенциалов сконфигурировано
для компенсации роста потенциала. Если все электрическое оборудование
связано вместе с одной и той же эффективной системой заземления, не будет
разницы в потенциале между точками в системе. Однако устройства,
подключенные к линиям, выходящим за пределы области равного потенциалы,
могут быть повреждены.
Несколько предметов влияют на вероятность удара. Частота удара
молнии возрастает вместе с квадратом высоты над средней местностью. Это
делает изолированные структуры особенно уязвимыми. Острые точки
увеличивают
вероятность
удара.
Сильное
электрическое
поле
ниже
заряженного облака создаст ионы в окрестности острых точек. Эти ионы
увеличивают проводимость воздуха и вероятность облачного разряда (удара)
[14].
Когда облако разряжается, точка удара молнии будет центром волны,
которая рассеивается в земле. Волна кажется круглыми кольцами вокруг точки
удара. Эти волны движутся наружу, очень похоже на рябь воды.
Скорость распространения волны соответствует скорости, с которой
энергия будет рассеиваться в землю. Сопротивление почвы и индуктивности
контура
заземления
образуют
постоянную
времени.
Волна
будет
распространяться быстрее с большей площадью контакта с грунтом и более
низким сопротивлением [1].
1.3 Параметры и характеристики промышленных молниеотводов
Для
нейтрализации
воздействия
атмосферных
электроразрядов
используется комплекс мероприятий и оборудований, таких как молниезащита.
Данные мероприятия обеспечивают безопасность людей, предотвращение
разрушений и пожаров зданий, различного оборудования.
10
Признаком классификации зданий и сооружений является характер
необходимых работ по молниезащите. Объекты делятся на три группы:
Первая категория – это опасные производственные объекты, где молния
может вызвать пожары, взрывы, серьезный ущерб и смерть (в местах, где
используются взрывоопасные и легковоспламеняющиеся материалы, работают
электростанции и подстанции).
Вторая категория – это взрывоопасные здания и сооружения, в которых
горючие и другие вещества хранятся в металлических или специальных
емкостях, то есть взрыв не вызовет серьезного повреждения и пожара
(хранение топлива, аммиачный холодильник, мельница).
Третья категория – объекты, которые могут быть опасны только при
прямых ударах молнии от пожара и разрушения (жилые здания, детские сады,
больницы, школы, трубы котельной и промышленные предприятия).
Здания
производств,
молниезащитой,
расположением
производится
которая
промышленные
определяется
предприятий.
посредством
Выбор
сооружения
назначением,
систем
специальных
обеспечиваются
защиты
расчетов.
строением
и
и
оборудования
Расчету
подлежит
количество ударов молнии, которые могут происходить каждые несколько
месяцев.
Для обеспечения непосредственной защиты зданий и сооружений от
ударов молний используется молниеотвод, который состоит из следующих
элементов:
1. приемника, регистрирующего разряд
2. заземлителя, отводящего ток в землю
3. токоотвода, который соединяет приемник и заземляющую систему.
В случае громоотвода разряд электричества пройдет через приемник и
обойдет защищаемый объект. Работа оборудования основана на высочайших
конструктивных особенностях и имеет хорошо заземленный удар молнии.
Молниеотводы делят на стрежневые и тросовые. Первый вариант
используется чаще, и использование кабельного оборудования ограничено
11
узкими конструкциями или объектами со многими подземными сооружениями,
которые могут помешать установке опор столбов.
Стержневые устройства могут быть:
− одиночными (антенны);
− двойными – с двумя отдельно размещенными стержнями;
− многократными – с тремя и больше стержнями, создающими общую
защитную зону.
Стержневые молниеприемники имеют длину от 200 до 1500 мм, площадь
сечения около 100 мм2.
Кабельный громоотвод также является одиночным, состоящим из кабеля
и двух поддерживающих его опор, а двойной состоит из двух одинаковых
устройств одинаковой высоты, установленных параллельно.
Объекты первой категории ниже 30 м оснащаются молниеотводами,
монтируемыми отдельно или непосредственно на здании, но изолированного от
него.
Сооружения
высотой
более
30 м
оборудуются
устройствами,
установленными не изолированно на самом здании.
Объекты второго типа защищены молниеотводами, расположенными на
здании.
Для
заземление
безопасности
металлической
третьей
группы
конструкций
кровли,
которая
используется
используется
в
качестве
молниеприемника.
Материал, используемый для изготовления молниеотвода – сталь. В
качестве устройства, принимающего удары молнии, используются различные
металлические конструкции: трубы, сетки и т.д., расположенные над
защищаемым объектом.
1.4 Электрические характеристики почв,
Электрические свойства почвы – это параметры искусственно созданных
электрических полей, которые подвержены влиянию распределения зарядов, к
примеру, неорганических ионов. Данные ионы являются итогом процесса
почвообразования. Электрические свойства почв очень важны для разделов
12
науки, таких как телекоммуникации, электротехника, геофизика и сельское
хозяйство. Существуют полуэмпирические диэлектрические модели для почв,
которые представляют действительную и мнимую части диэлектрической
проницаемости как функцию частоты. Методы измерения для получения
диэлектрических свойств грунтов описаны для разных полос частот от
нескольких кГц до нескольких ГГц. Конденсаторы с параллельными
пластинами широко используются для измерения диэлектрических свойств.
Способ линии передачи коаксиальной линии передачи может использоваться в
частотной области и временной области. Техника во временной области с
линиями передачи обычно называется рефлектометрией во временной
области, потому что она основана на измерении напряжения как функции
времени импульсов. Метод частотной области с линиями передачи основан на
измерении коэффициента отражения линии передачи. Способ линии передачи
описан с короткой нагрузкой и нагрузкой разомкнутой цепи, поскольку он
полезен для получения характеристического полного сопротивления и
электрической проницаемости среды внутри.
Почва имеет многообразный состав, в ней содержится множество
минеральных зерен. На элементах зерен образуются матрицы, минеральные и
органоминеральные, состоящие из активных центров, либо кластеров. С
помощью
катионов,
находящихся
на
данных
матрицах,
происходит
нейтрализация отрицательного заряда.
Вышесказанное определяет, что почва – носитель электрического заряда.
Коллоидная частица, находящаяся в почве, имеет двойной электрический
слой. При наложении поля на почву возникает электрокинетический
потенциал.
Электрический ток определённой силы возникает при движении ионов
(зарядов), которое осуществляется с помощью движения воды в почве. Его
улавливают электроды, помещенные в почву. Следовательно, почва имеет
стационарное электрическое поле, которое можно измерить. На потенциал
точки влияют содержание поглощенных катионов. Почвенный состав
13
увеличивает разность потенциалов в разных точках почвы, а высокая
влажность способствует выравниванию электрического поля в почве.
Объемная
плотность
характеризующим
заряда
естественное
является
основным
электрическое
параметром,
поле.
Напряжение
естественного поля измеряют в мВ, как разность потенциалов между двумя
точками, а сопротивление почвы с использованием наложения электрического
поля на почву. Обычно определяют сопротивление, которое суммирует все
сопротивления неоднородного участка почвы.
Электрическое
сопротивление
может
применяться
в
качестве
обобщенного показателя. Основываться это может только на выявленных
тесных связях его со свойствами почвы. В результате по потенциалу и
удельному электрическому сопротивлению (в зависимости от конкретных
условий) можно провести оценку большого спектра почвенных свойств.
Почкообразовательные
процессы
осуществляются
с
помощью
распределения плотностей подвижных электрических зарядов в почвах, где
каждый такой элементарный процесс будет работать на обеднение или
обогащение электрическими зарядами: горизонта, части или всей толщи
почвенного профиля [15-21].
1.5 Обзор современных методов расчета заземлителей
Система заземления должна учитывать низкий импеданс заземления, а
также низкое сопротивление.
Разработка процедур проектирования для заземления и защиты от
переходных
процессов
электрического
оборудования
основана
на
традиционных электрических концепциях, но применение может быть очень
сложным. Каждая установка уникальна. Местоположение, доступ, почвенные
условия, атмосферное взаимодействие и оборудование, подлежащее защите,
влияют на проект системы защиты. Разумные решения могут быть приняты при
наличии адекватной информации.
14
Заземлитель – это основной элемент ЗУ. Заземлитель представляет собой
одиночный заземляющий электрод или группу электродов (контур заземления),
находящихся в электрическом контакте с землей [2].
Для электрического соединения с землей электрических приборов, машин
и устройств используют ЗУ. Данное устройство состоит из металлического
электрода, находящегося под землей, а также проводника, который соединяет
заземление с заземляющей частью. Заземлители – стальные трубы, вертикально
забитые в землю, различные рельсы или стальные полосы и провода, которые
укладываются горизонтально. Стальной трос или цепь также может быть
использованы в качестве ЗУ. Электрическая система состоит из трех блоков:
рабочего заземления, защитного и заземления молниезащиты. Для обеспечения
безопасного заземления, путем уменьшения сопротивления проводника и
выравнивания
распределения
потенциала
в
области
цепи,
уменьшают
контактное и ступенчатое напряжения [5].
Конструкция и основные размеры заземлителей должны соответствовать
указанным на рисунках ниже (рис.1.1-1.3) [3].
Рис.1.1. Общий вид заземлителя
15
Рисунок 1.2 Стержни
Рис. 1.3. Зажим
16
Для предотвращения риска, связанного с поражением электрическим
током, передающимся на элементы оборудования и приборов, применяют
защитное заземление. Данный вид заземления характерен намеренному
заземлению непроводящих, металлических частей. Соответствующие части
находятся под высоким напряжением. Задача такого заземления состоит в том,
чтобы снизить напряжение относительно земли на металлических деталях,
которые находятся под напряжением. Причина напряжений может быть в
нарушении изоляции. Также данное заземление помогает предотвратить
поражение людей током при контакте с оборудованием.
Такие элементы оборудования, как корпуса насосов, конвейеров,
включателей и трансформаторов, корпуса кабельных коробок, опоры линий
электропередач,
осветительные
устройства
и
другие
комплектующие
электрооборудования могут оказывать сильное влияние на заземление.
Защитное
заземление
различных
электроустановок,
к
примеру,
стационарных или передвижных достигает порядков, более, чем 1 киловатт.
Общие заземляющие части машин и установок реализуется непрерывным
электрическим соединением между заземляющими проводами гибкого кабеля.
Через данные провода кабель соединяется с заземляющим проводником.
Заземляющая цепь в автономном режиме контролирует непрерывность ее
работы [6].
Конструкция общего ЗУ состоит из двух элементов: центрального и
локального. Первые ЗУ базируется на главной понижающей подстанцией.
Конструкция, изготовленная из стальных труб, угловой и круглой стали,
представляет из себя электроды, вертикально соединенные в цепь. Центральное
ЗУ скрыто под землей.
Такого вида заземлители в мерзлых породах представляют собой
скважину диаметров 150 – 250 миллиметров и значительно большей глубиной,
которая доходит до 50 метров. Скважину наполняют мелким грунтов, таким как
глина, песок, при этом перемешивая смесь с поваренной солью.
17
На пересечение шахтой поверхности грунта размещают отрезок трубы. К
отрезку приваривают вывод заземлителя.
Верхний конец электрода, опущенный в почву, присоединяют со
стальной полосой или круглой сталью, путем сваривания.
В качестве магистральных заземляющих проводов, прокладываемых на
опорах воздушных линий электропередачи, рекомендуется применять для
стационарных установок стальные однопроволочные и сталеалюминевые
провода.
Передвижные
электроустановки
и
пункты
их
электропитания,
выполненные на металлическом ходу (гусеницы, колеса, лыжи, опорные базы и
т.п.), имеют непосредственное соприкосновение с поверхностью грунта и
представляют собой естественные заземлители, включенные параллельно
искусственным заземлителям. Наибольшее сопротивление такая непрерывная
сеть будет иметь в наиболее удаленной от центрального заземлителя точке, т.е.
у наиболее удаленной электроустановки.
Расчет требуется для того, чтобы произвести правильный выбор
параметров искусственного заземлителя, при этом ЗУ должен соответствовать
требованиям электробезопасности [12].
При расчете ЗУ упор делают на значения токов замыкания на почву, либо
на корпус в сетях с напряжениями порядка 6 – 10 кВ. Расчет малых токов (0,4
кВ) не учитывается в связи с тем, что такие токи замыкания на землю
достаточно малы. Токи в однофазном замыкании на землю в сетях больших
напряжений, рассчитываются для выбора и настройки релейной защиты
однофазных замыканий. Также для определения величины допустимого
сопротивления защитного заземления. Величина тока однофазного замыкания
на землю в электроустановках определяется рядом факторов, из которых в
качестве основных выделяются протяженность воздушных и кабельных линий,
а также количество элементов подключенного электрооборудования [11].
Кроме того, ток однофазного замыкания определяется наличием
заземляющего провода, высотой подвеса над землей фазных и заземляющего
18
проводов воздушных линий, площадью поперечного сечения и расстоянием
между собой токоведущих и заземляющих проводников воздушных и
кабельных линий, типом подключенных элементов электрооборудования, т.е.
проводимостью относительно земли каждого из элементов электроустановки
[13].
Расчет ЗУ подразумевает собой расчет заземлителя, в виду того, что
заземляющие проводники принимаются по условиям механической прочности
и устойчивости коррозии. Не касается это лишь конструкций с выносным ЗУ.
Для них существует свой расчет. Считается последовательное сопротивление
линии, сопротивление заземлителя. Результат суммарного сопротивления не
должен превышать расчетное [7].
Расчет выполняется при наличие следующих характеристик:
- допустимое сопротивление распространения тока в земле;
- требования к конструкции искусственных заземлителей;
-
требования
к
расположению
заземлителей;
- величина токов замыкания.
вертикальных
и
горизонтальных
19
ГЛАВА 2. ТОК В ЗАЗЕМЛИТЕЛЕ
2.1 Общие положения
2.1.1 Зависимость от времени
Гармоническая
временная
зависимость
в
электродинамике
и
электротехнике от времени на языке комплексных амплитуд используется в
виде:
𝑈 𝑡 = 𝑈! cos 𝜔𝑡 + 𝜑 = 𝑅𝑒 𝑈! 𝑒 !" 𝑒 !"# = 𝑅𝑒 𝑈! 𝑒 !"# ,
где 𝑈! - комплексная амплитуда;
i – мнимая единица;
𝜔 = 2𝜋𝑓 - круговая частота;
f – циклическая частота.
Это означает, что для комплексных амплитуд операция взятия производной по
времени сводится просто к умножению на i𝜔.
Для распространяющихся волн зависимость комплексных амплитуд
полей расстояния R при этом выглядит таким образом:
exp −𝑖𝑘𝑅 ,
мнимые части комплексной диэлектрической проницаемости и волнового числа
k (характеризующие потери) должны быть меньше нуля и т.п.
Часто используют зависимость:
exp −𝑖𝜔𝑡 ;
иногда для получения мгновенного значения берут не реальную, а мнимую
часть.
Однако, при этом в электротехнике комплексное сопротивление,
например, индуктивности, выглядело бы так:
−𝑖𝜔𝐿.
Поэтому в электротехнике приходится действовать иначе. Все это создает
дополнительные трудности. Для теоретиков привычен такой формат (или
система
Гаусса
–
это
особенно
удобно
в
задачах
релятивистской
20
электродинамики). Практики в основном применяют наш формат и систему
СИ.
Мы так же здесь будем применять этот подход. Для сравнения с
результатами иных исследований надо в конечных формулах использовать
процедуру комплексного сопряжения (т.е. везде, где встречается мнимая
единица, надо менять у нее знак на противоположный)
2.1.2 Обозначения
Напряжения и ток как функции от времени будем обозначать как u(t) и
j(t); для гармонических во времени сигналов с круговой частотой 𝜔 применяем
метод комплексных амплитуд так, как отмечено ранее. Соответствующие
величины обозначаем как 𝑈! и 𝐼! (для упрощения записи – без точек над
буквами).
Постоянная распространения (волновое число) в воздухе:
𝑘! =
𝜔 2𝜋
=
,
𝑐
𝜆
где с – скорость света в вакууме;
𝜆 – длина волны в воздухе.
Если среда имеет относительную диэлектрическую проницаемость 𝜀 и
удельную проводимость:
𝜎=
1
𝜌!
(𝜌! – удельное сопротивление среды), то соответствующая комплексная
относительная диэлектрическая проницаемость равна:
𝜀=
𝜀 − 𝑖𝜎
,
𝜔𝜀!
а волновое число становится комплексным и равным:
𝑘=
!
𝜀 ! 𝑘!
;
где 𝜀! – электрическая постоянная.
Заземлитель представляет собой идеально проводящий стержень –
цилиндр кругового сечения длинной L радиуса a.
21
2.2 Задача электростатики
1) Между задачами электростатики в диэлектрике и задачами постоянных
токов в полупроводящей среде (так называемыми стационарными задачами)
много общего [4].
Итак,
сначала
электростатика.
Имеем
бесконечную
среду
с
относительной диэлектрической проницаемостью 𝜀 [10].
Для точечного заряда q потенциал:
𝜑=
!
!!!! !!
,
напряженность поля:
!
𝐸! = −𝑔𝑟𝑎𝑑𝜑 = !!!
! !!
!
.
Если заряд не точечный, а равномерно распределен по поверхности
сферы радиуса a, то на поверхности сферы:
!
𝜑! = !!!
!
!
, 𝐸! ! = !!!
!!
! !!
!
.
Поверхностная плотность заряда:
!
𝜎 = 𝜕!! = 𝜀𝜀! 𝐸!! = !!!! .
2) Возьмем более сложный случай – проводящий эллипсоид с зарядом q и
полуосями a, b и l (рис. 2.1). Имеется проводящий стержень длиной 2L с
зарядом
q,
помещенный
в
безграничную
среду
с
относительной
диэлектрической проницаемостью 𝜀. Требуется определить потенциал стержня
𝜑! . Точно такая задача не решается, однако в задачнике по электродинамике
В.В. Батыгина и И.Н. Топтыгина под номером 194 для сильно вытянутого
эллипсоида вращения с одной длинной полуосью L и двумя другими
одинаковыми полуосями, равными a. Ее решение для случая a=b=l приведено
выше [9]. А для случая цилиндра a=b, l>>a:
22
Рис. 2.1. Проводящий эллипсоид
𝜑! =
𝑞
𝑙𝑛
8𝜋𝜀! 𝜀 𝑙 ! − 𝑎!
𝜎 𝑥 =
𝑙 + 𝑙! − 𝑎!
𝑙 − 𝑙! − 𝑎!
𝑞
4𝜋𝑎! 𝑙
𝜎 0 =
𝑞
,𝜎 𝑙 =
4𝜋𝑎𝑙
𝐸!! =
≈
𝑥!
1
+
𝑎! 𝑙!
𝑞
4𝜋𝑎! 𝑙
𝑞
2𝑙
𝑙𝑛 .
4𝜋𝜀! 𝜀𝑙 𝑎
.
1
1
+ !
!
𝑎
𝑙
≈
𝑞
.
4𝜋𝑎𝑙
𝜎
𝑞
=
.
𝜀𝜀! 4𝜋𝜀𝜀! 𝑎𝑙
2.3 Задача постоянных токов
3) Рассмотрим такую же геометрию, но теперь среда проводящая с
удельной проводимостью 𝜎! (рис. 2.2), а в стержень извне по тонкому проводу
втекает постоянный ток I. Предполагаем, что второй полюс источника тока
подключен к бесконечно удаленной точке среды; таким образом, цепь
23
замкнута. С математической точки зрения эта задача аналогична предыдущей,
поэтому решение описывается такой же формулой:
𝚥 = 𝜎! 𝐸.
Если имеется сторонний ток I, то:
𝜑=
!
!!!! !
, 𝑗! =
!
!!! !
.
Если имеем проводящую сферу радиуса a в полупроводящей среде с
удельной проводимостью 𝜎! , то:
𝜑
= 𝜑! =
!!!
!
!!!! !
; 𝑗! =
!
!!! !
.
Для вытянутого эллипсоида:
𝜑! =
!
!!!! !
𝑙𝑛
!!
!
; 𝑗! =
!
!!"#
Из (формулы (2.1)) следует, что
𝐼=
4𝜋𝜎! 𝑙
𝜑 ,
2𝑙 !
𝑙𝑛
𝑎
а входное сопротивление z:
𝑧=
!!
!
=
!"
!!
!
!!!! !
.
Рис. 2.2. Полупроводящая среда
.
(2.1)
24
Таким образом, если приложить к этому идеально проводящему стержню
потенциал 𝜑! (относительно бесконечно удаленной точке в полупроводящей
среде с удельной проводимостью 𝜎! ), то потечет ток I:
4𝜋𝜎! 𝑙
𝜑 .
2𝑙 !
𝑙𝑛 𝑎
Очевидно, что если разрезать картину ровно посередине стержня
𝐼=
перпендикулярно ему, а отрезанную часть заменить воздухом, то результат
изменится ровно в два раза:
𝜑! =
𝐼
2𝑙
𝑙𝑛 .
2𝜋𝜎! 𝑙 𝑎
Чтобы перейти к напряжениям, а не потенциалам (неизвестно, от чего
отсчитанным), видоизменим задачу: над поверхностью земли параллельно ей
поместим проводящую бесконечную плоскость, обладающую потенциалом,
равным нулю. Тогда источник тока I, или, если угодно источник напряжения
𝑢 = 𝜑! , будет включен между приземленным концом стержня и проводящей
плоскостью.
В такой постановке сопротивление системы равно:
2𝑙
𝑙𝑛
𝑢
𝑎 .
𝑧= =
𝐼 2𝜋𝜎! 𝑙
Теперь рассмотрим задачу электродинамики.
Основная проблема состоит в том, как подключить источник и что такое
земля?
В статике можно говорить о потенциалах, так как эти задачи и поля
потенциальные. В электродинамике основное понятие – поле.
В электродинамике и теории антенн ставится задача так: имеется
проводящий стержень длиной 2l (рис. 2.3), радиуса a, и в середине вставляется
сторонняя ЭДС:
25
Рисунок 2.3. Проводящий стержень
Возбуждается ток I(x) (текущий по поверхности), составляющего
интегральное уравнение для тока из условия, что полное поле ( как создаваемое
ЭДС, так и возбужденным током) на поверхности стержня должно равняться
нулю (граничное условие).
2.4 Эквивалентная схема
Теперь разберемся с заземлителем.
26
Рисунок 2.4. Заземлитель
Электротехники задают уровень и форму напряжения (потенциалы)
который образуется в точке входа тока. Эквивалентная схема выглядит так:
Рисунок 2.6. Эквивалентная схема
У нас получился несимметричный вибратор.
27
𝜀
1
𝑍нв = 2 = 𝑍св,
𝐼 0
2
𝜀
𝑍св =
.
𝐼 0
Из симметричного вибратора, путем рассечения идеально проводящей
плоскостью
можно
получить
несимметричный.
Плоскость
проходит
перпендикулярно оси, через точку нулевого потенциала. Любая из половин
симметричного вибратора, учитывая проводящую плоскость, на практике
может рассматриваться как несимметричный вибратор.
2.5 Антенная задача
Итак, рассматриваем несимметричный вибратор:
𝑍=
1
𝑍 .
2 св
Из курса антенн известно (рис. 2.6) [8]:
Рисунок 2.6. Вибратор
𝑍св = 𝑅! − 𝑖𝑊𝑐𝑡𝑔𝑘𝐿 = 𝑅! + 𝑖𝑥! .
Для 𝑘𝐿 ≪ 1 → 𝑥! ≫ 𝑅! .
Поэтому:
𝑍св = −𝑖𝑊𝑐𝑡𝑔𝑘𝐿.
28
Пусть в воздухе помещен наш стержень длиной 2L радиуса а, в середине
которого подключена гармоническая ЭДС с комплексной амплитудой 𝑈! .
Требуется определить входное сопротивление Z или входной ток
𝐼=
𝑈
.
𝑍
Строгих решений здесь нет, однако антеннщики применяют здесь
приближение длинной линии, тогда:
𝑍 = −𝑖𝑊𝑐𝑡𝑔𝑘! 𝐿.
(2.2)
Такое приближение для так называемой однопроводной линии тем
лучше, чем меньше электрическая длина 𝑘! 𝐿.
Проблема, однако, состоит в том, что не очень понятно, что следует
понимать под волновым сопротивлением W в этом соотношении; линия-то
неоднородна
по
длине.
В
результате
долгих
исследований
(и
экспериментальных, и теоретических) принято использовать выражение (если
дело происходит в воздухе):
𝑊=
𝐿
𝜇! 𝑙𝑛 𝑎 − 1
𝐿
= 120 𝑙𝑛 − 1 .
𝜀!
𝜋
𝑎
Если вокруг среда с 𝜀 ≠ 1, то:
𝐿
𝜇! 𝑙𝑛 𝑎 − 1
.
𝜀! 𝜀
𝜋
𝑊=
Как видно, (в формуле (2.2)) отсутствует вещественная часть – так
называемое сопротивление излучения; при малой электрической длине она
существенно меньше мнимой части, а в рассматриваемой здесь задаче она тем
более не будет играть никакой роли.
Если
стержень
окружен
полупроводящей
средой
с
комплексной
относительной диэлектрической проницаемостью 𝜀 (а так же если в среде есть
потери с удельной проводимостью 𝜎), то логично ввести этот параметр (в
формулу (2.2)), тогда:
29
120
𝑍 = 𝑍 𝜔 = −𝑖
𝜀−𝑖
где 𝑘 =
Когда
𝜀−𝑖
!
!!!
имеется
𝐿
𝑙𝑛 − 1 𝑐𝑡𝑔𝑘𝐿 ,
𝑎
𝜎
𝜔𝜀!
𝑘! .
(2.3)
только
половина
закопанного
стержня,
то
при
подключении гармонической ЭДС точно также, как в задаче постоянных токов
входное сопротивление заземлителя определится так:
𝑍 𝜔 =
!(!)
!(!)
!"
= −𝑖
!!!
!
!
!!!
𝑙𝑛 − 1 𝑐𝑡𝑔𝑘𝐿.
(2.4)
!
Зная 𝑍 𝜔 и 𝑈(𝜔), можно найти ток I(𝜔), втекающий в заземлитель:
𝐼 𝜔 =
𝑈(𝜔)
.
𝑍(𝜔)
(2.5)
Для проверки рассмотрим ситуацию, когда 𝜔 стремится к нулю (и
соответственно модуль kL и tgkL стремятся к 0):
𝑍 𝜔 =
!!!"
!
!!!
!!!
!
!!!"
𝑙𝑛 ! − 1 𝑐𝑡𝑔𝑘𝐿 ≈
!
!
!!!
!!!
!
𝑙𝑛 ! − 1
!" !!!
!
!" !!
!
!!!
.
Или
𝑍 𝜔 =𝑍=
!
!
=
!"!!!
!
𝑙𝑛 ! − 1 =
!
!
!!"#
.
(2.6)
Сравнение с (формулой (2.3)) показывает, что здесь результат несколько
меньший. Это, по-видимому, связано с тем, что в задаче постоянных токов
рассматривается эллипсоид с переменным радиусом, изменяющимся от а до 0, а
не цилиндр постоянного радиуса а.
2.6 Импульсный ток
2.6.1 Импульсное напряжение
Предположим, что напряжение u(t) представляет собой импульс с малым
фронтом по сравнению с длительностью спада.
Его формулу наиболее просто описать в виде разности двух экспонент:
30
𝑓 𝑡 = exp −𝛼𝑡 − exp −𝛽𝑡 ,
причем 𝛽 ≫ 𝛼.
Так как при малых значениях аргумента 𝑥 ≪ 1 экспонента описывается
следующим образом:
𝑒
!!
𝑥! 𝑥!
= 1 − 𝑥 + − + ⋯ ≈ 1 − 𝑥,
2
6
то при 𝛼𝑡 < 𝛽𝑡 ≪ 1 функция будет выглядеть следующим образом:
𝑓 𝑡 = 1 − αt − 1 + 𝛽𝑡 = 𝛽 − 𝛼 𝑡.
Это обеспечивает крутой фронт импульса в его начале. Но довольно
быстро 𝛽𝑡 становится большим единицы и 𝑒 !!" перестает играть какую-либо
роль. Поэтому главным становится 𝑒 !!" - она и определяет медленный спад
импульса ( в сравнении с быстрым ростом 𝛽 − 𝛼 𝑡).
Максимум импульса в момент 𝑡! определяется экстремумом функции
𝑓 𝑡 .
Для построения 𝑡! надо положить
𝑓 ! 𝑡 = 0, т.е. −α𝑒 !!" + 𝛽𝑒 !!" = 0.
Максимальное значение импульса наступает при
𝛽
𝛼
𝑡! = ln
𝛽−𝛼
и равно
𝛽
𝛽
𝛼𝑙𝑛
𝛼 − 𝑒𝑥𝑝 −
𝛼 .
𝑏 = exp −
𝛽−𝛼
𝛽−𝛼
𝛼𝑙𝑛
Время 𝑡! следует считать длительностью фронта импульса.
Понятно, что импульс тянется теоретически до бесконечности. Поэтому
задается уровнем падения его от максимума в какой-то момент времени,
например в 2 раза.
Т.е.:
𝑓 𝑡! =
1
1
𝑓 𝑡! = 𝑏.
2
2
31
Или:
1
𝑏.
2
𝑒 !!!! − 𝑒 !!!! =
В принципе, задавшись значением 𝑡! и 𝑡! , можно найти параметры
𝛼
импульса
𝛽.
и
Однако,
решение
не
простое,
т.к.
тригонометрическое. Поэтому будем использовать приближения.
Итак, возьмем задачи 𝑡! = 10!! сек, 𝑡! = 2 ∙ 10!! сек = 50𝑡! .
Предположим, что 𝛽 ≫ 𝛼.
Тогда:
𝛽
𝛽
𝛼
𝛼
𝑡! = ln
≈ ln
.
𝛽−𝛼
𝛽
Величина ln
!
!
также будет большой, поэтому при 𝑡 = 𝑡! :
𝑒 !!!! =
!
!"#
!
!
𝑒 !
≈𝑒
! !
! !"
! !
!
≈ 𝑒 !!"! ≪ 1.
Поэтому:
𝑏≅
!
!"#
!
!
𝑒 !!!
!
≈𝑒
!
!"
!
!
! .
!
Т.к. ! ≫ 1, то:
−
!
!
!
!
!"
!
!
≪ 1 и 𝑏 ≅ 1 − ! 𝑙𝑛 ! ≈ 1.
Итак, если:
𝑓 𝑡! =
1
1
𝑏≈ ,
2
2
то:
𝑒 !!!! − 𝑒 !!!! ≈ 𝑒 !!!! =
1
,
2
𝛼𝑡! = 𝑙𝑛2 ≈ 0,7.
Тогда:
𝛼𝑡! = 𝛼𝑡!
𝑡!
= 𝛼𝑡! ∙ 50 = 0,7.
𝑡!
уравнение
32
0,7
≈ 0,014;
50
0,014 0,014
𝛼=
=
= 1400.
𝑡!
10!!
𝛼𝑡! =
Итак:
𝛽
𝛼 𝛽 𝑙𝑛 𝛼
𝛼𝑡! = 𝑙𝑛 =
= 0,014.
𝛽
𝛽 𝛼
𝛼
Последняя запись – это граничные уравнения относительно величины
!
!
.
Его решение:
𝛽
≈ 400.
𝛼
Тогда:
𝛽 = 400𝛼 = 560000.
Итак:
при 𝑡! = 10 мксек = 10!! сек,
𝑡! = 500 мксек = 50𝑡! , при 𝑓 𝑡! =
1
𝑓 𝑡! ,
2
где:
𝑓 = 𝑓 𝑡! ≅ 1, 𝑓 𝑡! =
1
,
2
получаем:
𝛼 = 1400 сек!! , 𝛽 = 560000 сек!! ,
𝑓 𝑡 = 𝑒 !!" − 𝑒 !!" .
𝑓 𝑡 - это нормированный импульс напряжения, где 𝑓 = 𝑏 ≈ 1. А вообще
задается максимальное значение напряжения 𝑈!"# .
Задавшись величиной максимального значения напряжения 𝑈!"# ,
получим выражение для импульсного напряжения, создающего ток в
заземлителе:
𝑢 𝑡 =
𝑈!"#
exp −𝛼𝑡 − exp −𝛽𝑡 .
𝑏
33
2.6.2 Спектры сигналов. Восстановление функции i(t)
Для получения 𝐼(𝜔) (по формуле (2.4)) найдем спектральную функцию
𝑈(𝜔) с помощью преобразования Фурье:
𝑈 𝜔 =
!
𝑢
!
𝑡 exp −𝑖𝜔𝑡 𝑑𝑡 =
!!"#
!!!
!
!"!! ! !!"(!!!)
.
Временная функция тока:
1
𝑗 𝑡 =
2𝜋
!!
𝐼(𝜔)𝑒
!"#
!!
1
𝑑𝜔 = 𝑅𝑒
𝜋
!
𝐼(𝜔) 𝑒
!"#
!
1
𝑑𝜔 = 𝑅𝑒
𝜋
!
!
𝑈(𝜔)𝑒 !"# 𝑑𝜔
.
𝑍(𝜔)
Тогда, применив обратное преобразование Фурье (к формуле (2.5)), получим:
𝑗 𝑡 =
(𝛽 − 𝛼)𝑈!"#
𝑅𝑒
𝐿
60𝜋𝑏(𝑙𝑛 𝑎 − 1)
!
!
𝜎
𝜀 − 𝑖 𝜔𝜀 𝑡𝑔
!
𝜎 𝜔
𝜀 − 𝜔𝜀 𝑐 𝐿 𝑑𝜔
!
𝛼𝛽 − 𝜔 ! + 𝑖𝜔 𝛼 + 𝛽
.
Т.к.
𝑒 !" − 𝑒 !!"
𝑠𝑖𝑛𝑥
𝑒 !" − 𝑒 !!"
1 − 𝑒 !!!!
2𝑖
𝑡𝑔𝑥 =
=
= −𝑖 !"
= −𝑖
,
𝑐𝑜𝑠𝑥 𝑒 !" + 𝑒 !!"
𝑒 + 𝑒 !!"
1 + 𝑒 !!!!
2
то получаем следующую формулу:
ο
[1 − exp( −i 2kL )] exp(iωt )dω
∞
ωε0
( β − α )U max
j (t ) =
∫−∞ [αβ − ω 2 + iω(α + β )][1 + exp( −i 2kL)] ] .
L
120πb(ln − 1)
ε −i
a
Или:
ο
[1 − exp( −i 2kL )] exp(iωt )dω
ωε 0
( β − α )U max
j (t ) =
Re ∫
.
2
L
[
αβ
−
ω
+
i
ω
(
α
+
β
)][
1
+
exp(
−
i
2
kL
)]
0
60πb(ln − 1)
∞
ε −i
a
Этот интеграл считается приближенно численно.
2.7 Колебательные процессы
Когда возможны колебательные процессы в рассматриваемой системе?
Вся
физика
определяется
выражении для тока:
особенностями
подинтегральной
функции
в
34
1
j (t ) =
2π
∞
U (ω )dω
Z (ω )
−∞
∫
Предположим, что спектр напряжения U(w) содержит достаточно
высокие частоты, позволяющие возбудить эти колебания. Для этого фронты
напряжения u(t) должны быть достаточно малыми. Итак, особенности в – это
частоты, при которых Z(w) обращается в 0. К сожалению, здесь нет чистых
полюсов подинтегральной функции; тогда можно было бы получить
относительно простые выражения для тока. Приходится иметь дело с
квадратными корнями, где есть точки ветвления и т.д.
Будем искать значения w, пусть и комплексные, из уравнения
Z (ω ) = −i
60
ε −i
σ
ωε0
(ln
L
− 1)ctgkL = 0 .
a
Или
− ictgkL =
1 + exp( −i 2kL )
= 0.
1 − exp( −i 2kL )
Один из корней (для наименьших по модулю значений комплексной
частоты) – это: 2kL=p. Подставляя сюда выражения для k0, e и оставляя в
стороне несложные математические преобразования, найдем значение для
корня:
ω = ω02 − α 2 + iα ,
где ω0 = 2πf 0 =
πc
,α= σ .
2εε 0
2 εL
Форма колебанйи во времени:
i(t ) = exp( −αt ) sin( ω02 − α 2 t ) ,
очень похожая на свободные колебания в последовательном контуре.
Очевидно, что при w0<a режим будет апериодическим. Возвращаясь к нашей
35
задаче, заметим, что в колебательном режиме (при σ <<
ε
120L
) для L=15 м и e=3
частота колебаний f0 составляет примерно 3 МГц, а период T=1/f0=0.3 мксек;
при этом s<<0.001 См/м, а удельное сопротивление грунта должно быть
существенно больше 1000 Ом*м.
Очевидно, что фронт импульса напряжения должен составлять доли мкс.
36
ГЛАВА 3. ПРИМЕРЫ РАСЧЕТОВ
3.1 Программа
Для построения графиков в среде fortran была разработана программа,
приведенная ниже:
j=cmplx(0.,1.)
t1=dalog(beta/alfa)
t1=t1/(beta-alfa)
b=dexp(-alfa*t1)-dexp(-beta*t1)
delt=tmaximp/nt
do 10 mt=1,nt
print *,m
t=delt*mt
f=dexp(-alfa*t)-dexp(-beta*t)
f=f/b
tt=t*1000.
write(2,1)tt,f
10
1
continue
format(2(1x,e11.4))
aa=60.*pi*b*(dalog(L/a)-1.)
bb=(beta-alfa)*umax/aa
ommax=10.*dsqrt(alfa*beta)
delt=tmax/nt
do 20 mt=1,nt
t=mt*delt
tt=t*1000.
delom=.1/t
n=int(ommax/delom)
e1=cmplx(0.,0.)
37
do 30 m=1,n
omega=delom*(m-.5)
eps=eps1-j*sigma*1.e12/(omega*8.85)
k=omega*eps/3.e8
e2=cdexp(-j*2.*k*L)
e3=cdexp(j*omega*t)
e4=alfa*beta-omega**2+j*omega*(alfa+beta)
e5=e4*(1.+e2)
e5=eps*(1.-e2)*e3/e5
e1=e1+e5
30
continue
e1=e1*bb*delom
cc=real(e1)
write(3,1)tt,cc
20
continue
stop
end
3.2 Построение графиков
3.2.1 График импульса напряжения
Для параметров земли, геометрии и т.п., характеристики которых
приведены ниже, построен график импульса напряжения (рис. 3.1):
500 - nt - число отсчетов по времени;
15000. - alfa (1/сек) - первый параметр импульса напряжения;
6.e6 - beta (1/сек) - второй параметр импульса напряжения;
2.e-4 - tmax (сек) - максимальное время просмотра тока;
1.e-3 - tmaximp (сек) - конец отсчета времени для просмотра импульса
напряжения;
2.e5 - umax (В) - максимальное значение напряжения;
3. - eps1 - относительная диэлектрическая проницаемость земли;
38
1.e-3 - sigma (См/м) - удельная проводимость земли;
.05 - a (м) - радиус заземлителя;
15.- L (м) - длина заземлителя.
Рисунок 3.1. Импульс напряжения
3.2.2 Графики тока
Для импульса напряжения (рис. 3.1) и параметров земли, геометрии и т.п.,
характеристики которых приведены ниже, построен график тока (рис. 3.2):
500 - nt - число отсчетов по времени;
15000. - alfa (1/сек) - первый параметр импульса напряжения;
6.e6 - beta (1/сек) - второй параметр импульса напряжения;
2.e-4 - tmax (сек) - максимальное время просмотра тока;
39
1.e-3 - tmaximp (сек) - конец отсчета времени для просмотра импульса
напряжения;
2.e5 - umax (В) - максимальное значение напряжения;
3. - eps1 - относительная диэлектрическая проницаемость земли;
1.e-3 - sigma (См/м) - удельная проводимость земли;
.05 - a (м) - радиус заземлителя;
15.- L (м) - длина заземлителя.
Рисунок 3.2. График тока
Параметры системы остались следующие L=15м, a=5см, e грунта равно 3.
Представлен график тока (рис. 3.3) для s=0.001См/м и длительности фронта
импульса напряжения, равной 1 мксек. При таких параметрах альфа и бета –
импульс тока апериодический. Параметры построения графика следующие:
40
500 - nt - число отсчетов по времени;
15000. - alfa (1/сек) - первый параметр импульса напряжения;
6.e6 - beta (1/сек) - второй параметр импульса напряжения;
5.e-4 - tmax (сек) - максимальное время просмотра тока;
1.e-5 - tmaximp (сек) - конец отсчета времени для просмотра импульса
напряжения;
2.e5 - umax (В) - максимальное значение напряжения;
3. - eps1 - относительная диэлектрическая проницаемость земли;
1.e-3 - sigma (См/м) - удельная проводимость земли;
.05 - a (м) - радиус заземлителя;
15.- L (м) - длина заземлителя.
Рисунок 3.3. График тока
41
Отличие данных для графика тока ниже (рис. 3.4) (приведены ниже) по
сравнению с предыдущим вариантом заключается в том, что удельная
проводимость грунта снижена. Здесь заметны затухающие колебания с
квазипериодом, примерно равным 0.3 мксек, что соответствует теории.
500 - nt - число отсчетов по времени;
150000. - alfa (1/сек) - первый параметр импульса напряжения;
6.e7 - beta (1/сек) - второй параметр импульса напряжения;
3.e-6 - tmax (сек) - максимальное время просмотра тока;
2.e-6 - tmaximp (сек) - конец отсчета времени для просмотра импульса
напряжения;
2.e5 - umax (В) - максимальное значение напряжения;
3. - eps1 - относительная диэлектрическая проницаемость земли;
1.e-4 - sigma (См/м) - удельная проводимость земли;
.05 - a (м) - радиус заземлителя;
15.- L (м) - длина заземлителя.
42
Рисунок 3.4. График тока
Для графика тока (рис. 3.5) параметры следующие:
500 - nt - число отсчетов по времени;
150000. - alfa (1/сек) - первый параметр импульса напряжения;
6.e7 - beta (1/сек) - второй параметр импульса напряжения;
3.e-6 - tmax (сек) - максимальное время просмотра тока;
1.e-3 - tmaximp (сек) - конец отсчета времени для просмотра импульса
напряжения;
2.e5 - umax (В) - максимальное значение напряжения;
3. - eps1 - относительная диэлектрическая проницаемость земли;
5.e-5 - sigma (См/м) - удельная проводимость земли;
.05 - a (м) - радиус заземлителя;
15.- L (м) - длина заземлителя.
43
Рисунок 3.5. График тока
Наконец, ниже (рис. 3.6) представлен график при тех же условиях, однако
при длительности фронта импульса напряжения, увеличенном в 10 раз – до 15
мксек:
500 - nt - число отсчетов по времени;
1500000. - alfa (1/сек) - первый параметр импульса напряжения;
6.e8 - beta (1/сек) - второй параметр импульса напряжения;
3.e-6 - tmax (сек) - максимальное время просмотра тока;
1.e-3 - tmaximp (сек) - конец отсчета времени для просмотра импульса
напряжения;
44
2.e5 - umax (В) - максимальное значение напряжения;
3. - eps1 - относительная диэлектрическая проницаемость земли;
5.e-5 - sigma (См/м) - удельная проводимость земли;
.05 - a (м) - радиус заземлителя;
15.- L (м) - длина заземлителя.
При этом на резонансной частоте системы уровень задающего сигнала
существенно уменьшился. Это привело к уменьшению размаха тока в
несколько раз.
Рисунок 3.6. График тока
Если фронт импульса напряжения еще уменьшить, то форма колебаний
может отличаться от синусоидальной.
45
3.3 О характере Z(w)
При w=0 сопротивление чисто активное и определяется соотношением (в
формуле (2.6)). Если частоты невелики (или длина L небольшая), когда модуль
kL существенно меньше единицы, величина Z имеет емкостной характер и
представляет
собой
параллельное
соединение
емкости
и
активного
сопротивления:
2πLεε 0
1
2πLσ
= Y (ω ) = g + ib =
+ iω
.
L
L
Z (ω )
ln − 1
ln − 1
a
a
Напротив, если модуль мнимой части kL значительно больше единицы,
т.е. сtgkL стремится к единице, величина Z определяется выражением:
120(ln
Z (ω ) =
ε −i
L
− 1)
a
.
σ
ωε0
Действительно, мнимая часть этого выражения больше 0, однако
зависимость ее от частоты (как, впрочем, и вещественная часть) имеет
странный характер.
46
ЗАКЛЮЧЕНИЕ
В данной работе объектом исследования были процессы тока в
заземлителе. Мы применили радиофизический подход. Мы использовали
аппарат антенн и исследовали при каких параметрах почвы, размерах возникает
колебательный процесс.
Все это полезно для тех, кто разрабатывает заземлители и занимается
созданием молниеотводов. Это дает возможность воспользоваться результатом
формул
для
исследования
тока
для
любых
импульсов,
параметров
заземлителей, земли и т.д.
В ходе работы по поставленной задаче, были рассмотрены современных
методы расчета заземлителей.
Была получена модель цилиндрического заземлителя, погруженного в
полупроводящую среду.
Исследован процесс распространения волн тока в заземлителе как в
несимметричном вибраторе. Разработана электродинамическая модель и
алгоритм расчета.
Рассмотрены переходные процессы при импульсном воздействии.
Рассчитаны уровни и формы токов.
Были получены аналитически формулы тока, на основе представления
заземлителя как электрического вибратора, погруженного в среду.
Показано, что колебательные процессы можно представить как аналог
колебательного контура.
Данные, полученные в этой работе можно использовать в следующих
случаях:
1. Когда система очень сложная – много горизонтальных проводящих
проводов стоят в разных местах и т.д. Может быть не 1 заземлитель, а 2 и
более.
2. Когда спадающий импульс тока, то все понятно. Электротехники часто
для ограничения тока ставят ограничительные диоды, которые не позволяют
47
каким-то токам переключаться. У них сложные системы, силовые. Когда
возникают немонотонные системы, не просто снижается, а ходит «туда-сюда»,
они хотят узнать при каких условиях это происходит и что в данном случае
нужно делать.
3. С помощью формул можно определить какие проводимости когда
возникают, на какой длине и каких частотах.
4. Так же это вызывает большой интерес. Как можно было ожидать, что
молния будет такой непериодической. Связано это с заземлителем, с его
длинной в первую очередь. Мы эквивалентную схему представляем в виде
резонансного контура.
48
СПИСОК ЛИТЕРАТУРЫ
1.
Боронин В.Н., Коровкин Н.В., Кривошеев С.И., Шишигин С.Л.,
Миневич Т.Г. Нетреба К.И. Математическое моделирование заземляющих
устройств при действии импульсных токов // Изв. Рос. акад. наук. Энергетика.
2013. № 6. С. 8089.
2.
Бургсдорф
В.В.,
Якобе
А.И.
Заземляющие
устройства
электроустановок. М.: Энергоатомиздат, 1987. – 400 с.
3.
ГОСТ 16556-2016 Заземлители для передвижных электроустановок.
Общие технические условия.
4.
Демирчян, К.С. Теоретические основы электротехники. Учебник
для вузов. 5-е изд. Том 2 / К.С. Демирчян, Л.Р. Нейман, Н.В. Коровкин. - СанктПетербург: Питер, 2009. – 432 с.
5.
Заземляющие
устройства
электроустановок
(требования
нормативных документов, расчет, проектирование, конструкции, сооружение):
справочник / Р.К. Борисов и др. — М.: Издательский дом МЭИ, 2013. - 360 а
6.
Карякин
Р.Н.
Заземляющие
устройства
электроустановок:
справочник / Р.Н. Карякин. – М.: Энергосервис, 2002. – 218 с.
7.
Коровкин, Н.В. Расчетные методы в теории заземления. Научно-
технические ведомости СПбГПУ / Н.В. Коровкин, С.Л. Шишигин // Изд-во
СПбГПУ. Вып. 1(166). - 2013. - С.74-79.
8.
Сазонов Д.М. Антенны и устройства СВЧ – Москва – 1988 – 434 с.
9.
Топтыгин И.Н. Современная электродинамика – Москва-Ижевск –
2002 – 763 с.
10.
Тихонов А.Н., Самарский А.А. Уравнение математической физики
– 1963 – 742 с.
11.
Черепанов,
А.В.
Расчет
заземлителей
с
учетом
частотной
зависимости удельного сопротивления грунта / А.В. Черепанов, С.Л. Шишигин
// Труды 4 всерос. науч. техн. конф. «Техно-ЭМС 2017»., М. – 2017. – С.110 –
112.
49
12.
Шишигин, С.Л. Математические модели и методы расчета
заземляющих устройств / С.Л. Шишигин // Электричество. - 2010. - №1. - C.1623.
13.
Электромагнитная
совместимость
в
электроэнергетике
и
электротехнике / А.Ф. Дьяков, И.П. Кужекин, Б.К. Максимов, А.Г. Темников;
под ред. А.Ф. Дьякова. - Москва: Изд. дом МЭИ, 2009. - 455 с.
14.
Якушева, Ю. С. Грозозащита и заземление. Молниеотводы,
грозозащитные тросы и разрядники // Актуальные проблемы энергетики :
материалы 70-й научно-технической конференции студентов и аспирантов /
Белорусский
национальный
технический
университет,
Энергетический
факультет. Секция 2: Электрические системы. – Минск : БНТУ. — 2014. — С.
104 – 105.
15.
Akbari M., Sheshyekani K., Alemi M. R. The effect of frequency
dependence of soil electrical parameters on the lightning performance of grounding
systems. // IEEE Trans. Electromagnetic Compatibility. — 2013. — vol. 55, N 4. —
P. 739-746.
16.
Cavka D., Mora N., Rachidi F. A comparison of frequency-dependent
soil models: application to the analysis of grounding systems. // IEEE Trans.
Electromagnetic Compatibility. — 2014. — vol. 56, N 2. — P. 177-187.
17.
Kherif O., Chiheb S., Teguar M., Mekhaldi A. Time-Domain modeling
of grounding systems impulse response incorporating nonlinear and frequencydependent aspects. // IEEE Trans. Electromagnetic Compatibility. — 2018. — vol.60,
N 4 — P. 907-916.
18.
Sekioka S. Frequency and current-dependent grounding resistance model
for lightning surge analysis. // IEEE Trans. Electromagnetic Compatibility. — 2018.
— vol.60, N 4. — P.907-916.
19.
Currents:
Visacro S., Rosado G. Response of Grounding Electrodes to Impulsive
An
Experimental
Evaluation.
//
IEEE
Trans.
Electromagnetic
Compatibility. — 2009. — vol. 51, N 1. — P.161-164.
20.
Visacro, S. HEM: A Model for Simulation of Lightning-Related
50
Engineering Problems / S. Visacro, A. Soares // IEEE Transactions On Power
Delivery. - Vol. 20, No. 2. - April 2005. - P.1206-1208.
21.
Zhou W., Liu Y., Li G., Yan G., Yang S., Li H., Su C. Q. Power
Frequency Interference and Suppression in Measurement of Power Transmission
Tower Grounding Resistance. IEEE Trans. Power Delivery, 2015, Vol.30, No.2,
P.1016-1023.
Отзывы:
Авторизуйтесь, чтобы оставить отзывСпасибо, достаточно интересная работа