Scientific journal
Modern problems of science and education
ISSN 2070-7428
"Перечень" ВАК
ИФ РИНЦ = 0,791

MODELLING OF PICTURES OF ELECTROMAGNETIC FIELDS ON THE BASIS OF TRANSIENTS

Vasileva O.V. 1
1 National Research Tomsk Polytechnic University
Работа посвящена сведению задач решения дифференциальных уравнений в частных производных к хорошо разработанным методам исследования нелинейных динамических систем, которые позволяют сопоставить уравнению в частных производных динамическую систему с последующим ее моделированием с помощью электротехнических устройств, аналого-цифровых комплексов. Рассматривается моделирование картин электромагнитных полей на основе переходных процессов без наличия ложных скачков при построении линий равного тока в схемотехнической среде MATLAB-Simulink. Новизна алгоритма заключается в сведении дифференциальных уравнений в частных производных Лапласа и Пуассона к динамической системе – системе нелинейных дифференциальных уравнений 1-го порядка на основе электротехнических схем замещения. Составлена структурная схема динамической системы, которая отображает динамику процесса, где элементами схемы являются операционные усилители, на основе которых может быть сопоставлена электротехническая схема замещения.
Work is devoted to the information of problems of the decision of the differential equations in private derivatives to well developed investigative techniques of nonlinear dynamic systems which allow to compare the equation in private derivatives with dynamic system to its subsequent modeling by means of electrotechnical arrangements, analog-digital complexes. Modeling pictures of electromagnetic fields on the basis of transients without availability of false jumps is considered at construction of lines of an equal current in environment of the circuitry MATLAB-Simulink. The novelty of the algorithm consists in the reduction of differential equations in partial derivatives of the Laplace and Poisson for a dynamical system – the system of nonlinear differential equations of the 1-st order on the basis of electrotechnical equivalent circuits. The block diagram of dynamic system which displays dynamics of process where elements of the scheme are operational amplifiers on the basis of whom the electrotechnical equivalent circuit can be compared is made.
Poisson´s equation.
differential equations
transient
MATLAB-Simulink
the block diagram
dynamic system
an electromagnetic field
lines of an equal current
modeling

Введение

Для визуального восприятия картин электромагнитных полей используется графическое построение пространственного распределения поля – как функции двух переменных . Это обычно пространственное распределение скалярного или векторного потенциалов полей [4-6]. Известно, что с помощью программ MathCAD и MATLAB в среде pdetool можно построить картину поля лишь для линий равного потенциала, причем программы, позволяющие построить картину поля для линий равного тока, отрисовывают контурные линии при наличии ложных скачков [1] (рис. 1, 2).

Это объясняется тем, что линии тока линейного заряда или линии потока линейного тока распределяются лишь в смысле главного значения, то есть в пределах от до (рис. 1а, 2а, 3а).

Рис. 1. Метод решения в частных производных для электростатического поля: а) фазовый портрет; б) пространственное распределение линий равного тока

Рис. 2. Метод решения в частных производных для магнитостатического поля: а) фазовый портрет; б) пространственное распределение линий равного потока

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

Теоретический анализ

Рассмотрим алгоритм для построения фазового портрета динамической системы в схемотехнической среде MATLAB-Simulink.

Рис. 3. Геометрическая интерпретация распределения линий тока линейного заряда: а) для уравнений в частных производных; б) для динамической системы

Новизна алгоритма заключается в сведении дифференциальных уравнений в частных производных Лапласа и Пуассона к динамической системе, то есть к решению системы дифференциальных уравнений 1 порядка в среде MATLAB-Simulink.

Алгоритм представляется в следующей последовательности:

1. Ввод данных: количество, координаты, величина зарядов (токов), количество изолиний равного потенциала (тока):

a) для электростатического поля;

b) для магнитостатического поля.

2. Решение дифференциальных уравнений в частных производных Лапласа и Пуассона, используя входные данные:

a) для линий равного потенциала;

b) для линий равного тока или потока.

3. Сведение решения дифференциальных уравнений в частных производных Лапласа и Пуассона к динамической системе.

4. Решение системы дифференциальных уравнений 1-го порядка, полученных с помощью динамической системы, в среде MATLAB-Simulink.

5. Формирование интерфейса блока входных данных:

a) выбор типа поля для решения системы дифференциальных уравнений 1-го порядка;

b) выбор количества координат, величин зарядов (токов);

c) выбор количества изолиний равного потенциала, тока (потока).

6. Выход: построение фазового портрета динамической системы.

Методика

