авторефераты диссертаций БЕСПЛАТНАЯ  БИБЛИОТЕКА

АВТОРЕФЕРАТЫ КАНДИДАТСКИХ, ДОКТОРСКИХ ДИССЕРТАЦИЙ

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

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

На правах рукописи

БУЛАТОВ Леонид Николаевич

АВТОМАТИЗАЦИЯ ПРОЕКТИРОВАНИЯ

АСИНХРОННЫХ МАШИН С ИСПОЛЬЗОВАНИЕМ

ПОЛЕВЫХ ДИНАМИЧЕСКИХ МОДЕЛЕЙ

Специальность:

05.13.12. – Системы автоматизации проектирования

(электротехника, энергетика)

Автореферат

диссертации на соискание учной степени

кандидата технических наук

Иваново – 2013 1

Работа выполнена в федеральном государственном бюджетном образовательном учреждении высшего профессионального образования «Ивановский государственный энергетический университет имени В.И. Ленина» (ИГЭУ).

Тихонов Андрей Ильич, доктор

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

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

Кроль Татьяна Яковлевна, кандидат технических наук, Ивановский центр информационных технологий – филиал ОАО «Электроцентромонтаж», директор Открытое акционерное общество «Яро

Ведущая организация славский электромашиностроительный завод» (ОАО «ELDIN»)

Защита состоится 25 июня 2013 г. в 11 часов на заседании диссерта ционного совета Д 212.064.02 при ИГЭУ по адресу: г. Иваново, ул. Раб факовская, 34, ауд. Б-301.

Отзывы на автореферат в двух экземплярах, заверенные печатью организации, просим присылать по адресу: 153003, г. Иваново, Рабфаковская, 34, Ученый совет ИГЭУ. Тел.: (4932) 38-57-12, факс: (4932) 38-57-01. Е-mail: uch_sovet@ispu.ru.

С диссертацией можно ознакомиться в библиотеке ИГЭУ, автореферат размещн на сайте www.ispu.ru.

Автореферат разослан 24 мая 2013 г.

Учный секретарь диссертационного совета Тютиков Д 212.064.02 Владимир Валентинович

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

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

До 80% всех электродвигателей, выпускаемых промышленностью, приходится на долю асинхронных двигателей (АД), поэтому разработка эффективных методов моделирования в целях совершенствования конст рукции этих машин является актуальной задачей. Одним из направлений совершенствования конструкции АД является улучшение их динамиче ских характеристик.

В настоящее время наблюдается тенденция вытеснения эксперимен тальных методов исследования опытных образцов электрических машин численным экспериментом с использованием точных математических моделей на основе теории поля и цепей. Наиболее распространенными в этом плане являются модели, построенные на основе метода конечных элементов. Преимуществом таких моделей является возможность учета особенностей конструкции исследуемого устройства, в том числе, в ди намических режимах работы. Универсальные системы инженерных рас четов, позволяющие формировать и исследовать полевые модели элек трических машин, являются необходимыми элементами современных систем автоматизированного проектирования (САПР) АД. Среди них популярны такие программные продукты, как Ansys, Maxwell, Comsol, ElCut и др. Несмотря на широкие возможности, такие продукты имеют ряд недостатков, в частности высокую стоимость, наличие излишнего с точки зрения конечного пользователя функционала, в то же время огра ничение возможностей систем изначально заложенным в них функцио налом, закрытость большинства алгоритмов, значительное время реше ния полевых задач и т.п. Это ограничивает возможности научного ис пользования таких пакетов.

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

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

Диссертационная работа выполнялась в Ивановском государственном энергетическом университете. Работа соответствует паспорту специаль ности 05.13.12 «Системы автоматизации проектирования (по отраслям)», так как содержит результаты разработки научных основ построения средств САПР, в частности разработки и исследования моделей, алго ритмов и методов для синтеза и анализа проектных решений в САПР АД.

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



Задачи, решаемые в диссертации:

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

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

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

4. Адаптация разработанных программных средств к использованию в современных высокопроизводительных кластерных системах.

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

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

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

Научная новизна.

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

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

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

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

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

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

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





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

- базовый вариант расчетной подсистемы САПР АД, основанный на модели динамических режимов АД, который может быть использо ван как для расчета новых асинхронных машин, так и для проведения научных исследований, в том числе и в рамках НИОКР.

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

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

