Решение полевых задач электротехники с помощью вторичных источников поля и рядов Фурье

Размещено на http://www.allbest.ru/

РЕШЕНИЕ ПОЛЕВЫХ ЗАДАЧ ЭЛЕКТРОТЕХНИКИ С ПОМОЩЬЮ ВТОРИЧНЫХ ИСТОЧНИКОВ ПОЛЯ И РЯДОВ ФУРЬЕ

Попов Борис Клавдиевич, к.т.н., доцент

Попова Ольга Борисовна, к.т.н., доцент

Кубанский государственный технологический университет, Краснодар, Россия, 350072, ул. Московская, 2, pbk47@mail.ru

Для полевой задачи приведено аналитическое решение контрольного примера, с помощью которого выведены условия правильного применения метода вторичных источников поля. Найдено численное решение в виде ряда Фурье той же самой задачи, но с использованием метода вторичных источников. Предложен метод учёта неоднородности среды, который так вводит вторичные источники, что аналитическое и численное решения совпадают

Ключевые слова: АНАЛИТИЧЕСКОЕ РЕШЕНИЕ, КОНТРОЛЬНЫЙ ПРИМЕР, РЯД ФУРЬЕ, МЕТОД ВТОРИЧНЫХ ИСТОЧНИКОВ ПОЛЯ, ЧИСЛЕННОЕ РЕШЕНИЕ

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

В настоящее время для расчёта электромагнитного поля в неоднородных средах применяется метод конечных элементов, являющийся развитием теории сплайнов и конечно-разностных методов. На наш взгляд метод конечных элементов обладает рядом недостатков. Один из существенных недостатков заключается в том, что надо запоминать характеристики поля во всех точках пространства. А это громадный объём информации, который необходимо обрабатывать и хранить. Как следствие, возникает потеря точности расчёта. Этому методу есть альтернатива ? метод вторичных источников поля, который получил широкое развитие в нашей стране в конце прошлого века. Он лишён указанного недостатка, так как решается, в основном, задача определения на поверхности раздела сред плотности зарядов, либо токов, имитирующих заданную неоднородность. Для реальных технических устройств поверхности раздела сред имеют ограниченные размеры. Поэтому объём информации о поверхностных источниках меньше чем о пространственном распределении поля. Однако и у этого метода есть недостаток, показанный ниже.

Рассмотрим правомерность применения метода вторичных источников поля в предложенном виде [1 - 4] на базе контрольного примера. Затем покажем возможность применения рядов Фурье для решения полевой задачи, используя понятие обобщённой функции для вторичных источников в виде поверхностных зарядов.

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

Рисунок 1. Эскиз рассматриваемой задачи

задача электротехника фурье

Имеются две бесконечные ферромагнитные области, обращённые друг к другу поверхностями 1 и 4. Магнитная проницаемость этих областей . Между этими областями расположена ферромагнитная бесконечная вставка с магнитной проницаемостью . Магнитный потенциал на поверхности 1 - , а на поверхности 4 - .

Определим распределение магнитного потенциала и напряжённости магнитного поля в областях , и (см. рисунок 1).

Между поверхностями 1 и 4 ток отсутствует, поэтому распределение магнитного потенциала будет описываться уравнением Лапласа для одномерного случая

(1)

Имеем три одномерные взаимосвязанные граничные задачи, подчинённые следующим краевым условиям:

поверхность 1 - , при ; (2)

поверхность 2 - , при ; (3)

поверхность 3 - , при ; (4)

поверхность 4 - , при , (5)

где , и - напряжённости магнитного поля в области (поверхность 2), области (поверхности 2 или 3) и области (поверхность 3), соответственно.

Для упрощения будем считать .

Решение уравнения (1) в общем виде будет иметь вид прямой линии

(6)

где и - постоянные интегрирования, определяемые из краевых условий (2 - 5).

Учитывая указанные краевые условия и выражение

