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

кандидата технических наук
Дарьин, Сергей Георгиевич
город
Томск
год
1991
специальность ВАК РФ
05.09.01
Автореферат по электротехнике на тему «Математическое моделирование и автоматизация расчетов магнитных полей электрических машин с произвольной конфигурацией магнитопровода»

Автореферат диссертации по теме "Математическое моделирование и автоматизация расчетов магнитных полей электрических машин с произвольной конфигурацией магнитопровода"

• р \ к. *

Томский ордена ОктяОрьской Революции и ордена Трудового Красного Знамени политехнический институт им. С.Ы.Кирова

На правах рукописи Дарьин Сергей Георгиевич

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

Специальность 05.09.01 - Электрические машины

Автореферат диссертации на соискание ученой степени кандидата технических наук

Томск 1991

Работа выполнена во Всесоюзном научно-исследовательском институте электромеханики (ВНШЭМ) и Всесоюзном научно-исследовательском проектно-конструкторском и технологическом институте электромашиностроения (ВНЖШШ).

доктор технических наук, профессор СтрельСицкий Э.К.

доктор технических наук, профессор Хорьков К. Л.

кандидат технических наук, Горисев С.А.

СКБ "Генератор", г. Москва

Защита состоится ц декабря 1991 года в 15 .00 часов на заседании специализированного совета К 063.80.01 по присуждению ученой степени кандидата технических наук Томского ордена Октябрьской Революции и ордена Трудового Красного знамени политехнического института имени С.М. Кирова в Актовом зале главного корпуса (634004, г.Томск, пр. Ленина, 30).

С диссертацией можно ознакомиться в библиотеке института.

Автореферат разослан ноября 1991 г.

Научный руководитель -

Официальные оппоненты -

Ведущее предприятие

Учений секретарь специализированного

совета, канд. техн. наук, доцент .. А.И.Чучалин

ЮМ^^--

■'•-.'Г..

* I

:.;./ 3

'! ОБЩАЯ ХАРАКТЕРИСТИК РАБОТЫ

• • '

АКТУАЛЬНОСТЬ ТЕ.!Ы. Современная экономическая ситуация, боз-никакдая при переходе к рынку заставило? разработчиков и производителей научно-технической продукт::!, в тем числе и в электротехнической промшленности, осваивать ковке виды изделий, сокращать сроки их разработки и освоения производства, всем! способам:! снижать езок издержки. В стране оау^зегся резкая нехватка различного рода электрических машин. Расширяется сфера их применения, становятся более разнообразными требования к техническим характеристикам, экономически« и зкеплуатадиотгкм показателям. Дгя удовлетворения этих требования все ча^з приходится искать нетрадиционные технические репэния.

Чтобы сократить сроки разработки и уменьшить затраты на до-■ водку необходимо еще ка стадии проектирования как молно более точно определить характеристики разрабатываемых устройств. Традиционные инженерные методики расчета, основы которых были заложены в 30 - 50-х годах, теряют свое значение. Эти методики, имеющие существенные ограничения по применению, . довольно часто дают результаты, расходящиеся с результатами испытаний даже по'основным показателям, что объясняется недостаточна.! учетом таких, факторов как реальная конфигурация капштапроводь, нелинейные магнитные свойства материалов, дискретная структура обмотск к т.д.

Повьтаиие гЭДяктивностя преектяо-кояегг.укторехях работ ноле? быть достигнуто только путс:< внедрения в практику проектирования совреке-нют методов математического моделирования процессов, имеющих место в электрических машинах. Широкому использованию таких методов в проекта« организациях препятствуя? слабость математического аппарата, требующая высокой квалификации пользов^елл, большой об-ьем рутшпюй раСота по подготовке исходна данных и анализу получаем»: результатов. Креме того, •эти метода для своей реализации предъявляй? ыиохио требования к вычислительной теиит-э.

В последнее время прооктно-кокструкторским организациям е?алп доступны соЕренвшше ».ятя- и горсэЕэяыте ЭК?.;, сблздаюгцке вкдагкч быстродействием, йолькш об-ьемхм оперативной и внешней мамяги, ги-рскгг-1 выбором яереферайгых устройств (графических дисплеев с логовкма устройствами, гр^слостроятслея, :,!пз.чич:ча принтеров). Использование такой техники открывает хороним церспек-гива по пош-пению точности и азтематкзации математического моделирспакия.

Учитывая ыгае иодтятк?, кежно гсойстатисовн-гь, что гадами'

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

ЦЕЛЬ РАБОТЫ. Целью работы является разработка на основе метода конечных элементов (МКЭ) алгоритмов к универсальных пакетов прикладных программ для моделирования электромагнитных процессов в электрических мааинах с разнообразной конфигурацией активной области.

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

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

