Задача растекания тока в проводящей среде

Задача. К сферическому заземлителю (электроду) радиуса м подводится ток А (рис. 6.3.1а). Заземлитель находится на расстоянии h = 6 м от границы раздела двух сред с разными удельными проводимостями (См/м) и (См/м). Исследовать картину поля.

а) б) в)

Рисунок 6.3.1. а) взаимное расположение шарового источника тока и двух проводящих сред; б) физическая модель анализа поля тока для 1-й среды;

в) то же для 2-й среды.

Задача решается методом отражений с учетом метода электростатической аналогии. Поле растекания токов обладает осевой симметрией (нет зависимости от угловой координаты), поэтому вместо декартовой системы координат принята цилиндрическая система - z, r.

Для расчета поля в 1-й среде вводится, кроме тока , фиктивный ток согласно (3.7):

.

При этом вторая среда замещается первой (см. рис. 3.1б). Для расчета поля во 2-й среде вводится фиктивный ток

,

а первая среда замещается второй (см. рис. 3.1в).

Поле сферического заземлителя подобно полю шарового заряда, что позволяет воспользоваться выражением (6.2.1), в котором заряд заменяется на ток , а диэлектрические свойства среды на проводимость :

.

Построим картину поля в пакете MathCAD. В поле рабочего файла указываем исходные данные задачи:

Вычисляем входящие в решение постоянные коэффициенты:

Определяем длины радиус-векторов до точек наблюдения, для чего предварительно следует выбрать систему координат (произвольно). В задаче начало отсчёта переменных r, z смещено влево на расстояние 3h:

Определяем функции потенциала в первой и второй средах используя метод наложения:

Определяем потенциал поверхности заземлителя (электрода):

Определяем потенциальную функцию для всего исследуемого пространства с учетом её поведения в первой и второй средах:

Так как на поверхности и внутри электрода потенциал постоянен и определяется ранее найденной величиной , окончательно доопределим потенциальную функцию как

Построим график изменения потенциальной функции вдоль оси z при r = 0, для чего предварительно зададим диапазон изменения переменной z:

График изменения потенциальной функции показан на рис. 6.3.2. Исследуемая функция непрерывна, включая границу раздела сред. Равенство потенциалов на границе эквивалентно выполнению граничного условия (2.6).

Построим график поверхности потенциальной функции и карту линий уровня, которая задана нами ранее как функция двух переменных, для чего указываем границы расчетной области:

Рисунок 6.3.2 - График изменения потенциальной функции вдоль оси z при r=0

Определяем, сколько точек следует отложить по координатным осям:

Введением дискретных аргументов i и j индексируем точки, где определяются значения функции:

Через операцию присваивания определяем значения двумерного массива , определяя его найденной потенциальной функцией ,

и строим график поверхности потенциальной функции (рис. 6.3.3).

При использовании массива для построения карты линий равного уровня следует скопировать график поверхности (рис. 6.3.3) или сразу на нем однократно щелкнуть правой кнопкой мыши, выбрать из всплывшего контекстного меню команду “Format…”. Далее в появившемся окне 3-D Plot Format подраздел General (рис. 1.16) следует сменить тип графика с Surface Plot на Contour Plot.

Для построения карты линий равного уровня можно использовать как массив (рис. 2.4а), так и воспользоваться встроенной функцией CreateMesh (рис. 2.4б):

Рисунок 6.3.3 - Пример построения поверхности потенциальной функции (Surface Plot).

Полученные графики (рис. 6.3.4) полностью идентичны. В большинстве случаев использование функции CreateMesh более удобно, так как в этом случае (по умолчанию) размерность осей в абсолютных единицах – метрах (рис. 6.3.4б), а на рис. 6.3.4а – определяется ранее заданным числом расчетных точек по координатным осям (по оси z определено 80 точек, по оси r - 40).

а) б)

Рисунок 6.3.4 - Карты линий равного уровня потенциальной функции, построенные с применением: а – массива , б - функции CreateMesh.

Если значения функции на линиях уровня, при типе графика Contour Plot, не выведены, то следует вызвать окно 3-D Plot Format, выбрать подраздел Special, столбец Contour Options, активировать пункт Numbered – щелкнув в квадратике рядом с ним (появится галочка) и далее “Применить” (рис. 6.2.17). Нажатие кнопки ОК завершает операцию.

Для построения в MathСAD линии уровня заданного значения следует определить заданную величину как постоянную и воспользоваться встроенной функцией знака sign. В качестве примера построим для определенной ранее потенциальной функции эквипотенциаль со значением 2 кВ:

При малом числе расчетных точек полученная эквипотенциальная кривая может иметь ступенчатый вид (рис. 6.3.5а), в этом случае следует увеличить число расчетных точек. При увеличении числа расчетных точек в 5 раз по сравнению с предыдущим случаем:

получим более приемлемый вид эквипотенциальной кривой (рис. 6.3.5б).

а) б)

Рисунок 6.3.5 - Примеры построения заданной линии равного уровня с малым (а) и с увеличенным (б) числом расчетных точек

Используя возможности пакета MathСAD, найдем с помощью оператора дифференцирования (панель Calculus (рис. 6.2.11), вызывается через меню View/Toolbars/Calculus) производные от потенциальной функции, которые (с учётом знака) равны проекциям вектора напряженности поля:

На рис. 2.6 и 2.7 приведены примеры построения графиков изменения компонент напряженности вдоль различных осей, из которых видно, что нормальная составляющая вектора на границе раздела сред изменяется скачком в соответствии с граничным условием (6.3.6) и изменяет знак при переходе с левого края электрода на правый.

Рисунок 6.3.6 - График изменения z-компоненты напряженности вдоль оси z при r = 0.

Рисунок 6.2.7 - График изменения r-компоненты напряженности вдоль оси r по границе раздела проводящих сред (z = 3h).

Составляющие вектора плотности тока могут быть найдены по формуле (6.3.2)

с учетом удельной проводимости каждой из сред.

В отличие от вектора (рис. 6.3.6) вектор плотности тока на границе раздела обладает свойством непрерывности (рис. 6.3.8), так как оба этих вектора для переменной r = 0 определены только своими нормальными составляющими (см. формулы 6.3.5 и 6.3.6).

Рисунок 6.3.8 - График изменения z-компоненты плотности тока

вдоль оси z при r = 0.