По результатам работы получено два акта внедрения: в ОАО "ELDIN" (г. Ярославль) и ОАО НИПТИЭМ (г. Владимир).

Апробация работы. Результаты работы докладывались на конферен циях: на международной научно-технической конференции «Современ ное состояние, проблемы и перспективы энергетики и технологии в энер гостроении» (Бенардосовские чтения, г. Иваново, ИГЭУ, 2009, 2011 гг.), на международной научно-технической конференции студентов и аспи рантов «Радиоэлектроника, электротехника и энергетика» (г. Москва, МЭИ, 2009, 2011, 2013 гг.), на региональной научно-технической кон ференции студентов и аспирантов «Энергия» (г. Иваново, ИГЭУ, 2008, 2009, 2012 гг.).

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

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

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

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

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

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

При этом время на проектирование каждой партии оказывается ограни ченным. Для того чтобы выдержать конкуренцию в этих условиях пред приятия вынуждены осваивать современные средства проектирования и моделирования (CAD (Computer-aided design)/CAE (Computer-aided engineering)), позволяющие отказаться от натурных исследований опыт ных образцов новой продукции в пользу численного моделирования.

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

Альтернативой современным CAD/CAE-системам, пытающимся вме стить в себя весь опыт численного решения задач математической физи ки, может служить модульный подход к построению САПР, при котором САПР организуется на базе некоторой среды, позволяющей создавать, отлаживать и запускать на выполнение подпрограммы, реализующие проектные операции. Задача, поставленная в данной работе, включает в себя разработку модели асинхронного двигателя, которая может быть использована в САПР АД на стадии НИОКР и поверочного расчета, и создание на ее основе расчетных подсистем САПР АД.

Используемые в САПР математические модели можно условно разде лить на три уровня: микроуровень (уровень полевых моделей), макро уровень (уровень цепных моделей) и метауровень (уровень структурных схем). В работе приводится обзор цепных моделей, применяемых для расчета динамических режимов АД.

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

Для расчета магнитного поля при создании комбинированной модели динамических режимов АД использовалась модифицированная версия разработанной в ИГЭУ библиотеки конечно-элементного моделирования магнитных полей EMLib, в которой для расчета магнитного поля реша ется уравнение Пуассона. Расчет нелинейной системы уравнений, к ко торой сводится решение задачи, осуществляется итерационным методом.

Значения магнитных потенциалов в узлах модели на каждой итерации рассчитываются по формуле A t+1 = A t - t. Dt, (1) где At+1, At – векторы магнитных потенциалов в узлах конечно элементной сетки на следующей и текущей итерациях соответственно;

t – коэффициент релаксации;

D – вектор приращений магнитных потен циалов.

Вектор D вычисляется по формуле J t D t = - F t. (2) где J – матрица Якоби;

F – вектор невязок.

Использование собственного средства расчета поля позволяет разра батывать модели и подсистемы САПР, не зависящие от крупных про граммных продуктов.

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

Для сокращения времени расчета динамических режимов АД с ис пользованием комбинированной модели может быть применена техноло гия параллельных вычислений на кластерных системах. Появление кла стерных систем относят к 90-м годам XX века. В общем случае основ ным параметром классификации параллельных систем является наличие общей (SMP) или распределенной (MPP) памяти. В данной работе при менялась технология MPP.

Вторая глава посвящена разработке платформонезависимой библио теки конечно-элементного моделирования магнитного поля.

Ядром подсистемы расчета динамических режимов работы АД, а также наиболее сложным ее компонентом с точки зрения разработки яв ляется библиотека конечно-элементного моделирования магнитных по лей. Такая библиотека была разработана в ИГЭУ. Однако она недоста точно адаптирована к использованию в моделях динамических режимов электрических машин, а также к расчетам на кластерной технике. Поэто му было принято решение о переносе алгоритмов библиотеки EMLib на кроссплатформенную основу с одновременным устранением указанных недостатков. В качестве системы программирования использовалась сво бодно распространяемая система CodeBloks, позволяющая создавать приложения на языке С++. С ее помощью была создана автономная биб лиотека MKE_cpp, имеющая следующие отличия от библиотеки EMLib:

1) MKE_cpp написана на базовой версии языка С++, что расширяет возможности ее использования;

2) В MKE_cpp оптимизированы вопросы выделения и освобождения памяти;

3) MKE_cpp более компактна.

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