разработать' алгоритмы и программные средства для структур ного моделирования разнообразных электромашинных объектов в дву мерной и трехмерной постановке с использованием интерактивной машинной графики;

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

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

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

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

1,!ЕТ0ДЫ ИССЛЕДОВАНИЯ. Основным методом является численш эксперимент. При разработке, алгоритмов использовались злемен* теории электромагнитного поля, вычислительные методы алгебры, ; риационное исчисление. Физические эксперименты выполнялись в I

ганизациях-разработчиках электрических машин.

НАУЧНАЯ НОВИЗНА работы:

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

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

3. Создана база типовых фрагментов (элементов) для формирования моделей различных электрических малин.

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

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

6. Разработаны интерактивные процедуры обработки результатов моделирования.

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

ПРАКТИЧЕСКАЯ ЦЕННОСТЬ работы заключается в том, что разработанные на базе метода конечных элементов пакеты прикладных программ позволяют проводить точные расчета и анализ электромагнитных процессов в электрических малинах с разнообразной конфигурацией активной части, в том числе оригинальных конструкций, для которых отсутствуют методики расчета. Кроме того, пакеты могут быть использованы для обоснования принимаемых допущений и корректировки инженерных методик. Высокая автоматизация процессов подготовки данных и анализа результатов при использовании пакетов дает возможность пользователю сосредоточиться на сути задачи, не требуя глубокого знания теории метода и вычислительной техники. Структура пакетов и принятая организация данных позволяет использовать пакеты в САПР различных электрогехшиеских устройств в качестве подсистемы высокоточных расчетов. Высокая наглядность пакетов позволяет использовать их в учебном процессе в вузах при изучении численных методов расчета поля в электрических машинах.

РЕАЛИЗАЦИЯ РЕЗУЛЬТАТОВ РАБОТЫ. Разработанные пакеты программ расчета магнитного поля использовались во Всесоюзном научно-

исследовательском проектно-конструкторском и технологическом институте электромашиностроения при проведении научно-исследовательских работ по теме "Совершенствование методов расчета и методов оптимального проектирования активной части асинхронных 'двигателей" и Всесоюзном научно-исследовательском институте электромеханики по теме "Создание криогенных электрических машин переменного тока для судового электропривода". Упомянутое выие программное обеспечение внедрено и использовалось в Научно-исследовательском центре "Завода имени Владимира Ильича", НПО РОТОР г. Черкассы, ряде других предприятий и вузов в исследовательских и учебных целях.

АПРОБАЦИЯ РАБОТЫ. Основные положения и научные результаты диссертационной работы докладывались и обсуждались на: заседании секции № 2 Научного Совета АН СССР по теме "Состояние вопроса разработки программных продуктов по расчету электромагнитных полей" (г. Ленинград, 1937 г.), VIII Всесоюзной научно-технической конференции "Состояние и перспективы совершенствования разработки, производства и применения низковольтных электродвигателей переменного тока" (г. Владимир, 1986 г.).

ПУБЛИКАЦИИ. По теме диссертации опубликовано 8 печатных работ.

СТРУКТУРА И ОБЪЕМ ДИССЕРТАЦИИ. Диссертационная работа состоит из введения, пяти разделов, заключения, содержит 136 страниц основного текста, 37 рисунков, 8 таблиц, включает список использованной литературы из 155 наименований и приложения на 33 страницах.

ОСНОВНЫЕ ПОЛОЖЕНИЯ, ПРЕДСТАВЛЕННЫЕ К ЗАЩИТЕ. .

1. Алгоритмы и пакеты прикладных программ для численного моделирования магнитных полей в электрических машинах с произвольной конфигурацией магнитопровода.

2. Решение с помоац>:о разработанных пакетов широкого круга двумерных и трехмерных полевых задач.

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

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

ВО ВВЕДЕНИИ обоснована актуальность теш, сформулирована цель работы, характеризуется научная новизна и практическая ценность результатов диссертационной работы.

3 ПЕРВОМ РАЗДЕЛЕ рассмотрено состояние вопроса в области

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

С этих позиций проанализированы цепные, полевые и комбинированные методы моделирования. Показаны преимущества численных полевых методов, в том числе метода конечных элементов, перед другими методах®. Использование средств современной вычислительной техники'при реализации МКЭ позволяет эффективно решать многие задачи моделирования.■

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

На основе проведенного анализа выявлены и поставлены конкретные задачи по реализации метода конечных элементов в виде пакетов прикладных задач для математического моделирования.

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

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

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

Б

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

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

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

