автореферат диссертации по электротехнике, 05.09.05, диссертация на тему:Метод электрических цепей Киргофа в задачах расчета стационарных магнитных полей

На правах^рсописи 0\

Байрамкулов Казим Нюрахматович

МЕТОД ЭЛЕКТРИЧЕСКИХ ЦЕПЕЙ КИРХГОФА В ЗАДАЧАХ РАСЧЕТА СТАЦИОНАРНЫХ МАГНИТНЫХ ПОЛЕЙ

05.09.05 — Теоретическая электротехника

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук

1 О ИЮН 2010

Новочеркасск - 2010 г.

004603931

Работа выполнена на кафедре «Прикладная математика» ГОУ ВПО «Южн Российский государственный технический университет (Новочеркасск политехнический институт)» и в лаборатории энергетики и электротехники Южно научного центра РАН

Научный руководитель: доктор технических наук, профессор

Астахов Владимир Иванович

Официальные оппоненты: доктор технических наук, профессор,

Басан Сергей Николаевич

кандидат технических наук, доцент Бурцев Юрий Алексеевич

Ведущая организация: Таганрогский технологический институт

ФГОУ ВПО «Южный федеральный университет»

Защита диссертации состоится «25» июня 2010 г. в 1400 ч. на заседай диссертационного совета Д 212.304.01 при ГОУ ВПО «Южно-Российск государственный технический университет (Новочеркасский политехничес институт)» в аудитории № 107 главного корпуса по адресу: 346428, г. Новочерка Ростовской обл., ул. Просвещения, 132.

С диссертацией можно ознакомиться в научно-технической библиотеке Г ВПО «Южно-Российский государственный технический универс! (Новочеркасский политехнический институт)». С текстом автореферата мо; ознакомиться на сайте ЮРГТУ (НПИ) www.npi-tu.ru .

Автореферат разослан «¿£» 2010 г.

Ученый секретарь диссертационного совета

Колпахчьян П.Г.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ Актуальность темы. Рост мощностей электротехнических устройств, тенденция уменьшения их габаритов и экономии дорогостоящих материалов, ведущая к режимам функционирования близким к предельным, остро ставят вопрос о повышении точности расчета электромагнитных полей в этих устройствах. Детальное исследование особенностей распределения поля в элементах устройств, невозможно без численных расчетов полей с использованием вычислительной техники. Поэтому актуальной проблемой теоретической электротехники является разработка эффективных методов моделирования и расчета магнитного поля, пригодных для сложных форм границ раздела сред с неоднородными, анизотропными и нелинейными свойствами.

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

Об актуальности использования электрических цепей Кирхгофа (ЭЦК) для расчета магнитных полей и характеристик электротехнических устройств свидетельствует значительное число диссертаций, монографий и статей, посвященных этой тематике, в том числе работы акад. АН СССР Л.Р. Неймана, акад. АН СССР и РАН К.С. Демирчяна, акад. РАН Я.Б. Данилевича, чл. корр. РАН П.А. Бутырина, А.Б. Новогородцева, М.А. Шакирова, Кияткина Р.П., С.М. Аполлонского, A.B. Иванова-Смоленского, Э.В. Колесникова, Э.А. Мееровича, JI.A. Цейтлина и др., работы профессоров ЮРГТУ (НПИ) (Л.Ф. Коломейцев, В.И. Астахов, А.Н. Ткачев, Г.К. Птах), а также работы зарубежных авторов (lin Hur, Hamid. A. Toliyat, Chun Y-D., Lee J. и др.). Но, несмотря на обилие работ в данном направлении, вопросы экономизации расчетов ЭЦК остаются актуальными. Целью работы является разработка эффективных алгоритмов метода ЭЦК для расчета стационарных магнитных полей с последующей реализацией на ЭВМ. Поставленная цель предполагает решение следующих задач:

1. Развить метод энергетического баланса (ЭБ) получения модели в виде ЭЦК как наиболее перспективного и обладающего возможностью оценки погрешности моделирования в терминах энергетических характеристик исследуемых устройств.

2. Создать на графах ЭЦК необходимый математический аппарат для формулирования и исследования аналогов краевых задач математической физики.

3. Адаптировать методы теории потенциала к графам ЭЦК.

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

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

1. Предложен способ моделирования методом ЭБ, отличающийся от известных

возможностью учета источников поля в элементе дискретизации. Предложен оригинальная технология получения модели в виде бесконечной ЭЦК методо энергетического баланса.

2. Впервые, на графах ЭЦК введен необходимый математический аппарат дл формулирования и исследования аналогов классических краевых зада математической физики. А именно, введены скалярные и векторные функци определены аналоги алгебраических операций и операций векторного анализ поставлены и исследованы на корректность в евклидовых пространствах аналог основных внутренних и внешних краевых задач математической физики дл уравнений Лапласа-Пуассона.

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

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

5. Создан пакет программ для расчета методом ЭЦК стационарных магнитны полей, позволяющий рассчитывать ЭЦК с бесконечным графом.

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

Практическая ценность работы. На основе полученных математических моделе и алгоритмов разработан пакет программ, предназначенный для расче стационарных и квазистационарных магнитных полей методом ЭЦК неоднородных, анизотропных или нелинейных средах с комплексн проницаемостью. Результаты работы используются во Всероссийском научн исследовательском и проектно-конструкторском институте электровозостроен (ВЭлНИИ) (г. Новочеркасск), Hl Iii «Эметрон», также пакет программ может най применение в электромашиностроении, трансформаторостроени приборостроении, и в вузах на специальностях «Прикладная математика «Теоретическая электротехника» и «Электромеханика» при выполнении курсовых дипломных работ.

Основные положения, выносимые на защиту:

1. Способ моделирования фрагмента области, удобный при наличии источник поля в этом фрагменте

2. Технология получения модели в виде бесконечной ЭЦК метод энергетического баланса

3. Математический аппарат для формулирования и исследования аналог классических краевых задач математической физики на графах ЭЦК

4. Адаптированные, к задаче дифракции, методы теории потенциала, позволяющ получать модели в виде СЛАУ минимальной размерности

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

6. Программный пакет, предназначенный для расчета стационарных магнитных

