Bulletin of student science of the Department of Information Systems and Programming
Scientific Journal
RUSENG

Engineering
DEVELOPMENT OF THE SOFTWARE FOR CALCULATION VECTOR MAGNETIC POTENTIAL OF THE SECOND PART ELECTROMAGNETIC COILS
Popov Boris Klavdievich 1, Popova Olga Borisovna 2, Panesh Davlet Muratovich 3

1. Department of Power Supply of Industrial Enterprises of the Kuban State Technological University
2. Department of Information Systems and Programming of the Kuban State Technological University
3. Student in Masters

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

Наиболее важным аспектом первоначального исследования проектируемого устройства является его математическая модель. Так ка удачный выбор её приводит к успешному результату проектирования.

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

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

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

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

В результате развития вычислительной техники появились системы автоматизированного проектирования (САПР) различных технических устройств. С развитием САПР стало широко применяться математическое моделирование.Современные ЭВМобладают памятьюбольшого объёма, большим быстродействием и развитой периферией. В результате развития САПР появилисьновыеаналитические и численные методы расчета электромагнитных и тепловых полей.

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

Одними из первых аналитических методов являются методы Роговского и Рота. Здесь решение получается в виде одинарных или двойных рядов Фурье.Метод применим для прямоугольных границ и для магнитной системы с µ=const.

Метод конформных преобразований является более действенным аналитическим методом, чем методы Роговского и Рота, так как позволяет учесть влияние границ более сложной формы,чем прямоугольник, но он применим только для двухмерных полей и магнитных систем с µ=0.

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

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

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

Следует также отметить метод вторичных источников поля, являющийся развитием членом-корреспондентом Гринбергом Г. А. метода, предложенного ещё Максвеллом. Метод позволяет решать практически любые задачи технической электродинамики, но требует большого объема оперативной памяти машины и большого быстродействия, а также глубоких математических знаний.

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

На наш взгляд наиболее интересен метод, объединяющий аналитические и численные методы. Ранее, одним из авторов был применён метод вторичных источников поля с использованием рядов Фурье и метода Фурье для решения задач в двухмерных полях [2 ‒ 14]. Как показала вычислительная практика при числе членов ряда около 100, явление Гиббса не наблюдалось, и результаты расчётов были близки к точным результатам. Это стало понятно после сравнения результатов расчёта по предлагаемой методике с результатами расчёта контрольных примеров.

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

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

Решение, как и для двухмерного случая, будем искать для следующей конструкции. В сплошной среде с µ=const образуем полость, заполненную воздухом, в виде параллелепипеда. В середину этого параллелепипеда поместим цилиндрическую катушку. Наша задача будет заключаться в том, что надо будет рассчитать электромагнитное поле в любой точке пространства.

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

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

Векторный магнитный потенциал  выразим в виде тройного ряда Фурье

В результате проведённых исследований нами была выведена формула для определения

где  − плотность тока в катушке;

        l, d − размеры полости;

       mk, nk, ut ‒ номера гармоник.

 

Формула (4) может служить для расчёта поля в диаметральной плоскости, параллельной оси z.Мы ограничились определением GZh,k,t , так как была поставлена задача расчёта векторного магнитного потенциала второй части электромагнитной катушки.

Из формулы (4) видно, что интеграл от z не имеет аналитического решения. Поэтому для определения данного интеграла применим формулу в виде тройного шаблона. Исследование использования тройного шаблона в нашей задаче показало его эффективность и точность, так как удалось решить проблему численного интегрирования периодических функций. Решение оказалось в пределах допустимой точности.

С использованием приведённых формул и описанного выше подхода были составлены алгоритмы, а по ним программы вычисления GZh,k,t [15, 16]. В результате проведённых вычислений было выявлено, что работа программы занимает незначительное время, не превышающее пары минут. Картина поля, построенная по результатам вычислений, полностью соответствует физическим представлениям, описанным в литературе. Отсюда вывод, что данная методика подходит для решения трёхмерных задач в неоднородных средах и для любых видов катушек.