В ТРЕТЬЕМ РАЗДЕЛЕ рассмотрены вопросы реализации пакета программ расчета двумерного магнитного поля. Проведен анализ структур пакетов программ,■использующих метод конечных элементов. По функциональному назначению модули пакетов разделены на три части: модули подготовки данных для расчета, модули решения систем уравнений МКЭ и модули анализа результатов. Выбрана структура пакета, где модули выполняющие разные функции существуют в виде самостоятельных единиц, вызов которых производится с помощью программы -монитора (рис. 1). Связь между программами осуществляется через файлы данных с индивидуальными именами, размещаемыми в базе данных пакета. Работы по ведению базы осуществляются с помощью специальной программы, позволяющей оперативно просматривать содержимое ее разделов (быстрое воспроизведение контуров и сети на экране дисплея), удалять и перемещать их в архив.

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

о

АЦПУ ГРАФОПОСТРОИТЕЛЬ

Рис. 1

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

Рис. 2

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

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

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

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

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

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

Разработаны вопросы реализации процедур обработки результатов расчета поля МКЭ. Графические процедуры: вывод контура и сети области, картины поля, закраска по индексам сред и индукциям, кривых распределения индукции на экран, графопостооитель и принтер позволяют оценивать корректность отдельных этапов моделиро-вагля (генерацию сети, задание признаков материалов, граничных условий, нагрузок), документировать результаты расчета. Вычислительные процедуры позволяют определять локальные и интегральные параметры поля: потенциалы, индукции, потоки, потокосцепления обмоток, силы, моменты, НДС.

Программное обеспечение написано на языке F0RTRAN-77 в трех версиях: для стандартного комплекта персонального компьютера IBM PC XT/AT с одним дисплеем, принтером; для расширенного комплекта

Î2

компьютера IEM PC с дополнительным графическим дисплеем и графопостроителем (размерность решаемой задачи до 2 тысяч расчетных точек); для мини-ЭВМ VAX-11 (размерность задачи до 10000 расчетных точек). В первой верст графическая и цифровая информация выводится на один экран, разделенный на три части: графическое окно, окно меню и окно ввода-вывода численных значений. В других версиях пакета графическая информация выводится на отдельный графический дисплей. Функции работы с клавиатурой, экраном, мышью и принтером реализованы на языке С. Предусмотрено использование Графопостроителей CALCOMP, DMP, ДИГИГРАФ, HEWLETT-PACKARD.

В ЧЕТВЕРТОМ РАЗДЕЛЕ рассматриваются вопросы реализации пакета программ расчета трехмерного магнитного поля. В качестве расчетной функции используются скалярный магнитный потенциал при анализе потенциальных полей и векторный магнитный потенциал для вихревых полей. Число расчетных точек до 18 тысяч. Пакет состоит из трех частей: модуля подготовки трехмерной конечно-элементной модели, общего для скалярного и векторного потенциалов, программ определения значений потенциала в узлах сети (независимые программы для векторного и скалярного потенциалов) и модуля анализа результатов расчета. Режимы работы: интерактивный на этапах подготовки данных и анализа результатов и пакетный на этапе расчета. Расчет может производиться в активном процессе и как фоновая задача в многопользовательском режиме работы ЭВМ. Передача данных между модулями осуществляется через файлы на дисках. Выбранная структура пакета обеспечивает рациональное использование графи- . ческих устройств и времени центрального процессора мини-ЭВМ.

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

Метод заключается в следующем:

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

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

сколько слоев трехмерных конечных элементов;

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

На рис. 3 приведен пример двумерной сети и построенной на ее основе трехмерной сети половины полюсного деления синхронного микродвигателя С.аксонометрическая проекция).

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

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

Разработанное программное обеспечение написано на языке F0RTRAN-77+ в операционной среде DOS VMS-11. Использован графический дисплей MEGATEK и плоттер CALC0MP.

Рис. 3

В ПЯТОМ РАЗДЕЛЕ рассмотрены примеры применения разработанных пакетов программ для расчета двух- к трехмерных магнитных полей в электрических машинах. Из большого числа решенных задач отобраны наиболее характерные, показывающие возможности пакетов при моделировании широкого класса устройств с самой разнообразной конфигурацией: расчет угловой характеристики управляемого реактивного двигателя, анализ поля в воздушном зазоре линейного синхронного двигателя, расчет поля возбуждения синхронного микродвигателя с постоянными магнитами в двух- и трехмерной постановке, расчет трехмерного поля статической модели крлотурбогенерато-ра (сверхпроводящая седлообразная обмотка б цилиндрическом ферромагнитном экране), анализ распределения поля в воздушном зазоре синхронной электрической машины с электромагнитной редукцией частоты вращения.

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