(7)

определим постоянные интегрирования в выражении (6) для всех областей.

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

область - , (8)

(9)

область - , (10)

(11)

область - , (12)

(13)

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

Напомним метод вторичных источников поля [1]. Сначала определяют электромагнитное поле в однородном пространстве от системы независимых источников поля. Затем на границах раздела сред с разными характеристиками вносят систему распределённых поверхностных зарядов либо токов (вторичные источники поля), зависящую от поля исходных источников и характеристик сред. Определяют распределение поля от этой системы вторичных источников поля. Складывают поля от исходных источников поля и вторичных источников поля и получают искомую картину электромагнитного поля. Величины вторичных источников поля на поверхности раздела сред определяют, исходя из граничных условий и теоремы Гаусса для зарядов или закона полного тока для токов.

Найдём, используя решение (6) уравнения (1) и граничные условия (2) и (5), распределение магнитного потенциала в рассмотренной выше конструкции при отсутствии магнитной области .

(14)

Учитывая выражение (7), вычислим напряжённость магнитного поля

(15)

Как было сказано выше, определим конфигурацию поля, которое нужно сложить с полем по формулам (14) и (15), чтобы получить решения (8 - 13) для неоднородной среды.

Для вычисления напряжённости , имитирующей вторичные источники поля, вычтем из решения (8 - 13) решение (14):

область - ; (16)

область -

(17)

область -

(18)

Результирующие графики приведёны на рисунке 2.

Рисунок 2. Зависимости для .

Определим зависимости для напряжённостей магнитного поля, которые должны обеспечиваться предполагаемыми вторичными источниками поля (см. рисунок 3).

Рисунок 3. Зависимости для .

Эти напряжённости получены, учитывая выражения (7) и (16 - 18), а именно:

область - ; (19)

область - ; (20)

область - . (21)

Рассмотрим, как учитывают среду в виде ферромагнитной области , в соответствии с методом вторичных источников поля другие авторы [1]. Вводится дополнительная нормальная составляющая напряжённости магнитного поля на границе раздела сред. В области суммарная нормальная составляющая напряжённости поля равна , а в области - [1]. Учитывая граничные условия (3)

(22)

После несложных преобразований определяем

(23)

Проанализируем полученные результаты. Во-первых, видно из (19) и (20), что в точном решении краевой задачи, полученном аналитически, дополнительная составляющая нормальной напряжённости поля , направленная в немагнитную среду, по абсолютной величине не равна такой же составляющей, направленной в магнитную среду.

А также видно, что значение , соответствующее точному решению, не равно , введённому с помощью известного метода вторичных источников поля. Это видно из выражений (19), (20), (21) и (23).

Решим рассматриваемую задачу с помощью метода, предложенного в [1], применяя ряды Фурье.

В соответствии с этим методом на поверхности 2 сосредоточен поверхностный заряд , а на поверхности 3 - .

Для определения и рассмотрим рисунок 4 и применим теорему Гаусса [5]

(24)

где - вектор внешней нормали единичной поверхности;

- заряд, заключённый в области .

Рисунок 4. Поле вторичных источников.

Учитывая (23), (24) и рисунок 4, получим [1]

(25)

(26)

Уравнение (24) в дифференциальной форме будет иметь вид [5]

(27)

Применяя выражение (7), получим одномерное уравнение Пуассона

 (28)

Для рассматриваемого случая уравнение (28) принимает вид

(29)

где и - обобщённые функции, учитывающие поверхностное распределение зарядов (25), (26).

В данном случае мы вынуждены применить аппарат обобщённых функций, так как толщина области распределения поверхностных зарядов равна нулю и обычное разложение в ряды Фурье без рассмотрения предельных переходов невозможно.

Решение уравнения (29) будем искать в виде нечётного ряда Фурье

(30)

Подставим (30) в (29), умножим левую и правую части полученного выражения на и проинтегрируем от до .