Bibliography:
1. Popov B.K. Matematicheskoye modelirovaniye magnitnykh sistem gruzozakhvatnykh ustroystv na postoyannykh magnitakh: monografiya/ B.K. Popov, O.B. Popova, A.I. Popov; FGBOU VO «KubGTU», Krasnodar, 2017, 100 s.
2. Popov B.K. Issledovaniye primenimosti metoda vtorichnykh istochnikov prireshenii krayevykh zadach elektromagnitnogo polya/ Kuban. gos. tekhnol. un-t. – Krasnodar, 2003. – 15 s. – Dep. v VINITI 11.04.03, № 698 – V2003.
3. Popov B.K. Issledovaniye vliyaniya razmerov neodnorodnosti sredy na velichinu vtorichnykh istochnikov pri reshenii krayevykh zadach elektromagnitnogo polya/ Kuban. gos. tekhnol. un-t. – Krasnodar, 2003. – 9 s. – Dep. v VINITI 06.06.03, № 1103 – V2003.
4. Popov B.K. Opredeleniye zavisimosti vtorichnykh istochnikov polya ot srednego znacheniya napryazhonnosti polya na granitse razdela sred / Kuban. gos. tekhnol. un-t. – Krasnodar, 2004. – 4 s., Dep. v VINITI 08.06.2004, № 970 – V2004.
5. Popov B.K. Issledovaniye primenimosti ryadov Fur'ye v metode vtorichnykh istochnikov polya dlya resheniya ploskoparallel'noy neodnorodnoy polevoy zadachi/ Kuban. gos. tekhnol. un-t. – Krasnodar, 2004. – 6 s., Dep. v VINITI 12.10.04, № 1592 – V2004.
6. Popov B.K. Vyvod formuly dlya opredeleniya velichiny i znaka vtorichnykh istochnikov polya/ Kuban. gos. tekhnol. un-t. – Krasnodar, 2004. – 4 s., – Dep. v VINITI 08.12.04, № 1962 – V 2004.
7. Popov B.K. Issledovaniye primenimosti integrala Fur'ye v metode vtorichnykh istochnikov polya dlya resheniya neodnorodnoy ploskoparallel'noy polevoy zadachi/ Kuban. gos. tekhnol. un-t. – Krasnodar, 2004. – 5 s., – Dep. v VINITI 15.12.04, № 2001 – V2004.
8. Popov B.K. Privedeniye beskonechnogo chisla posledovatel'nykh priblizheniy k konechnomu analiticheskomu vyrazheniyu pri reshenii polevoy zadachi s pomoshch'yu integrala Fur'ye/ Kuban. gos. tekhnol. un-t. – Krasnodar, 2005. – 7 s., – Dep. v VINITI 10.10.05, № 1290 – V2005.
9. Popov B.K. Raschot metodom Rogovskogo magnitnogo polya dvukh shin pryamougol'nogo secheniya so vstrechno napravlennymi tokami v pryamougol'noy oblasti, okruzhonnoy ferromagnetikom s beskonechnoy magnitnoy pronitsayemost'yu/ Kuban. gos. tekhnol. un-t. – Krasnodar, 2007. – 11 s.: il., Dep. v VINITI 23.03.07, № 285 – V2007.
10. Popov B.K. Raschot metodom Rotamagnitnogo polya dvukh shin pryamougol'nogo secheniya so vstrechno napravlennymi tokami v pryamougol'noy oblasti, okruzhonnoy ferromagnetikom s beskonechnoy magnitnoy pronitsayemost'yu/ Kuban. gos. tekhnol. un-t. – Krasnodar, 2007. – 7 s.: il., Dep. v VINITI 23.03.07, № 287 – V2007.
11. Popov B.K. Vyvod formuly dlya opredeleniya velichiny i znaka vtorichnykh istochnikov polya v vide poverkhnostnykh tokov/ Kuban. gos. tekhnol. un-t. – Krasnodar, 2007. – 7 s., Dep. v VINITI 23.03.07, № 286 – V2007.
12. Popov B.K. Raschot magnitnogo polya dvukh shin pryamougol'nogo secheniya so vstrechnonapravlennymi tokami v kusochno-neodnorodnoy pryamougol'noy oblasti, okruzhonnoy ferromagnetikom s beskonechnoy magnitnoy pronitsayemost'yu/ Kuban. gos. tekhnol. un-t. – Krasnodar, 2007. – 9 s., Dep. v VINITI 02.05.07, № 485 – V2007.
13. Popov B.K. Resheniye polevykh zadach elektrotekhniki s pomoshch'yu vtorichnykh istochnikov polya i ryadov Fur'ye/ B.K. Popov, O.B. Popova; Politematicheskiy setevoy elektronnyy nauchnyy zhurnal Kubanskogo gosudarstvennogo agrarnogouniversiteta (Nauchnyy zhurnal KubGAU) [Elektronnyyresurs]. – Krasnodar: KubGAU, 2013. – №09(093). – IDA [article ID]: 0931309040. – Rezhim dostupa: http: //ej.kubagro.ru/2013/09/pdf/40.pdf.
14 Popov B.K. Uchot neodnorodnosti sredy pri raschote magnitnogo polya/ B.K. Popov, O.B. Popova; Politematicheskiy setevoy elektronnyy nauchnyy zhurnal Kubanskogo gosudarstvennogo agrarnogo universiteta (Nauchnyy zhurnal KubGAU) [Elektronnyy resurs]. – Krasnodar: KubGAU, 2014. – №01(095). – IDA [article ID]: 0951401052. – Rezhimdostupa: http: //ej.kubagro.ru/2014/01/pdf/52.pdf.
15. Popov B.K. Vychisleniye chlena ryada Fur'ye dlya vektornogo magnitnogo potentsiala katushki po osi Z/ B.K. Popov, O.B. Popova; SVIDETEL'STVO o gosudarstvennoy registratsii programmy dlya EVM 2014619337, zaregistrirovano15.09.14.
16. Popov B.K. Modifitsirovannoye vychisleniye chlena ryada Fur'ye dlya vektornogo magnitnogo potentsiala katushki po osi Z/ B.K. Popov, O.B. Popova; SVIDETEL'STVO o gosudarstvennoy registratsii programmy dlya EVM 2014619542, zaregistrirovano 18.09.14.

Reference

Popov Boris Klavdievich, Popova Olga Borisovna, Panesh Davlet Muratovich DEVELOPMENT OF THE SOFTWARE FOR CALCULATION VECTOR MAGNETIC POTENTIAL OF THE SECOND PART ELECTROMAGNETIC COILS // Вестник студенческой науки кафедры информационных систем и программирования. – 2017. – № 03;
URL: vsn.esrae.ru/en/3-11 (Date Access: 29.06.2025).


Embed on your website or blog

Viewed articles

Today: 2097 | Week: 2098 | Total: 3040


Comments (0)