В основе модели динамических режимов работы АД лежит поиск за висимости приращений потокосцеплений обмоток от приращений токов в обмотках и изменения пространственной конфигурации модели. В мо дели АД с короткозамкнутым ротором приращение потокосцепления k го контура с током k определяется по формуле di k d n (3) d k k j.

j1 i j dt dt dt k k k k Для определения необходимо рассчитать, при, i j i j ращения потокосцеплений k при изменении токов в контурах ij и угла поворота ротора на величину ij и. Для этого используется модель магнитного поля исследуемого двигателя. При расчете частных произ водных от потокосцеплений контуров с током принимается, что p const. С учетом этого члены матрицы J и вектора F в уравнении (2) могут быть определены по следующим формулам:

N 1 j bpibpj cpicpj Api 3 jpSp, Fi (4) p 1 4Sp p j1 i1N 1 1 p Ni bpi bpj cpi cpj Jij. (5) 4Sp p i1N p j1 N где N – количество элементов конечно-элементной сетки;

b, c – коэф фициенты, зависящие от координат узлов конечно-элементной сетки;

Sp, jp, p, Api – соответственно площадь, плотность тока, магнитная прони цаемость и векторные магнитные потенциалы в узлах p-го элемента.

При условии p const (2) является системой линейных уравнений и окончательные значения приращений векторного магнитного потенциала вычисляются в ходе первой итерации, против 10 и более итераций при стандартном конечно-элементном расчете. Дополнительное ускорение расчета достигается благодаря отсутствию необходимости обращения к библиотеке магнитных характеристик материалов.

Точность расчета момента через магнитную энергию/коэнергию при анализе динамических режимов оказалась неудовлетворительной. По этому было принято решение рассчитывать электромагнитный момент по результатам конечно-элементного полевого расчета с использованием формулы суммарного момента магнитных сил, действующих на тела, заключенные в заданном объеме. Для плоскопараллельной задачи фор мула для расчета момента имеет вид 1 NE X H y Y H x X Bx Y By L L, M (6) 2 R 1 X By Y Bx X Bx Y By где X,Y – координаты центра отрезка L, соединяющего точки пересече ния элемента конечно-элементной модели выбранной токоограничи вающей поверхностью;

NE – число элементов, пересекаемых токоограни чивающей поверхностью;

Hx, Hy, Bx, By – компоненты векторов напря женности и индукции поля в точке X,Y, которые могут быть получены из конечно-элементной модели.

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

Модель строится на основе системы дифференциальных уравнений d i u Ri.

L (9) dt Механическая составляющая представлена уравнением d M. (10) dt J Здесь L – матрица индуктивностей контуров;

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

R – матрица сопротивлений контуров;

i – вектор мгновенных значений токов в контурах;

– угловая скорость вращения ротора;

J – момент инерции ротора;

M – сумма моментов на валу ротора.

При формировании системы уравнений в модели используются схемы замещения, показанные на рис. 1–3.

u1 es rs Ls1 В приведенной на рис. 3 схеме за мещения паз ротора разбит на ns сло is1 ев. Это позволяет учесть эффект вы u2 es rs Ls 2 теснения тока в пазах короткозамкну того ротора. При этом количество уравнений, описывающих ротор, уве u3 is 2 личивается в ns раз. Соответственно, es увеличивается и число обращений к конечно-элементной модели на каж Ls 3 rs 3 дом шаге итераций.

Рис. 1. Схема замещения статора LRZ 1 rRZ 1 LRZ rRZ rR1 LR 2 rR LR rRZ rR 3 rRZ rR1 rR iRZ iR 2 i RZ i R i RZ eR 3 eRZ eR1 eR 2 eRZ L'RZ 1 r 'RZ 1 L'RZ r 'RZ L'R1 r 'R1 L'R 2 r 'R Рис. 2. Схема замещения короткозамкнутого ротора LR1 rR r rR rR1n rR11 rR12 rR 2 n R iR11 iR12 iR1n eR eR eR1n L'R1 r 'R eR11 eR12 eR 2 n Рис. 3. Схема замещения ротора с разбиением пазов на слои Величина ns зависит от формы пазов и размеров магнитной системы машины.

На основании приведенных схем замещения формируется система дифференциальных уравнений:

d 1 d di1 di u1 u 3 i1rs1 dt (i1 i 2 )rs3 dt (L s1 L s3 ) dt L s3 dt, u u i r d 2 (i i )r d 3 (L L ) di 2 L di1, 2 3 2 2 dt s2 s3 s 1 2 s dt dt dt d i 0 i di LRi L'Ri Ri i r i r N nsi N nsi i Ri2 2ns R 2ns r, dt i1..ns dt di Ri d i (9) 0 iins rRins2 ii rRi2 iins rRins2 L Ri L'Ri, dt i3 ns...N ns dt d i 0 iins rRins2 ii rRi2 i 2ns rR 2ns LRi L'Ri Ri di, dt i N ns1...N dt d, dt d M M M н, dt J rR1i rR 2i rR1i r 'R1i, если j 1, rRji rRji rR ( j1)i rRji r 'Rji, если 1 j Z2, (10) rRZ2i rR1i rRZ2i r 'RZ2i, если j Z2.

В матричном виде данная система уравнений имеет вид d L x U, (11) (14) dt где x – вектор искомых величин.

Шаг по времени в модели должен определяться из условия p p dt max (12), 2f Z2 4f Z где Z – число зубцов ротора в случае, когда оно равно или превышает число зубцов статора, или число зубцов статора в другом случае;

p – число пар полюсов АД;

f – частота.

Основные трудности при расчете динамических режимов АД с ис пользованием разработанной модели связаны с расчетом матрицы L.

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

m m В модели с прямым обращением к полевой модели и, соот, i ветственно, элементы матрицы L определяются с использованием серии конечно-элементных расчетов магнитного поля на каждом шаге итера ции. Главными достоинствами данной модели являются точность и уни версальность, а основным недостатком – низкая скорость расчета.

В модели с предварительным расчетом матрицы потокосцеплений элементы матрицы L определяются из полученных предварительно зави симостей i = f(Fx, Fy, ), где Fx и Fy - составляющие результирующей намагничивающей силы машины. В этой модели приняты следующие допущения:

1) обмотки АД создают синусоидальную намагничивающую силу;

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

При расчете переходного процесса данным методом приращения по токосцеплений по времени определяются по формуле d m m dFx m dFy m d 3 Z2 ns m dFxj m dFyj m d Fy dt dt dt Fx dt Fy dt j 1 Fx dt dt (13) 3 Z 2 ns di j m d 3 Z2 ns di j m d m m K xj Lmj K yj, dt dt dt j 1 Fx Fy dt j Приращения потокосцеплений определяются из полученных на этапе предварительного расчета зависимостей i = f(Fx, Fy, ), аппроксимиро ванных трехмерными сплайнами.

Второе допущение серьезно снизило точность расчета, поэтому был разработан модифицированный вариант метода, в котором на этапе предварительного расчета определяются зависимости i = f(FxS, FyS, FxR, FyR, ), где FxS, FyS, FxR, FyR – составляющие НС статора и ротора со ответственно по осям X и Y. Это повышает время предварительного рас чета и объем хранимых данных, однако позволяет получить точность, близкую к точности модели с прямым обращением к полевой модели.

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

Было произведено сравнение результатов, полученных с использова нием моделей, с опытными данными (рис. 4), а также с результатами, полученными с использованием современных систем. Сравнение показа ло, что созданные модели позволяют получать достоверные результаты, не уступая по скорости расчета системам, имеющимся на рынке.

б) а) Рис. 4. Зависимость скорости вращения ротора от времени: a – модель с непосредст венным обращением к полевой модели;

б – модель с предварительным расчетом мат риц потокосцеплений;

1 – расчетные данные;

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

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

Четвёртая глава посвящена разработке структуры и подсистем САПР АД. Основные подсистемы САПР АД показаны на рис. 5.

Подсистема проектного расчета Оптимизационная Подсистема поверочного расчета подсистема Электромагнитный расчет Механический расчет Тепловентиляционный расчет Подсистема генерации технической документации Рис. 5. Подсистемы САПР АД Созданный в рамках данной работы математический аппарат был ис пользован для создания подсистемы поверочного расчета САПР. В на стоящее время большое распространение получают программные про дукты, обладающие широкими интерфейсными возможностями, а также встроенной средой программирования. К данным продуктам относится программа Microsoft Excel, имеющая встроенный язык программирова ния Visual Basic. Она была выбрана в качестве основы для создаваемых средств САПР.

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

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