Левая часть уравнения (29) примет вид

(31)

Правая часть уравнения (29) преобразуется к выражению

(32)

При вычислении интеграла (32) использовалась идея интегрального соотношения для обобщённых функций [6].

Приравняем (31) и (32) и, учитывая выражение найдём

(33)

На рисунке 5 представлено решение (линия 4) в виде ряда Фурье.

Рисунок 5. Сравнение численного и аналитического решений

Полученное решение рассчитано на ЭВМ с использованием выведенных формул (30) и (33). Сходимость ряда хорошая, уже при 20 членах ряда все углы острые. Явление Гиббса не наблюдается.

Видно, что наше решение, полученное аналитически, и решение по известному методу вторичных источников поля [1] не совпадают.

Из результатов, приведённых на рисунке 5 видно, что уравнения с помощью рядов качественно решаются корректно. Однако заряды и определены неправильно, что говорит о неправомерности выводов, приведённых в [1].

Попробуем решить эту проблему, учитывая намагничивающую среду следующим образом. В отличие от рисунка 4 и выражений (22) и (23) среда учитывается двумя разными векторами (воздух) и (ферромагнетик) по формуле

(34)

При учёте влияния границы знаки и учитываться не будут. Эти знаки получатся в процессе решения.

Из (34) следует

(35)

Используя выражения (6), (15) и (35) найдём распределение магнитного потенциала для вторичных источников поля в области

(36)

Аналогично найдём распределение магнитного потенциала для вторичных источников поля в области .

(37)

Аналогично получим результат и для области .

 (38)

Используя выражения (36), (37), (38) и свойство непрерывности магнитного потенциала

, при (39)

и , при , (40)

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

(41)

Решая систему (41), вычислим

(42)

(43)

Подставляя (42) и (43) в (36), (37) и (38), получим выражения для магнитного потенциала вторичных источников поля:

область - ; (44)

область - ; (45)

область - . (46)

Применяя (7) к (44) ? (46), получаем выражения для напряжённости магнитного поля от вторичных источников:

область - ; (47)

область - ; (48)

область - . (49)

Мы видим, что решения (44) - (49), полученные с помощью предложенной нами системы учёта ферромагнитной среды, полностью соответствуют аналитическому решению (16) - (21), полученному для контрольного примера классическим методом. Это подтверждает правильность предложенной нами системы учёта ферромагнитной среды.

Определяя и с помощью выражений (47) - (49) и применяя формулы (30) и (33) нами было найдено численное решение в виде ряда Фурье, которое полностью совпало с аналитическим решением, что также говорит в пользу предложенной системы.

Выводы. Нами показано, что существующая методика учёта неоднородности среды с помощью вторичных источников поля противоречит аналитическому решению, полученному классическим методом решения краевых задач. Нами предложена методика введения вторичных источников поля, позволяющая ликвидировать указанное противоречие. Показано, что возможно применение рядов Фурье в методе вторичных источников поля.

Список литературы

Бинс К. Анализ и расчёт электрических и магнитных полей: Пер. с англ./ К. Бинс, П. Лауренсон. - М.: Энергия, 1970. - 376 с.

Бессонов Л.А. Теоретические основы электротехники. - М.: Высш. шк., 1973. - 752 с.

Тозони О.В. Метод вторичных источников в электротехнике. - М.: Энергия, 1975. - 296 с.

Тозони О.В. Расчёт трёхмерных электромагнитных полей/ О.В. Тозони, И.Д. Маергойз. - Киев: Технikа, 1974. - 352 с.

Нейман Л.Р. Теоретические основы электротехники/ Л.Р. Нейман, К.С. Демирчян. - М.; Л.: Энергия, 1966. - Т. 2 407 с.

Колобов А.М. Избранные главы высшей математики. - Минск: Высш. шк., 1965. - 224 с.

Размещено на Allbest.ru