Ка генерацию сети конечных элементов для приведенных примеров потребовалось от 5 до 20 минут в интерактивном режиме работа. Время счета одной итерации по нелинейности для задачи размером около 2000 узлов на ПЭЕМ с тактовой частотой 12 МГц - 8 минут, на мини-ЭВМ 0.6 минуты чистого процессорного времени. Результаты расчета сравниваются с экспериментальными. Максимальная погрешность при определении индукции составляет б процентов (трехмерные задачи), при вычислении моментов - э процентов.

3 ПРИЛОЖИ' приведет« структуры данных описания сети конечных элементов и данных описания составных блоков, принятые в разработанных пакетах и алгоритм генерации сети составных блоков, таблицы коэффициентов интерполяционных полиномов кривых намагничения часто используемых ферромагнитных материалов, документы, подтверждающие внедрение и практическое использование результатов работы.

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

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

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

Tl.. Q _ 7

образок:

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

2. Разработай; пакеты прикладных программ, обеспечивающие генерацию конечно-элементной модели устройства, задание свойств материалов, граничных* условий я источников возбуждения поля, формирование и реаюние' систем уравнений, интерпретацию получаемых результатов в' Еиде картин поля, значений тех или иных параметров в задаваемых точках, графиках изменения параметров и их гармонического анализа. В качестве источников возбуждения поля могут использоваться объемные, линейные и точечные токи, постоянные магниты.

3. Разработаны'алгоритмы и процедуры интерактивного построения конечноэлементных моделей двухмерных и трехмерных объектов с визуализацией результатов на дисплее, графопостроителе, мозаичном принтере. Разработаны средства анализа поля как в локалышх областях так и во ьсей машине. .Модули подготовки данных и анализа результатов пакетов без существенных изменений могут быть использованы для других задач, описываемых уравнениями Лапласа и Пуассона. ■ ~

Л. Проведен сравнительный анализ различных методов решения систем линейных уравнений МКЭ при использовании ПЭВМ с ограниченны;,: объемом оперативкой памяти (исключения по Гауссу, тройной факторизации, Холецкого, сопряженных градиентов). Разработаны программы, реализующие ети метода, проведено испытание программ на типоеых задачах расчета поля в электрических машинах. Выбран метод тройкой факторизации, обеспечивающий наибольшую скорость решения при высокой точности и приемлемыми требованиям:! к. оперативной памяти.

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

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

гэ

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

Основные положения диссертации опубликованы в следующих работах:

1. Распределение электромагнитного поля в зубцовом делении ротора / Горелик Л.В., Ларьки А.Г., Дарьин С.Г., Попов П.Г. // Исследование электрических машин: Тр. ВНИПТИЭМ .- Владимир, 1981. С. 185-187."

2. Дарыш А.Г'., Дарьин С.Г. Применение метода проводимостей зубцовых контуров для расчета электромагнитных полей в электрических машинах,- М.: Информэлектро, 1985 .- 32 с.

3. Дарыш С.Г. Пакет программ расчета электромагнитного поля электрических машин методом конечных элементов // Разработка и освоение производства асинхронных двигателей серии АИ: Тр. ВНИПТИЭМ. Владимир, 1686 С. 31-39.

4. Архипов A.B., Дарьин А.Г., Дарьин С.Г. Диалоговый рехнм при аппроксимации характеристик намагничения ферромагнитных материалов // Исследование параметров и проектирование электрических машин: Тр. ВНИПТИЭМ Владимир, 1987 .-С. 32-40.

!). Дарьин С.Г., Максимов E.H. ' Расчет"трехмерных магнитных полей в электрических машинах методом конечных элементов // Автоматизация проектирования и производства асинхронных двигателей единых серий: Тр. ВНИПТИЭМ .- Владимир, 1988 .- С. 38-47.

6. Дарьин А.Г., Дарьин С.Г. Расчет электромагнитных полей методом конечных элементов на персональном компьютере // Состояние и перспехтивы совершенствования разработки, производства и применения низковольтных электродвигателей переменного тока: Тез. дохл. / VIII Всесоюэ. науч.-техн. кокф. 14-18 ноября 1988 г. Владимир, 1988 С. 78.

7. Дарыш А.Г., Дарьин С.Г., Стрельбицкий Э.К. Пакет программ расчета магнитных полей методом конечных элементов // Системы автоматизированного проектирования в электротехнической промышленности: Тр. ВНИИЭМ Н., 1989 Т. SO .- С. 18-24.

8. Горелик Л.В., Дарьин С.Г., Талклинскнй А.Р. К вопросу об энергетических показателях асинхронных двигателей с сегментированным магн'итопроводом // Электротехника .- 1989 te 3 .- С. 10-1".