полей методом ЭЦК в неоднородных, анизотропных или нелинейных средах. Апробация работы. Теоретические положения и практические результаты диссертационной работы докладывались и обсуждались на: 54, 55 научно-технических конференциях профессорско-преподавательского состава, научных работников, аспирантов и студентов ЮРГТУ (НПИ) (Новочеркасск, 2005, 2006); V, VI школах-семинарах «Математическое моделирование, вычислительная механика и геофизика» для студентов, аспирантов и молодых ученых Юга России (Ростов-на-Дону, 2006, 2007); всероссийской выставке-ярмарке «ИННОВ-2007» (Новочеркасск, 2007); III и IV ежегодных конференциях студентов и аспирантов базовых кафедр ЮНЦ РАН - 2007, 2008 (Ростов-на-Дону, 2007, 2008); V всероссийской научно-практической конференции студентов, аспирантов и молодых ученых "Молодежь XXI века - будущее российской науки", диплом за второе место. (Ростов-на-Дону,

2007); Ш всероссийской школе-семинаре «Математические методы и биомеханика в современном университете» (пос. Дивноморское, 2007); международной конференции «Lyapunov Memorial Conference 2007» (Харьков, Украина, 2007); международной конференции «Теория операторов. Комплексный анализ и математическое моделирование» (Волгодонск, 2007); международных конференциях «52, 53, 54 Internationales Wissenschaftliches Kolloquium» (Ilmenau, Germany, 2007, ¿008, 2009); всероссийской конференции «Воронежская зимняя математическая школа С.Г. Крейна ВЗМШ-2008» (Воронеж, 2008); международной конференции «XII International Scientific Kravchuk Conference» (Киев, Украина,

2008); международных конференциях «Days on Diffraction'2008, 2009» (Санкт-Петербург, 2008, 2009); межрегиональной выставке «Информационные технологии в технике и образовании» (Новочеркасск, 2010).

Публикации. По теме диссертации опубликовано 18 печатных работ, общим объемом 5,5 печатных листа, в том числе 1 статья в рецензируемом научном журнале, рекомендованном ВАК.

Структура и объем работы. Диссертация состоит из введения, четырех глав, заключения и списка литературы. Основной объем работы составляет 141 стр. и включает 52 рис. Список цитируемой литературы содержит 54 наименования.

СОДЕРЖАНИЕ РАБОТЫ

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

Первая глава «Электрическая цепь Кирхгофа как модель среды, несущей поле», посвящена получению модели области с полем в виде ЭЦК. Рассмотрены традиционный метод и метод ЭБ получения модели в виде ЭЦК, проведено их сравнение. Предложен новый способ метода ЭБ позволяющий учитывать источники в моделируемом объекте.

При расчете электротехнических устройств, для повышения качества моделирования, проводится дискретизация расчетной области. Затем, исходя из условий задачи, моделируется каждый элемент дискретизации, объединением моделей для элементов дискретизации получают модель всей области. В работе, новым способом метода ЭБ, получается модель многополюсного участка - элемент дискретизации в общем виде. Предполагается, что на элементе AV:

rot# -0=> H = -gradtp, div£ = 0, B = \iH, B = Ybk, di\bt =-)—r

mes(AF)

и/тев^наЯ,,

Ькп=\ — --

[О на , 9 = 1, V, «у * А: = 1,у.

Затем рассматриваются интегралы

\\\\iidV = - \\\ь&ъ&ч>с1У =-г ДОфс^ - = ¡к - Фк), к = ¡Я

ду дк тез^дк) у

где <ру =-г ГГГфс/К, с другой стороны,

тез (А К)

\\\ЬкНс1У=\\\Ьк^Вс1У= = ¿ =

ДК ДГ дс р=1 Р=1

На рис. 1 приведена соответствующая ЭЦК.

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

На рис.2, изображены, для примера, сечения БК катушки несущей ток J (L -условная перегородка), и намагничивающих тел 5~ (5" = 5," и^").

Ч>1 Г11 Г12 I»

o-i-I^?—en----

ф2 Г21 Г22 Г2г |1

О-Ч—[сз^—■• —сэ—I

s* !

Рис.2

Рис.1

К типичным были отнесены следующие условия:

- сечения катушки и намагничиваемых тел занимают конечную часть плоскости и, таким образом, могут быть включены в круг SR некоторого радиуса R;

- первичные источники (токи катушки) и вторичные источники поля (микротоки) уравновешены (токи в сечениях SK и SJ, у = 1,2 равны нулю);

- среда, окружающая намагничиваемые тела (область S* на рис. 1.4, а), изотропна и однородна, то есть имеет магнитную проницаемость ц = ц+ = const;

- материал намагничиваемых тел характеризуется ограниченным, симметричным, положительно-определенным тензором ц~ второго ранга, причем его компоненты -функции положения точки в сечении S~.

{8 на 9 — — — —

к' 5=ёг5;сИу5 = 0;5 = цЯ; ji = JT = О вне SK,

на S ;

О)

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

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

~ц, О" О ц2

Н = Ц+ на5+; в; = В;, Щ =НГ на Г; \b(M)\-j^Q{M гг),

где учтено, что источники поля локализованы и уравновешены.

Задача (1) была названа задачей дифракции стационарного магнитного поля на плоскости 5 = S' u . Некоторые из свойств ее решения:

- напряженность Я0 первичного поля соленоидальна, то есть divH0 = 0 на S;

- напряженность поля реакции намагничиваемых тел Я* = Я - Н° потенциальна;

- на S* поле напряженности (Я) и индукции (В) отличаются лишь масштабом;

- в силу соленоидальности индукции имеем ^Hndl = 0 для VL с S*;

l

- благодаря уравновешенности источников, выполняется [|Яс// = 0 для VZ,, вне SR;

L

- в слое единичной толщины, основанием которого является плоскость S, энергия магнитного поля конечна, то есть <00.

s

На практике, как правило, задачу (1) преобразуют к скалярному виду, вводя потенциал <р напряженности или функцию потока у индукции магнитного поля.

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

В работе рассматривается преобразование задачи для магнитного поля на плоскости к ЭЦК использованием функции потока. Полагается, что катушки с током имеют магнитную проницаемость ц = ji+, и нас интересует поле вне их сечения. В такой постановке свойство соленоидальности индукции В магнитного поля, входящее в (1), обеспечивается представлением В через функцию потока у, то есть В = [grady, ёг]. Затем проводится триангуляция плоскости равномерным разбиением на треугольники вида AS,,AS2 (см. рис.3). Из соображений простоты, считается, что Н =ёхНх +ёуНу; Нх,Ну— const на AS,, и, что магнитная проницаемость в пределах отдельного элемента дискретизации является постоянным тензором второго ранга

О "I

ц = , то есть = const, \а)У = const на AS,. Отмечается также, что

. Руу}

однородная напряженность на элементе дискретизации не может иметь источников

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

Рис.3 Рис.4 Рис.5

В работе, рассматривается возможность моделирования, методом ЭБ, элементов дискретизации с меняющейся магнитной проницаемостью (см. рис.4). Практическая значимость последнего - возможность моделирования областей с малыми зазорами без дополнительного учащения дискретизации. К примеру, для элемента дискретизации приведенного на рис.4 получается модель, показанная на рис.5.

К каждому элементу дискретизации, для получения модели в виде ЭЦК, применяется метод ЭБ. Сшиванием моделей для элементов триангуляции, получается модель всего объекта в целом, показанная на рис.6.

Вторая глава «Краевые задачи на графе электрической цепи», посвящена созданию на графах ЭЦК необходимого математического аппарата для формулирования и исследования аналогов классических внутренних и внешних краевых задач математической физики.

Рассматривается пространственный или плоский направленный граф,

состоящий из ветвей, /--г—т-д---Г—"1—-д_____

соединённых в контуры. Бесконечным, называется граф, содержащий

бесконечное число ветвей и узлов — точек соединения ветвей. Узлы Мк, к = 1,2,.. конечного или

бесконечного графа О нумеруются, а их множество ]