Структура разработанной подсистемы поверочного расчета САПР по казана на рис. 8. Базовый модуль подсистемы, созданный на основе Microsoft Excel 2007, координирует взаимодействие всех остальных мо дулей. В этом модуле производится подготовка динамической модели к расчету. Параметрический генератор так же выполнен на базе Microsoft Excel 2007, и для удобства совмещен с ним в одну книгу Excel.

Рис 6. Входные параметры параметрического генератора Рис. 7. Результат работы параметрического генератора Модуль динамического расчета, реализованный в виде отдельной библиотеки AD_SLVR.dll, выполняет расчет двигателя с использованием комбинированной динамической модели. Модуль обработки данных служит для сбора данных расчета и представления их в удобном для пользователя виде.

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

Этот модуль обеспечивает базу для создания новой подсистемы САПР – подсистемы генерации быстродействующей модели АД.

Подсистема проектного расчета АД Данные проектного расчета Подсистема поверочного расчета САПР АД Данные проектного Данные для создания расчета, MKE-модели информация о модели Базовый модуль Данные Параметрический динамической Модуль обработки MKE-модель генератор модели данных Модуль динамического Модуль сбора данных Хранилище расчета АД для ускоренной модели данных Графический Данные расчета процессор Рис. 8. Подсистема поверочного расчета САПР АД Целью разработки этого модуля было создание механизма формиро вания быстродействующей динамической модели АД, обладающей дос тоинствами численных моделей на основе метода конечных элементов.

Для этого была использована модель с предварительным расчетом мат риц потокосцеплений. Зависимости i = f(FxS, FyS, FxR, FyR, ) для каж дой обмотки собираются с использованием подсистемы поверочного расчета.

Массив данных, из которых фор Рабочая зона Fy мируются сплайны, может быть очень велик. Для сокращения объема данных, используемых в каждый мо F мент времени, могут быть использо ваны сплайны, построенные на осно ве малой части полученных в ходе Fx предварительного расчета данных (рис. 9). При этом модель должна ав томатически смещать границы рабо чей зоны в случае приближения к ним рабочей точки. Это обеспечивает высокое быстродействие модели без Рис. 9. Сокращение объема используе потери точности.

мых данных.

ОСНОВНЫЕ РЕЗУЛЬТАТЫ РАБОТЫ 1. Математический аппарат полевой динамической модели асин хронной машины, основанный на использовании результатов конечно элементного расчета магнитного поля при расчете разветвленной элек трической схемы замещения машины.

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

3. Модифицированная платформонезависимая библиотека конечно элементного моделирования магнитного поля.

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

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

6. Расчетная подсистема системы автоматизированного проектирова ния асинхронных машин, основанная на разработанном методе модели рования динамических режимов.

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

Публикации по теме диссертации Научные статьи, опубликованные в изданиях по списку ВАК 1. Булатов, Л. Н. Разработка динамической модели асинхронной маши ны с использованием результатов конечно-элементного расчета / Л. Н.

Булатов, А. И. Тихонов // Вестник ИГЭУ. – 2012. – Вып. 4. – С. 32-34.

2. Булатов, Л. Н. Динамическая модель асинхронной машины с прямым обращением к конечно-элементной модели магнитного поля / Л. Н.

Булатов, А. И. Тихонов // Вестник ИГЭУ. – 2013. – Вып. 1. – С. 40-44.

Публикации в других изданиях 1. Булатов, Л. Н. Разработка кластерной полевой динамической модели асинхронного двигателя / Л. Н. Булатов, А. И. Тихонов, Д. В. Руб цов // Энергия 2008: тез. докл. региональной науч.-техн. конф. студен тов и аспирантов / Иван. гос. энерг. ун-т. – Иваново, 2008. – С. 24–25.

2. Булатов, Л. Н. Разработка независимой библиотеки конечно элементного моделирования / Л. Н. Булатов, А. И. Тихонов // Энергия 2009: тез. докл. региональной науч.-техн. конф. студентов и аспиран тов / Иван. гос. энерг. ун-т. – Иваново, 2009. – Т. 2. – С. 25–26.

3. Тихонов, А. И. Моделирование эффекта вытеснения тока в пазах ро тора асинхронного двигателя при пуске с использованием многопро цессорной техники / А. И. Тихонов, Л. Н. Булатов // Тезисы докл. ме ждунар. науч.-техн. конф. (XV Бенардосовские чтения) / Иван. гос.

энерг. ун-т. – Иваново, 2009. – С. 24–25.