Исследование статических электрического и магнитного полей будет осуществляться с помощью динамической системы и последующим ее моделированием аналоговой техникой в схемотехнической программе MATLAB-Simulink.

Под динамической системой мы понимаем электромеханическую систему, которая описывается дифференциальными уравнениями [2; 3].

Каждая динамическая система имеет свой фазовый портрет, которому соответствуют свои особые точки (это точки положения равновесия).

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

Из методов исследования можно выделить следующие:

1. Методы дифференциальной геометрии и теории автоматического регулирования.

2. Методы математической физики – уравнения Лапласа, Пуассона.

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

, (1)

где - комплексный электрический потенциал, - функция потенциала, - функция потока, - объемная плотность зарядов, - электрическая постоянная, - диэлектрическая проницаемость.

Комплексный потенциал описывает совокупность силовых и эквипотенциальных линий поля, которые ортогональны по отношению друг к другу, то есть картину поля.

Решение уравнения Пуассона (1) для электростатического поля:

, (2)

где - величина заряда, - регуляризирующий коэффициент.

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

Для плоского электростатического поля имеем систему уравнений:

, (3)

где - плотность зарядов, при решении приняли .

Уравнение (3) означает потенциальность поля:

. (4)

Постоянство потенциала поля при изменении координат описывается уравнением:

. (5)

Из уравнения (5) получаем в параметрической форме систему уравнений:

. (6)

Теперь получим уравнения для линий напряженностей электрического поля, для этого рассмотрим уравнение (4):

(7)

Постоянство линий потока при изменении координат описывается уравнением:

. (8)

Из уравнения (8) получаем в параметрической форме систему уравнений:

. (9)

Таким образом, имеем две динамические системы:

, . (10)

Экспериментальная часть

Пусть заданы декартовые координаты трех линейных зарядов:

.

Используя уравнение (10), получаем обыкновенные дифференциальные уравнения для трех линейных зарядов, полученные с помощью динамической системы (рис. 5а):

, (11)

. (12)

На рисунке 4 представлено всплывающее окно блока входных данных, позволяющее пользователю, не зная программирования, получить на выходе интересующую картину статического электромагнитного поля в фазовой плоскости и на основе частных структурных схем (рис. 4а, б) сопоставить электротехническую схему замещения [1].

Блок входных данных состоит из следующих пунктов:

  1. Тип поля: электростатическое поле, магнитостатическое поле.
  2. Количество зарядов / токов: 1-6.
  3. Координаты по оси x: -3 – 5.
  4. Координаты по оси y: -4 – 6.
  5. Величина зарядов / токов: -3 – 3.
  6. Количество изолиний равного потенциала: 5, 10, 13, 20, 23, 30.
  7. Количество изолиний равного тока: 5, 10, 13, 20, 23, 30

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

Картина поля становится наиболее наглядной и информативной при построении линий равного тока, которые позволяют судить о траектории движения заряженных частиц в электростатическом поле. Распределение линий равного потенциала и линий равного тока взаимоперпендикулярны в пространстве. Образованная этими линиями площадь криволинейных прямоугольников, показывающих распределение адиабатического инварианта, в частности энергии электростатического поля и заряда, имеет величину емкости С. Густота и разреженность линий позволяют судить об интенсивности распределения энергий в пространстве.

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

Аналогичные рассуждения можно сделать и по магнитному полю [1].

Рис. 4. Всплывающее окно блока входных данных и картина поля в фазовой плоскости (д): а) количество зарядов; б) координаты и величина заряда; в) количество изолиний равного потенциала; г) количество изолиний равного тока

Рис. 5. Фазовый портрет, полученный с помощью динамической системы

Результаты

В данной статье приведено описание метода сведения решения дифференциальных уравнений в частных производных Лапласа и Пуассона к решению системы обыкновенных нелинейных дифференциальных уравнений первого порядка с помощью динамической системы.

Заменены реальные процессы аналоговой техникой через операционные усилители с помощью программы схемотехники MATLAB-Simulink. Получен фазовый портрет электромагнитных полей на основе моделирования динамической системы аналоговой техникой.

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

Рецензенты:

Лукутин Б.В., д.т.н., заведующий кафедрой ЭПП ЭНИН ФГБОУ ВПО «НИ ТПУ», Национальный исследовательский Томский политехнический университет, г. Томск.

Канев Ф.Ю., д.ф.-м.н., профессор кафедры ЭСиЭ ЭНИН ФГБОУ ВПО «НИ ТПУ», Национальный исследовательский Томский политехнический университет, г. Томск.