обозначается М. Для ветвей вводится двойная нумерация, указанием номеров узлов, к которым

она присоединена, и Рис.6

Ж тя(А10)

ориентации. Например, обозначение означает, что данная ветвь соединена с узлами Мк, Мр и ориентирована единичным вектором , направленным вдоль ветви от Мк к Мр. Если парам узлов соотнесены две или более ветви, чтобы различать, их пометили дополнительными символами, например, /¿,,/J, и т.д. Множество ветвей графа обозначено L.

На Г2 определяются скалярные {<р} и векторные {а} дискретные функции. Каждому узлу Мк соотносится число <рк (функция времени <рк (t)), а объединение последних обозначается ф и называется скалярной дискретной функцией ф или, короче , скаляром 9 на Q. Каждой ветви соотносится число а^ (функция времени akp{t)), затем берется акр = l^a^, их объединение обозначается а и

называется векторной дискретной функцией а или, короче, вектором а на Q . Числовые значения скалярных дискретных функций в узлах и векторных — на ветвях любой конечной части графа полагаются ограниченными. В случае, когда значения дискретных функций - есть функции времени, полагается, что множество точек на временной оси, в которых указанные значения бесконечны, имеют меру нуль.

Равенства ф = 0 и а = 0 понимаются в том смысле, что все значения данных дискретных функций равны нулю, поэтому 0 и 0 уместно называть нулевыми, соответственно, скалярной и векторной дискретной функцией. Если значения ф одинаковы во всех узлах графа и равны константе С (функции времени C(t)), то такая функция называется постоянной и обозначается как ф = С. Затем

_ _ (ф* + Ф J

определяются операции фу = ф^, в V Мк е М, фа = а^-^ , на V е L,

__ i

ab = ~2_1 aqpbqp> в VMj б М . Через s обозначается некоторое сечение графа, то есть 2 Ч

множество точек замкнутой поверхности S, ориентированной единичной внешней нормалью «, являющихся одновременно точками графа, и определяется поток вектора а через s

равенством ^ ña = ^ ~ñl¡^atp, где s под знаком

S I

суммы означает, что суммирование ведется по точкам ветвей j/^j, принадлежащим сечению s, а

произведение = ±1. Дивергенция векторной дискретной функции a на О определяется как diva = ^Гйд, в VMq е М. Примеры сечений s и s? Рис.7

показаны на рис.7. Окрестностью узла Mq называется совокупность ветвей, сходящихся в нем, и узлов, соединенных этими ветвями с Mq, а также введены

обозначения: М?- множество узлов в окрестности узла Мц (на рис.7 выделены дополнительным кружком); множество ветвей в окрестности узла Мд. В этих

обозначениях получим, что = в УМЦ еМ(Мр еМ5). Градиент

ч

скалярной дискретной функции ф на О определяется ка ^ск? = Т^Ыу = Т°(ур -фДнаУ/^ еЬ.

В этих обозначениях на П выполняются тождества ётра =§гаёфа + ф(11уа ^(Цуа = из которых следуют ][](gradфgrach|^+фД\|/) = ^йфgrad\y

Ма I М8 I

^(фAy-^|/Дф) = ^и(фgrad\i/-vlígгadф), аналоги первого и второго тождеств Грина

М5 I

здесь М8 под знаком суммы означает, что суммирование ведется по узлам принадлежащим части графа О, охваченной поверхностью Б, порождающе сечение е. символом А, для краткости, был обозначен оператор divgrad, то ест Аф = d¡vgradф = - Ф,), в \/Мч е М . /„

Располагая на графе оператором А, можно рассматривать аналоги уравнений Лапласа-Пуассона, уравнения , р

теплопроводности, волнового уравнения и др., а на фрагментах графа - аналоги классических краевых задач

математической физики для этих уравнений в совпадающих обозначениях. еШуг^^

В работе рассматриваются р типичные варианты краевых условий, в) имеющие практическое значение для электрических сеток. В частности, Рис.8

электрический потенциал ф? узлов сетки, обозначенных на рис.8, а темны

кружком, при заданных е1,е2,...,еи известен. А именно,

%=еч> ? = (условия Дирихле), (

На рис.2.5, б при заданных известны

~— = г1 , ^ = 1,2,(условия Неймана), (

дпч

если фрагмент графа электрической сетки, отвечающий ее плоской части, отнести О". Для сетки, изображенной на рис.8, в, в тех же узлах выполняется И. Я(ч"

ф" + ——— = е , ¡7 = 1,2,(смешанное условие), (

г опч

а для рис.8, г -

1 dgr = 1 Эф+ r~ dnq г+ dnq'

ДФ = 0 вне fiR; = |ф(л/?)|-^->0(1/ра)

<7 = 1,2.....N (условия склейки), (5)

где ± в индексах соотносятся фрагментам Q* графа электрической сетки. Будем различать (классифицировать) краевые задачи по типу основных краевых условий.

При исследовании и решении краевых задач удобно рассматривать исходные и искомые дискретные функции как элементы некоторых евклидовых пространств. В этой связи на конечном связном графе £2 рассматривается класс С(М) скалярных дискретных функций с ограниченными значениями на М. Для постановки и исследования внутренних краевых задач вводятся следующие пространства:

М LMM

\1/2

/2>0(М): ve/2(м) и 2> = 0, /'„(М):^), ^GradyGradiH, J Yfimdyf

М L 4L

Для постановки и исследования внешних краевых задач рассматривается класс С(М) ограниченных дискретных функций, подчиняющихся условиям

1 Эф

г* дп

Для внешних краевых задач вводятся следующие пространства:

4(М): (V4)= X vi |H| = (w)1/2,/2.0(M): ^/2(М),^0наГк, =

М^иГ, MBr

T2l(M): (V47 =^Gradi)/Gradi ||V||3 ,^eC(M)

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

Теорема 1. Решение внутренней задачи Дирихле для уравнения Пуассона единственно.

Теорема 2. Внутренняя задача Неймана для уравнения Пуассона

= на М", ~-=f, М еГ, <7 = 1,2,...,N, (6)

dnq

разрешима не всегда. Необходимым условием разрешимости является выполнение

£f+2/=O, (7)

М- г

Если последнее имеет место, то любая пара решений может отличаться лишь на постоянную.

Теорема 3. Решение внутренней задачи для уравнения Пуассона с условием

смешанного типа при Re

/ N,

г

~R„

V ?

>0 единственно.

Теорема 41. Необходимым условием разрешимости задачи Аф = -Б на М~ У М+,

1 5Ф+ 1 ар- _

дп г~ дп

г ап г

на конечном графе П=О- и У Г является равенство

и* ' м- ' г

Если оно выполняется, то решение задачи (8) при

Яе

г -\ г

> 0 и (или) 1т

(9)

(10)

определяется с точностью до постоянной.

Теорема 5. Пусть ф- гармоническая дискретная функция на конечном фрагменте О* с границей Г некоторого графа. Тогда

Г* дп

(П)

Теорема 6 (о среднем). Пусть дискретная функция ф подчиняется уравнению Лапласа в узле Мд е М. Тогда в этом узле она равна средневзвешенному своих значений в узлах, составляющих окрестность Шд. А именно,

я>, =-Тт2>(?>*)ф*> М* еМ?> 02)

где п{д)~ число ветвей в ; п(с},к)-число ветвей, соединяющих Мч с Мк. Теорема 7 (о максимуме и минимуме). Пусть ф- отличная от постоянной, гармоническая на связном фрагменте графа, дискретная функция. Тогда наибольшее (наименьшее) значение ф не может достигаться в узлах данного фрагмента. Теорема 8. Решение уравнения Пуассона в С(М) на конечном графе О. при свободном члене из 120 существует и определяется с точностью до постоянной дискретной функции.

Теорема 9. Задача Дирихле для уравнения Пуассона с однородным краевым условием на конечном фрагменте О1 графа О при ограниченном свободном члене корректна в пространстве 1\ а.

Теорема 10. Задача Дирихле с неоднородным краевым условием для уравнения Пуассона на конечном фрагменте О1 графа П при ограниченных исходных данных имеет в 1\ решение и притом единственное.

1 В данной теореме связность Ь" или I/ не предполагается.

Теорема 11. Решение внутренней задачи Неймана для уравнения Пуассона (6) в 1\ при выполнении условия (7) существует и определяется с точностью до постоянной дискретной функции.

Теорема 12. Задача Ду = -И на , + — = / на Г для уравнения Пуассона

сп

со смешанным граничным условием на конечном фрагменте графа при ограниченных дискретных функциях Б и /, вещественных и положительных значениях дискретной функции % на Г разрешима и притом единственным образом в ¡2.

Теорема 13. Задача дифракции Ду = -Б на М\Г, а.,,- ига Г, на

а——= дп дп

конечном графе С2 со смежной границей Г его фрагментов £2+, О" при ограниченных дискретных функциях Б на М\Г, /, е на Г, а также Яе(а/р)>0 разрешима в классе дискретных функций с конечной суммой Дирихле, причем ее решение определяется с точностью до произвольной постоянной функции на М .

Теорема 14. Решение уравнения Пуассона в С(М) Ду = -У7, е /2 0 (М) на правильном бесконечном графе существует, единственно и представимо в форме объемного потенциала

у(м?)= £ Р{мру{мч,мр),м^ш (13)

м п„\м.

Теорема 15. Внешние задачи Дирихле для уравнений Лапласа и Пуассона в постановках

Д^ = -^наМ\ Аф = 0 на М+,

\ = Ск наГ4, ^ = = Щ Ф = У-С,наГ„Х^ = 0,А = Щ

Ту сп оп

разрешимы в /2'0(М) и притом единственным образом при любой ограниченной у на Г и любом свободном члене из

Третья глава «Методы теории потенциала на графе ЭЦК», посвящена адаптации методов теории потенциала к задаче дифракции на графе ЭЦК. что позволило рассчитывать ЭЦК с бесконечным графом, решая конечномерное уравнение.

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

соединяющей узлы Мч и Мр еМ, М- множество неустранимых узлов цепи (на

рис.6 выделены светлыми кружками). Затем, исходя из первого закона Кирхгофа, ставится краевая задача

ЕЫр)-^)} 1-— -г+у(?),мгем,им„