4. Булатов, Л. Н. Разработка кластерной полевой модели асинхронного двигателя / Л. Н. Булатов, А. И. Тихонов // Радиоэлектроника, элек тротехника и энергетика: тез. докл. XV междунар. науч.-техн. конф.

студентов и аспирантов. – М.: МЭИ, 2009. – Т. 2. – С. 10–11.

5. Булатов, Л. Н. Динамическая модель асинхронного двигателя с корот козамкнутым ротором / Л. Н. Булатов, А. И. Тихонов // Вестник науч но-промышленного общества. – М.: Алев-В, 2011. – Вып. 1. – С. 35–40.

6. Булатов, Л. Н. Комбинированная динамическая модель асинхронного двигателя / Л. Н. Булатов, А. И. Тихонов // Радиоэлектроника, элек тротехника и энергетика: тез. докл. XVII междунар. науч.-техн. конф.

студентов и аспирантов. – М.: МЭИ, 2011. – Т. 2. – С. 7–8.

7. Булатов, Л. Н. Варианты модели расчета переходных процессов в асинхронной машине с обращением к расчету поля / Л. Н. Булатов // Тезисы докл. междунар. науч.-техн. конф. (XVI Бенардосовские чте ния) / Иван. гос. энерг. ун-т. – Иваново, 2011. – Т. 3 – С. 148–151.

8. Булатов, Л. Н. Динамическая модель асинхронного двигателя с ко роткозамкнутым ротором / Л. Н. Булатов, А. И. Тихонов // Инноваци онные проекты молодых ученых за 2010 г.: cб. отчетов / Иван. гос.

энерг. ун-т. – Иваново, 2011. – C. 135–140.

9. Булатов, Л. Н. Динамическое моделирование асинхронного двигателя машин с использованием конечно-элементного расчета / Л. Н. Була тов // Инновационные проекты молодых ученых за 2011 г.: cб. отчетов / Иван. гос. энерг. ун-т. – Иваново, 2012. – Т.2 – C. 66–73.

10. Булатов, Л. Н. Оптимизация комбинированной динамической модели асинхронной машины / Л. Н. Булатов, А. И. Тихонов // Энергия 2012:

тез. докл. региональной науч.-техн. конф. студентов и аспирантов / Иван. гос. энерг. ун-т. – Иваново, 2012. – Т. 4. – С. 66–71.

11. Булатов, Л. Н. Расчет пуска асинхронного двигателя с прямым обра щением к конечно-элементной модели магнитного поля / Л. Н. Була тов, А. И. Тихонов // Радиоэлектроника, электротехника и энергетика:

тез. докл. XIX междунар. науч.-техн. конф. студентов и аспирантов. – М.: МЭИ, 2013. – Т. 2. – С. 141.

Свидетельства на программные продукты 1. Свидетельство о государственной регистрации программы для ЭВМ.

Платформонезависимая библиотека конечно-элементного моделирова ния магнитного поля / А. И. Тихонов, Л. Н. Булатов. – № 2011614852;

22.06.2011. – М.: Федеральная служба по интеллектуальной собствен ности, патентам и товарным знакам.

2. Свидетельство о государственной регистрации программы для ЭВМ.

Программа расчета динамических режимов асинхронной машины с уче том эффекта вытеснения тока в пазах ротора с непосредственным обра щением к библиотеке конечно-элементного расчета поля / А. И. Тихонов, Л. Н. Булатов. – № 2011615132;

30.06.2011. – М.: Федеральная служба по интеллектуальной собственности, патентам и товарным знакам.

БУЛАТОВ Леонид Николаевич АВТОМАТИЗАЦИЯ ПРОЕКТИРОВАНИЯ АСИНХРОННЫХ МАШИН С ИСПОЛЬЗОВАНИЕМ ПОЛЕВЫХ ДИНАМИЧЕСКИХ МОДЕЛЕЙ

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата технических наук Подписано в печать 20.05.2013г. Формат 60х841/ Печать плоская. Усл.печ.л.1.16.

Тираж 100 экз. Заказ № 112.

ФГБОУ ВПО «Ивановский государственный энергетический университет имени В.И.Ленина»

Отпечатано в УИУНЛ ИГЭУ 153003, г. Иваново, ул. Рабфаковская,

 

Похожие работы:





 
2013 www.netess.ru - «Бесплатная библиотека авторефератов кандидатских и докторских диссертаций»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.