м, \ гчр )

->0|1

(15)

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

Используется фундаментальное решение Я(р,д) уравнения Лапласа на О., определенное условиями

(16)

при фиксированном узле Мр е М (узел Мр. вынесен в б.у.у).

Дискретная функция в терминах ЭЦК имеет ясную физическую

интерпретацию. А именно, на ЭЦК с единичными сопротивлениями ветвей рассматриваемого графа Я(р,д) есть потенциал узла Мч в условиях, когда в узел Мр «впрыскивается» извне единичный ток, а потенциал последнего узла

принимается равным единице по соображениям удобства.

Рассматривается дискретный аналог объемного потенциала уравнения Лапласа на графе О. А именно,

м

Прямая подстановка (17) в уравнение = еМ, показывает, что

последнее удовлетворяется. А если носитель М^ функции Р(д) имеет конечно число узлов и

Е^М, а«

то выполнится условие в б.у.у. В работе показывается, что правая часть уравнени (14) обладает свойством (18). Согласно изложенному функция

г ^

У(я) = -Ек(<1'Р)

ЕМа)-у(Р)}

(19

Мч &М,Мр еМрМл тМр,Мк еМ_,, удовлетворяет всем условиям задачи (14), (15

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

г+

1--

+£*(<?,¿УД*), (20

Му

Мд еМ,Мр е М5,М„ еМр,Мк еМ,. С другой стороны, сужение (20) на М^ есть СЛАУ

Т = АЧ'+В, (21)

( г А

где =^{Л(д,р)-Цд,Н)} 1-— ,ь9 =£*(?,?,Р = 1,2..(Шп(Ма).

м, \ гр1, ) Му

Поскольку матричный оператор СЛАУ конечномерен, то вопрос о ее разрешимости сводится к вопросу о единственности решения.

В работе показывается, что рассматриваемая задача и СЛАУ эквивалентны. Следовательно, доказанная ранее единственность решения задачи (14), (15) обеспечивает единственность решения СЛАУ (21), а значит, и их разрешимость. Таким образом, доказана

Теорема 16 Решение задачи (14), (15) в С(М) (введенном для исследования внешних краевых задач) на правильном бесконечном графе существует, единственно и представимо в виде (19).

В случае кусочно-однородной изотропной среды представления (19), (20) существенно упрощаются и получается СЛАУ минимальной размерности.

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

Наиболее сложными и важными с точки зрения приложений являются нелинейные задачи. В работе рассматривается случай, когда гистерезис отсутствует, но материал сердечника обладает нелинейными свойствами. В работах Маергойза И. Д. рассмотрены несколько вариантов итерационного метода расчета поля в такой среде и обоснована их сходимость. Существенно, что оценки и неравенства из работ Маергойза И. Д., лежащие в обосновании сходимости, сохраняют силу и в том случае, если обновление по кривой намагничивания значений магнитной проницаемости выполнять не на каждом шаге уточнения решения СЛАУ. Это означает, что итерационный процесс может быть разделен на два подпроцесса. В первом из них магнитная проницаемость фиксирована и решается СЛАУ, во втором - по найденному с достаточной точностью решению СЛАУ выполняется коррекция значений магнитной проницаемости.

Отсюда следует ряд полезных выводов. Во-первых, медленно сходящиеся последовательные приближения при решении СЛАУ можно заменить значительно более быстрым алгоритмом метода исключения Гаусса. Во-вторых, способ получения самой СЛАУ не имеет принципиального значения, поэтому вместо использовавшегося в работе Маергойза И.Д. допустимо применение метода ЭЦК-

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

р = ц(|#|). Способ использования этих зависимостей по кривой намагничивания

представлен на рис.9.

Я*

в tgaw = ц("> /

я<">, wl("4

вУ

0 W? и

а) б)

Рис.9

а), б) - соответственно, области сильного, слабого насыщения. Причем, разделение магнитопровода на эти участки оказывается достаточным проделать по результатам первой итерации, которая выполняется при ц = const . Как показали численные эксперименты при нарушении указанной рекомендации наблюдается знакопеременный характер приращений магнитной проницаемости на каждом шаге итераций. Отметим также, что при пользовании ЭЦК, на каждом шаге нами рассчитываются одновременно Н и В. При использовании указанного правила получим изменения магнитной проницаемости, показанные на рис.10, слева - в точке М,, справа — в точке М2 (см. рис. 11).

Четвертая глава «Программная реализация», посвящена созданному программному пакету, в ней приводится его краткое описание и внешний вид программных модулей. Отладка и контроль программного пакета и метода проводились на модельной задаче и сравнением с существующими пакетами программ по расчету стационарных магнитных полей. В качестве контрольного | примера, предназначенного для отладки программного пакета, использовалась плоская квадратная область с магнитным диполем, висящим над ее центром на высоте, много меньшей длины стороны квадрата. В такой системе можно считать, что диполь висит над идеально проводящим полупространством. В такой! постановке функцию потока индукции магнитного поля легко определить, аналитически методом зеркальных отображений.

В,Тл|

0.5

2 3

В,Тл

0.15

0.05

250

450 H, А/м

Рис.10

V

80 Н, Мм

Щ

Рис.11

Рассчитанные с помощью программного пакета (на графике под номером 1) и аналитически (под номером 2) значения функции потока на верхней стороне прямоугольника изображены на рис.12. На участке, представленном на рис.12, полученные результаты зрительно не различимы. При численном решении использовалось разбиение верхней стороны на 440 элементов, при этом относительная невязка численного и аналитического решений, не превысила 1.01 %.

V'i0\ Вб/м

Рис.12 Рис.13 Рис.14

Также, количественный и качественный контроль созданного программного пакета, проводился с помощью широко известного конечно-элементного пакета БЕММ 4.0. Полученные результаты сравнивались, как в случае кусочно-однородных сред, так и в случае нелинейности, относительная погрешность не превысила 4%.

На рис.13, для наглядной демонстрации масштабирования приведен расчет электромагнита и, слева направо, два шага масштабирования. Примеры расчетов, для области с малым зазором и нелинейной характеристикой материала, приведены на рис.14.

ОСНОВНЫЕ РЕЗУЛЬТАТЫ РАБОТЫ

1. Развит метод ЭБ - математически строгого получения модели в виде ЭЦК. А именно, предложена оригинальная технология получения модели в виде бесконечной ЭЦК методом энергетического баланса. Предложен новый способ моделирования методом ЭБ, удобный при наличии источников поля в моделируемом объекте. Проведено сравнение методов получения схем замещения, традиционного и метода ЭБ. Рассмотрены типичные для практики электромагнитных расчетов варианты краевых задач, используемые при нахождении стационарного магнитного поля. Показано использование функции потока и потенциала в преобразовании задачи для магнитного поля на плоскости к ЭЦК.

2. На графах ЭЦК введен необходимый математический аппарат для формулирования и исследования аналогов классических краевых задач математической физики. А именно, введены скалярные и векторные функции, определены аналоги алгебраических операций и операций векторного анализа, поставлены и исследованы на корректность в евклидовых пространствах аналоги основных внутренних и внешних краевых задач математической физики для уравнений Лапласа-Пуассона.

3. Применительно к задаче дифракции на графе ЭЦК адаптированы методы теории потенциала позволяющие свести краевую задачу дифракции на графе к матричному уравнению минимальной размерности. Обоснована применимость

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

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

5. На основе полученных математических моделей и алгоритмов разработан пакет программ «KCNM», предназначенный для расчета стационарных и квазистационарных магнитных полей методом ЭЦК в неоднородных, анизотропных или нелинейных средах с комплексной проницаемостью. Предусмотрена возможность масштабирования фрагмента области. Он может найти применение в электромашиностроении, трансформатороспроении, приборостроении, а также в вузах на специальностях «Прикладная математика», «Теоретическая электротехника», «Электромеханика» при выполнении курсовых и дипломных работ.

ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ По материалам диссертационного исследования лично и в соавторстве

опубликовано 18 научных работ (общим объемом 5,5 печатных листа), в том числе

одна в ведущем научном журнале, рекомендованном ВАК. В опубликованных

научных работах достаточно полно отражены основные научные результаты.

Основные работы, опубликованные по теме диссертации: Издания рекомендованные ВАК:

1. Байрамкулов К. Н.-А., Астахов В. И. Расчет магнитного поля в среде с неоднородными и анизотропными свойствами на основе электрической цепи Кирхгофа//Изв. вузов. Электромеханика. - 2010. - №1. - С. 3-11. (1.04/0.73) (вклад соискателя: технология получения ЭЦК как модели для среды, с неоднородными и анизотропными свойствами, несущей магнитное поле)

Остальные:

2. Байрамкулов К.Н.-А., Астахов В.И. Расчет стационарного магнитного поля методом электрических цепей КирхгофаУ/Студенческая научная весна - 2005: Сб. науч. тр. аспирантов и студентов ЮРГТУ (НПИ)//Новочеркасск ЮРГТУ. -2005. - С. 234-236. (0.13/0.07) (вклад соискателя: разработка математической модели в виде ЭЦК и реализация алгоритмов вычислений)

3. Астахов В.И., Байрамкулов К. Н.-А. Внутренние краевые задачи на графе электрической цепи//Математические методы в физике, технике и экономике. -Новочеркасск: Редакция журнала «Изв. вузов. Электромеханика», 2006. - С. 337. (1.99/1.69) (вклад соискателя: исследование внутренних и внешних краевых задач на графах ЭЦК)

4. Байрамкулов К.Н.-А., Астахов В.И. Решение задачи дифракции методом граничных уравнений на электрических сетках//Студенческая научная весна — 2006: Материалы 55-й науч. - техн. конф. студентов и аспирантов ЮРГТУ (НПИ). - Новочеркасск, 2006. - С. 98-100. (0.15/0.075) (вклад соискателя: разработка математической модели в виде ЭЦК и реализация алгоритмов вычислений)

5. К.Н.-А. Байрамкулов, В.И. Астахов О масштабировании фрагмента области при расчете магнитного поля в кусочно-однородной среде методом теории цепей//Математическое моделирование и информационные технологии. -Новочеркасск: Редакция журнала «Изв. вузов. Электромеханика», 2007. - С. 5760. (0.22/0.18) (вклад соискателя: разработка математической модели в виде ЭЦК и реализация алгоритмов вычислений)

6. Байрамкулов К. Н.-А., Астахов В.И. Расчет магнитного поля методом граничных уравнений на графе электрической цепи//Труды Южного научного центра РАН. - Ростов-на-Дону: Изд. ЮНЦ РАН, 2007. - Т. 2. С. 72-79. (0.61/0.45) (вклад соискателя: разработка математической модели в виде ЭЦК и реализация алгоритмов вычислений)

7. Байрамкулов К.Н.-А. Математическое моделирование и расчет статических и стационарных полей методом граничных уравнений на графах электрических цепей// Математическое моделирование, вычислительная механика и геофизика. Труды V Школы-семинара, Ростов-на-Дону, 18-21 декабря 2006, г. Ростов-на-Дону: Изд. «ЦВВР», 2007. - С. 41-43. (0.18)

8. К. N.-A. Bayramkulov, V.l. Astakhov The method of the boundary equations in problems of computing static and stationary fields on the topological graph (Метод граничных уравнений в задачах расчета статических и стационарных полей на топологических графах)//52. Internationales Wissenschaftliches Kolloquium Technische Universität Ilmenau, September 10 - 13. 2007. P. 169-170. (0.13/0.07) (вклад соискателя: разработка математической модели в виде ЭЦК и реализация алгоритмов вычислений)

9. К. N.-A. Bayramkulov Mathematical modeling and computing of static and stationary fields by a method of the boundary equations on columns of electric circuits(MaTeMaTH4ecKoe моделирование и расчет статических и стационарных полей методом граничных уравнений на электрических цепяху/LYAPUNOV MEMORIAL CONFERENCE. International Conference on the occasion of the 150th birthday of Aleksandr Mikhailovich Lyapunov: Book of abstracts. - Kharkiv: Verkin Institute for Low Temperature Physics and Engineering of NASU. 2007. - P. 10-11. (0.06)

10. K. N.-A. Bayramkulov The modeling of fragment of area with magnetic field by Kirchhoff electric circuit network (Моделирование электрической цепью Кирхгофа фрагмента области с полему/Information Technology and Electrical Engineering - Devices and Systems, Materials and Technology for the Future: 54 Internationales Wissenshaftiches Kolloquium, Germany, Ilmenau, 7-10 September 2009/Techniche Universität Ilmenau. - Ilmenau, 2009. - P.195-196. (0.08) Личный вклад соискателя в работы, опубликованные в соавторстве, состоит в:

технологии получения ЭЦК как модели для среды, с неоднородными и анизотропными свойствами, несущей магнитное поле [1]; исследовании внутренних и внешних краевых задач на графах ЭЦК [3]; разработка математической модели в виде ЭЦК и реализация алгоритмов вычислений [2,4,5, 6, 8].

Байрамкулов Казим Нюрахматович

Метод электрических цепей Кирхгофа в задачах расчета стационарных магнитных полей

Автореферат

Подписано в печать 21.05.2010. Формат 60x84 1/16 Бумага офсетная. Печать цифровая. Усл. печ. л. 1,0. Уч.-изд. л. 1,6. Тираж 130 экз. Заказ «Я-АШ

Отпечатано в ИД «Политехник» 346428, г. Новочеркасск, ул. Просвещения, 132 тел., факс (863-5) 25-53-03

ВВЕДЕНИЕ

Глава 1 Электрическая цепь Кирхгофа как модель среды, несущей поле

1.1 Методы моделирования электрическими цепями

1.1.1 Традиционный метод

1.1.2 Метод энергетического баланса

1.2 Постановка задачи расчета магнитного поля в присутствии намагничиваемых тел

1.3 Функция потока в преобразовании задачи для магнитного поля на плоскости 22 1.3.1 Моделирование малых зазоров

1.4 Потенциал в преобразовании задачи для магнитного поля на плоскости

1.5 Потенциал поля реакции в преобразовании задачи для магнитного поля на плоскости

1.6 Особенности моделирования ЭЦК в трехмерной постановке

1.7 Сравнение методов 37 Выводы по главе

Глава 2 Краевые задачи на графе электрической цепи

2.1 Скалярные, векторные дискретные функции и операции над ними

2.2 Уравнения и тояадества на графе и его фрагментах

2.3 Краевые условия и внутренние краевые задачи 54 и/ 2.4 Пространства дискретных функций

2.5 Свойства гармонических дискретных функций

2.6 Фундаментальное решение уравнения Лапласа и первая краевая задача

2.7 Вторая и третья краевые задачи

2.8 Задача дифракции

2.9 Внешние краевые задачи

Выводы по главе

Глава 3 Методы теории потенциала на графе ЭЦК

3.1 Краевая задача на графе ЭЦК

3.2 Фундаментальное решение уравнения Лапласа

3.3 Переход от краевой задачи к матричному уравнению II рода

3.4 Обоснование применимости метода последовательных приближений

3.5 О способе вычисления фундаментального решения

3.6 Масштабирование как способ экономизации

3.7 Случай нелинейных свойств материала сердечника 115 Выводы по главе

Глава 4 Программная реализация

4.1 Краткое описание пакета программ

4.2 Контроль разработанного программного пакета

4.3 Примеры расчетов 128 Выводы по главе 4 133 ЗАКЛЮЧЕНИЕ 134 СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ

Непрерывный рост мощностей электротехнических устройств, тенденция уменьшения их габаритов и экономии дорогостоящих материалов, ведущая к режимам функционирования близким к предельным, остро ставят вопрос о повышении точности расчета электромагнитных полей в этих устройствах. Детальное исследование особенностей распределения поля в элементах устройств, невозможно без численных расчетов полей с использованием вычислительной техники [1—3]. Поэтому актуальной проблемой теоретической электротехники является разработка эффективных методов моделирования и расчета магнитного поля, пригодных для сложных форм границ раздела сред с неоднородными, анизотропными и нелинейными свойствами. На сегодняшний день создан целый ряд численных методов анализа магнитного поля: методы интегральных уравнений, конечных разностей, конечных элементов, граничных элементов, декомпозиционнные методы, методы сведения полевых задач к цепным задачам.

Использование схем замещения при выполнении инженерных расчетов делает последние более наглядными, однако длительное время электрические цепи Кирхгофа (ЭЦК) служили лишь математической основой для физического моделирования полей различной природы на так называемых электрических сетках [4-8]. Соответствующие ЭЦК, находили, используя аналогию между конечно-разностным представлением дифференциальных уравнений в частных производных и соотношениями между током и напряжением на элементах ЭЦК [4, 5, 9]. Но, появившиеся в последнее время мощные пакеты прикладных программ по расчету сложных электрических и электронных схем, вкупе с бурным развитием компьютерной техники, стимулируют использование ЭЦК в качестве математической модели физических полей, предназначенной для компьютерной реализации. Об этом свидетельствуют работы российских, к примеру в [10—15], и зарубежных авторов, таких как [16-18]. Вместе с тем, при решении типичных для практики электромагнитных расчетов задач методом ЭЦК, зачастую получаются цепи с десятками или сотнями тысяч ветвей и узлов, или даже с их бесконечным числом, расчет которых остается проблемой для современных ЭВМ, а попытки снизить её остроту представляются актуальными.

Об актуальности построения новых схемных моделей для расчета магнитных полей и характеристик электротехнических устройств свидетельствует значительное число диссертаций, монографий и статей, посвященных этой тематике, в том числе работы акад АН СССР и РАН К.С. Демирчяна, акад. АН СССР JI.P. Неймана, акад. РАН Я.Б. Данилевича, А.Б. Новогородцева, М.А. Шакирова, Кияткина Р.П., С.М. Аполлонского, чл. корр. РАН П.А. Бутырина, А.В. Иванова-Смоленского, Э.В. Колесникова, Э.А. Мееровича, JI.A. Цейтлина, В.А. Филина, Г.Н. Цицикяна и др., работы профессоров ЮРГТУ (НПИ) (Коломейцев Л.Ф., В.И. Астахов, Г.К. Птах, А.Н. Ткачев) , а также работы зарубежных авторов (Jin Hur, Hamid. A. Toliyat, Jung-Pyo Hong, Chun Y-D., Lee J., Wakao S., M.V.K. Chari, H.B. Dwight, R.L. Ferrari, B. Hague, H. Kaden, C.I. Mocanu, P.P. Silvester, M. Stafl, R.L. Stoll, J. Turowski и др.). Несмотря на обилие работ в данном направлении вопросы экономизации расчетов ЭЦК остаются актуальными.

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

Цель данной работы состоит в экономизации математического аппарата для расчета стационарного магнитного поля на основе ЭЦК, с последующей реализацией на ЭВМ. Поставленная цель предполагает решение следующих задач: 1. Развить метод ЭБ [10, 11] получения модели в виде ЭЦК как наиболее перспективного и обладающего возможностью оценки погрешности моделирования в терминах энергетических характеристик исследуемых устройств.

2. Создать на графах ЭЦК необходимый математический аппарат для формулирования и исследования аналогов краевых задач математической физики.

3. Адаптировать методы теории потенциала к графам ЭЦК.

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

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

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

В третьей главе рассмотрены методы теории потенциала на графе ЭЦК. Предложен способ получения фундаментального решения уравнения Лапласа и приведено представление дискретной функции на графе через фундаментальное решение. Выполнено сведение краевой задачи дифракции на графе к матричному уравнению второго рода, обоснована применимость метода последовательных приближений к его решению. Описано масштабирование как способ экономизации расчета ЭЦК. Рассмотрен случай нелинейных свойств материала сердечника.

В четвертой главе описан созданный на основе построенной теории программный пакет, приведены результаты контрольных расчетов. Также рассмотрены примеры моделирования технических задач с использованием созданного программного пакета.

Основные результаты диссертационной работы опубликованы в [19—35]. Она была апробирована на следующих конференциях:

1. Конференция студентов и аспирантов ЮРГТУ (НПИ) 2005 и 2006 годов.

2. V Школа-семинар «Математическое моделирование, вычислительная механика и геофизика» для студентов, аспирантов и молодых ученых Юга России, декабрь 2006 г., г. Ростов-на-Дону, Россия.

3. «52. Internationales Wissenschaftliches Kolloquium», сентябрь 2007 г., г. Ильменау, Германия.

4. Y Всероссийская научно—практическая конференция студентов, аспирантов и молодых ученых "Молодежь XXI века - будущее российской науки", май

2007 г., г. Ростов-на-Дону, Россия (диплом за второе место).

5. Третья ежегодная конференция студентов и аспирантов базовых кафедр ЮНЦ РАН, апрель 2007 г., г. Ростов-на-Дону, Россия (диплом за второе место).

6. III Всероссийская школа-семинар «Математические методы и биомеханика в современном университете», мая -июнь 2007 г., пос. Дивноморское, Россия.

7. Lyapunov Memorial Conference 2007, июнь 2007 г., г. Харьков, Украина.

8. Теория операторов. Комплексный анализ и математическое моделирование, август-сентябрь 2007 г., г. Волгодонск, Россия.

9. VI Школа-семинар «Математическое моделирование, вычислительная механика и геофизика» для студентов, аспирантов и молодых ученых Юга России, декабрь 2007 г., г. Ростов-на-Дону, Россия.

10. Days on Diffraction'2008, июнь 2008 г., г. Санкт-Петербург, Россия.

11. Воронежская зимняя математическая школа С.Г. Крейна ВЗМШ-2008, январь

2008 г., г. Воронеж, Россия.

12. «53. Internationales Wissenschaftliches Kolloquium», сентябрь 2008 г., г. Ильменау, Германия.

13. Четвертая ежегодная конференция студентов и аспирантов базовых кафедр ЮНЦ РАН - 2008, апрель 2008 г., г. Ростов-на-Дону, Россия {диплом за третье место).

14. Days on Diffraction'2009, май 2009 г., г. Санкт-Петербург, Россия.

15. «54. Internationales Wissenschaftliches Kolloquium», сентябрь 2009 г., г. Ильменау, Германия.

Основной объем работы составляет 141 стр. и включает 52 рис. Список цитируемой литературы содержит 54 наименования.

Основные результаты работы формулируются в следующем виде:

1. Развит метод ЭБ — математически строгого получения модели в виде ЭЦК. А именно, предложена оригинальная технология получения модели в виде бесконечной ЭЦК методом энергетического баланса. Предложен новый способ моделирования методом ЭБ, удобный при наличии источников поля в моделируемом объекте. Проведено сравнение методов получения схем замещения, традиционного и метода ЭБ. Рассмотрены типичные для практики электромагнитных расчетов варианты краевых задач, используемые при нахождении стационарного магнитного поля. Показаны особенности использования функции потока и потенциала в преобразовании задачи для магнитного поля на плоскости к ЭЦК.

2. На графах ЭЦК введен необходимый математический аппарат для формулирования и исследования аналогов классических краевых задач математической физики. А именно, введены скалярные и векторные функции, определены аналоги алгебраических операций и операций векторного анализа, поставлены и исследованы на корректность в евклидовых пространствах аналоги основных внутренних и внешних краевых задач математической физики для уравнений Лапласа-Пуассона.

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

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

5. На основе полученных математических моделей и алгоритмов разработан пакет программ «KCNM», предназначенный для расчета стационарных и квазистационарных магнитных полей методом ЭЦК в неоднородных, анизотропных или нелинейных средах с комплексной проницаемостью. Предусмотрена возможность масштабирования фрагмента области. Разработанный пакет программ контролировался на тестовой задаче и сравнением с известными пакетами. Он может найти применение в электромашиностроении, трансформаторостроении, приборостроении, а также в вузах на специальностях «Прикладная математика», «Теоретическая электротехника», «Электромеханика» при выполнении курсовых и дипломных работ.

ЗАКЛЮЧЕНИЕ

В диссертационной работе предложен эффективный алгоритм метода ЭЦК для расчета стационарных магнитных полей в неоднородных, анизотропных или нелинейных средах.

1. Демирчян К.С., Чечурин В.Л. Машинные расчеты электромагнитных полей. М.: Высш. шк., 1986. 240 с.

2. Астахов В.И. Математическое моделирование инженерных задач в электротехнике. Новочеркасск: Новочерк. гос. техн. ун-т, 1994 с. 192.

3. Демирчян К.С., Бутырин П. А. Моделирование и машинный расчет электрических цепей. М.: Высш. шк., 1988. 335 с.

4. Гершгорин С.А. Об электрических сетках для приближенного решения дифференциального уравнения Лапласа. // Журнал прикладной физики. 1929, т. 6, вып. 3-4, стр. 3-29.

5. Волынский Б. А. и Бухман В. Е. Модели для решения краевых задач. Физматгиз, М., 1960. 451 с.

6. Чудаков А.Д. Электрические моделирующие сетки и их применение. М.: Энергия,1968. 136 с.

7. Шимони К. Теоретическая электротехника /Пер. с нем. М.: Мир, 1964. 773 с.

8. Гутенмахер Л.И. Электрические модели. К.: Техшка, 1975. 176 с.

9. Р. Курант, К. Фридрихе, Г. Леви, О разностных уравнениях математической физики//УМН, 1941, № 8, С.125-160.

10. Астахов В.И. Моделирование цепями Кирхгофа электротехнических устройств. Режим переменного тока // Изв. вузов. Электромеханика, 1999. № 3. С.105-120.

11. Астахов В.И. Моделирование цепями Кирхгофа электротехнических устройств. Режим постоянного тока //Изв.вузов. Электромеханика. 1998. № 5-6. С. 95-108.

12. Иванов-Смоленский А. В., Кузнецов В. А. Универсальный численный метод моделирования электромеханических преобразователей и систем

13. Электричество. 2000, N 7, с. 24-33

14. Шакиров М.А., Кияткин Р.П. Схемы замещения для систем левитации и тяги напеременном токе//Электротехника- 1999- №8 С. 11-19.

15. Jin Hur, Hamid. A. Toliyat, Jung-Pyo Hong Dynamic analysis of linear induction motors using 3-D equivalent magnetic circuit network (EMCN) method/ZElectric Power Components and Systems, vol. 29, pp. 531-541, 2001.

16. J. Perho, Reluctance network for analyzing induction machines Ph.D. Dissertation, Dept. Electrical Eng., Helsinki University of Technology, Finland, 2002.

17. Chun Y-D., Lee J., Wakao S., Magnetic Field Analysis Using 3D Equivalent Magnetic Circuit Network Method// Papers of Technical Meeting on Static Apparatus, IEE Japan, vol.SA-02, pp.83-88, 2002.

18. Байрамкулов K.H.-A., Астахов В.И. Расчет магнитного поля методом граничных уравнений на графе электрической цепи // Труды южного научного центра РАН.-2007. Том 2. С. 72-79.

19. Байрамкулов К.Н.-А., Астахов В.И. Расчет стационарного магнитного поля методом электрических цепей Кирхгофа//Студенческая научная весна 2005: Сборник научных трудов студентов и аспирантов ЮРГТУ (НПИ)/ Новочеркасск ЮРГТУ,2005. - С.234-236.

20. Байрамкулов К.Н.-А., Астахов В.И. Решение задачи дифракции методом граничных уравнений на электрических сетках // Студенческая научная весна -2006: Материалы 55-й науч. техн. конф. студентов и аспирантов ЮРГТУ (НПИ), 2006. - С.98-100.

21. Байрамкулов К.Н.-А. Метод граничных уравнений в задачах расчета статических и стационарных полей на топологических графах // Третья ежегодная научная конференция студентов и аспирантов базовых кафедр

22. Южного научного центра РАН: Тезисы докладов (Ростов-на-Дону, 5-24 апреля 2007г.). Ростов-на-Дону: Изд-во ЮНЦ РАН, 2007. С 114-115.

23. Байрамкулов K.H.-A. Масштабирование фрагмента области при расчете магнитного поля в кусочно-однородной среде методом теории цепей // Труды воронежской зимней математической школы С.Г. Крейна 2008. Воронеж: ВорГУ, 2008. С. 41-44.

24. Байрамкулов К. Н.-А., Астахов В. И. Расчет магнитного поля в среде с неоднородными и анизотропными свойствами на основе электрической цепи Кирхгофа / Известия вузов. Электромеханика. 2009. №6

25. К. N.-A. Bayramkulov The modeling of fragment of area with magnetic field by kirchhoff electric circuit network // Days on difraction'2009. International conference. Saint Petersburg, May 26-29, 2009. Abstracts. P. 17.

26. Астахов В.И. О вариационном методе расчета магнитных полей // Изв. вузов. Электромеханика. 2003. №2. с.З—17.

27. Астахов В.И., Кочубей Т.В., Шапошников К.С. Метод ортогональных проекций в задачах расчета стационарных магнитных полей // Труды южного научного центра РАН. -2007. Том 2. С. 51-72.

28. Астахов В.И. О допустимости идеализации границ поляризуемых тел и некоторых энергетических соотношениях для стационарного магнитного и электростатического полей // Изв. вузов. Электромеханика. 2000. №1. с.3-14.

29. Тамм И. Е. Основы теории электричества. М.: Наука, 1989. 504 с.

30. Ректорис К. Вариационные методы в математической физике и технике. М.: Мир, 1985. 590 с.

31. Кочин Н.Е. Векторный анализ и начала тензорного исчисления / М.: Изд. АН СССР, 1961.426 с.

32. Фихтенгольц Г.М. Основы математического анализа. М.: Наука, Т.2, 1964. 440с.

33. Ахиезер Н.И., Глазман И.М. Теория линейных операторов в гильбертовом пространстве. Харьков: Высш. шк., 1977. 315 с.

34. Полна Г., Сеге Г. Изопериметрические неравенства в математической физике. М.: Физматиз,1962. 336 с.

35. Михлин С.Г. Линейные уравнения в частных производных. М.: Высш. шк., 1977. 431 с.

36. Г. Нуссбаумер Быстрое преобразование Фурье и алгоритмы вычисления сверток /Пер. с англ. М.: Радио и связь, 1985. 246с.

37. Владимиров В7СГУравнения математической~физики. М.: Наука, 197675281^

38. Тихонов А.Н., Самарский А.А. Уравнения математической физик. М.: Наука, 1977. 735 с.1. ЗУ- (цу

39. Маергойз И. Д. Итерационные методы расчета статистических полей в неоднородных, анизотропных и нелинейных средах. Киев, Наук, думка, 1979. 210с.

40. Бахвалов Н. С. Численные методы. М.: Наука, 1973. 632с.

41. Калиткин Н.Н. Численные методы. М.: Наука, 1978. -512с.

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

43. Андре Анго Математика для электро- и радиоинженеров. М.: Наука, 1965. 780с.

44. Бронштейн И.Н., Семендяев К.А. Справочник по математике для инженеров и учащихся втузов. М.: Наука, 1986. 544 с.