Structural analysis in ANSYS of the axial flux permanent magnet generator

Full Text

Низкоскоростные торцевые синхронные генерато- и имеют высокотехнологичную конструкцию и низры находят применение для нужд малой энергетики. кую материалоемкость [1]. В торцевых генераторах Ввиду своих малых осевых размеров и возможности магнитный поток направлен параллельно оси вращеработы на низких оборотах их используют для преоб- ния (аксиально), в то время как в цилиндрических разования энергии движения ветро- или гидротурби- перпендикулярно оси вращения (радиально). ны в электрическую энергию. Торцевые генераторы Существует несколько различных вариантов кон-отличаются от традиционных цилиндрических син- струкций торцевых генераторов, неизменным остается хронных генераторов направлением магнитного потока то, что во всех топологиях аксиально намагниченные 63 Вестник СибГАУ. 2014. № 3(55) магниты вмонтированы в круглый диск ротора. Для варианта конструкции с одним ротором и двумя статорами магниты располагаются с двух сторон диска ротора, который размещен между статорами. В такой компоновке на ротор действуют две силы, одинаковые по величине и противоположные по направлению, эти силы компенсируются в случае, если ротор помещен строго посередине между статорами. Оба диска статора также испытывают на себе действие сил. В компоновке с двумя роторами и одним статором два зеркально расположенных ротора с постоянными магнитами, вмонтированными с одной стороны дисков ротора, установлены на вал. Статор располагается между двумя дисками ротора и испытывает на себе силы, с которыми роторы притягиваются к статору. Если статор расположен не строго посередине между роторами, то действует результирующая сила. В многоступенчатых топологиях количество роторов на один больше количества статоров. При этом в конструкции используются два типа роторов: внешний (постоянные магниты вмонтированы только с одной стороны диска ротора) и промежуточный (постоянные магниты вмонтированы с двух сторон диска ротора). В процессе сборки генератора будет действовать максимальная притягивающая сила между статором и ротором. Действие этой силы необходимо учитывать при проведении структурного анализа электрической машины. Для выбора оптимальной конструкции необходимо провести ряд проверочных расчетов для нескольких компоновок машин. Существуют аналитические зависимости, связывающие максимально допустимые деформации с конструктивными размерами. Однако эти выражения не дают возможности быстро и с большой точностью получить результаты. Поэтому чаще используются численные методы решения, например, метод конечных элементов [2]. Программный продукт ANSYS позволяет с достаточным уровнем достоверности оценить влияние притягивающей силы на конструкцию электрической машины, в частности, активную ее часть - статор и ротор. Дальнейший расчет будет проводиться в среде ANSYS. Среда конечно-элементного анализа ANSYS позволяет получить на выходе из решения файл с расширением «log», хранящий в себе все данные расчетной модели с возможностью их корректировки, написанный на языке APDL. Это позволяет разработать программный продукт для управления параметрами решения. Такой подход позволяет сделать построение расчетной модели автоматизированным, что существенно сокращает затрачиваемое время и не требует высокой квалификации пользователя при работе в среде ANSYS. В некоторых случаях работа в среде ANSYS требует не только достаточных навыков, но и большой внимательности, что затрудняет получение решения в короткие сроки, тогда когда время подчас является критическим фактором. Следует отметить, что программный продукт является частью автоматизированной среды проектирования торцевого генератора [3; 4]. Интерфейс разработанного программного продукта с введенными данными показан на рис. 1. В программу задаются требуемые геометрические размеры ротора, статора, вала, типы элементов в среде ANSYS, параметры материала, величина силы магнитного притяжения, характеристики подшипников. Также возможна корректировка размерности сетки. Все эти переменные являются параметрами, которые можно изменять и в автоматизированном режиме перестраивать расчетную модель. А регулировка плотности сетки позволяет получить результаты с требуемой точностью. Программный продукт может быть использован в том числе как инструментарий для проведения исследований влияния размерности сетки на точность получаемых результатов при проведении расчетов в различных дисциплинах, таких как статика, динамика и анализ электромагнитных процессов. В программе используется буквенный код обозначения номера объема, который позволяет в выходном файле наглядно отслеживать нумерацию создаваемых объемов, а также при дальнейших изменениях конструкции позволяет вставлять новые построения, избегая при этом полной корректировки программного кода. Формируемый программой log-файл содержит в себе процедуру решения задачи в среде ANSYS, которая по шагам описывается ниже. Программа препроцессинга геометрии в ANSVS Параметры [ *.l&q файхі.І Размеры вала, м Размеры ступииьі, м Средняя сига тгокения, Н Тип элемента 2D R1 L1 R1 L1 F 1 PLANE 182 0.019 0.055 0.048 0,014 700 R2 L2 R2 L2 [sOLID 185 -*-1 0.021 R3 0.064 L3 0.050 R3 0.030 L3 Подшипники Кг - радиальная жесткость Н/м Параметры материала 0.023 0.024 0.080 0.020 100000 Модуль упругости Юнга (Н/м2) R4 L4 R4 L4 Kz - осевая жесткость Н/м 210000000000 0.026 0.026 0.048 0.006 34000 Коэффициент Пуассона 0.28 Плотность (кг/мМ R5 L5 Cv - кяэф рассеиван. энергии 0.038 RS 0.006 L6 Размеры писка ротора, Р.1 Размеры магнитов,м R1 0.2 0.033 0.096 0.110 0.056 R7 L7 R2 L1 0.023 0.024 0.198 0.004 Расчет 1 L1 Количество в четверти 0.018 б Рис. 1. Интерфейс программного продукта 64 Математика, механика, информатика Выполнение конечно-элементного анализа в среде ANSYS подразумевает выполнение следующих этапов [5]: 1) определение типов элементов; 2) задание свойств материалов; 3) создание расчетной параметризованной модели; 4) разбиение на конечно-элементную сетку; 5) приложение нагрузок и ограничений; 6) выполнение расчета; 7) просмотр результатов. Для решения задачи использованы следующие элементы: плоский PLANE182 и объемный SOLID185. Элемент PLANE182 предназначен для двухмерного моделирования плоско-деформируемых или осесимметричных задач. Элемент состоит из четырех узлов, имеющих 2 степени свободы - перемещение вдоль оси X и Y. В данной задаче плоский элемент используется для разбиения на сетку поверхностей модели с целью последующего экструдирования, что позволяет полностью контролировать процесс создания регулярной объемной сетки с заданными характеристиками. Элемент SOLID 185 предназначен для трехмерного моделирования. Элемент состоит из восьми узлов, имеющих 3 степени свободы - перемещение вдоль осей X, Y и Z [6]. В задаче используется модель изотропного материала со следующими характеристиками: Модуль Юнга (Young’s Modulus) = 2,1е 11 (Па); коэффициент Пуассона (Poisson’s ratio (minor)) = 0,28; плотность (Mass density) = 7700 (кг/м3) [7]. Для построения расчетной модели конструкция генератора была несколько упрощена. Поскольку искомой является величина деформации диска ротора, статоры заменены действием сил магнитного притяжения по обеим сторонам диска ротора. Неизменными в расчетной модели (рис. 2) остаются вал, ступица, диск ротора и магниты. Рис. 2. Расчетная модель торцевого генератора с одним ротором и двумя статорами Расчетная модель имеет две плоскости симметрии и ось вращения. Для расчета методом конечных элементов необходимо создание сеточной модели, т. е. замена твердотельной модели набором конечных элементов. В основу разбиения сетки положено разбиение торца вала. С учетом полученных элементов генерируется трехмерная сетка для всей четверти модели (рис. 3). Параметризация геометрии позволяет организовать процесс поиска нужного узла или другого объекта, в том числе в зависимости от координат, и реализовать оптимизацию по заданным критериальным ограничениям и параметрам. Такой принцип построения позволяет задавать нагрузку в каждом узле и выполнять требования к качеству сетки. Сетка проверяется на отсутствие вырожденных углов, на точность стыковки узлов сетки в сегментах модели и др. Так как диаметр вала в несколько раз меньше диаметра диска ротора, появилась необходимость в применении шаблона соединения частей сетки, без которого размеры элементов увеличивались бы с увеличением диаметра, что привело бы к огрублению результатов. Рис. 3. Конечно-элементная сетка Путем варьирования количества ступеней генератора создается расчетная модель генератора большей мощности (рис. 4). б Рис. 4. Расчетная модель двух- (а) и трехступенчатых (б) генераторов а 65 Вестник СибГАУ. 2014. N 3(55) Помимо количества ступеней существует возможность создать расчетную модель активной части торцевого генератора, варьируя следующими параметрами: 1) количество магнитов; 2) одностороннее или двустороннее приложение сил в зависимости от наличия в конструкции статора с одной стороны или с двух сторон от диска ротора; 3) геометрические параметры и др. Сила магнитного притяжения на постоянных магнитах для анализа напряженно-деформированного состояния активной части генератора может быть получена из трех источников: из натурного эксперимента, по расчетной теоретической модели, из численного анализа электромагнитных процессов в генераторе. Для примера приведен случай использования экспериментальных данных силы магнитного притяжения, а именно, максимального и минимального значения силы для торцевого генератора мощностью 1 кВт с одним ротором и одним статором (рис. 5). Рис. 5. Приложение силы притяжения на магниты ротора Непосредственное значение сил рассчитывалось по следующей методике. Для случая, когда известны значения максимальной и минимальной сил притяжения, значение силы в узлах определяется по формуле Fi = 2 • ( Fmax + Fmin ) + 1 • ( Fmax - Fmm) • Cos ^ , где Fi - сила притяжения i-го магнита, Н; Fmax - максимальная сила притяжения, Н; Fmin - минимальная сила притяжения, Н; P - число пар полюсов. Для случая, когда известна средняя сила притяжения, значение силы в узлах магнитов определяется по формулам F = Fo +AF • cos ( 2p1 ), AF = F0 • R • sin a, где a - начальный угол перекоса диска под влиянием технологических погрешностей; F0 - средняя сила, Н. Ограничениями для вала являются подшипниковые опоры. Существует возможность как задать их жестко, так и задать радиальную и осевую жесткость. Как показали результаты исследования, выбор способа задания подшипниковых опор существенно не влияет на суммарные деформации. Структурный анализ позволяет получить на выходе значения напряжений, деформаций и перемещений в расчетной модели. Поскольку ставится задача проверки соблюдения зазора между статором и ротором, анализируются значения суммарных перемещений под действием силы магнитного притяжения. По результатам решения суммарные перемещения для одного диска ротора (рис. 6) без учета начального перекоса диска составляют не более 0,016 мм (рис. 6, а). С учетом начального перекоса величина суммарных перемещений возрастает и составляет не более 0,048 мм (рис. 6, б). НОШ. >СШК* Л) б Рис. 6. Суммарные перемещения для одного диска ротора При анализе конструкции, состоящей из двух роторов, на валу суммарные перемещения составили не более 0,001 мм без учета начального перекоса и не более 0,048 мм с учетом начального перекоса дисков роторов. Полученная расчетная модель в среде ANSYS позволяет получить полную картину напряженно-деформированного состояния расчетной модели. Разработанный программный продукт как часть автоматизированной среды проектирования торцевого генератора существенно упрощает процесс проведения решения, поскольку требует ввода лишь входных параметров без ввода команд на языке APDL. Такое автоматизированное построение модели существенно сокращает время, затрачиваемое на построение расчетной модели, и исключает случаи некорректного написания команд для среды ANSYS. В дальнейшем этот программный продукт существенно упростит проектирование и проверку новых конструкций торцевых генераторов.

About the authors

Ekaterina Valer’evna Karpenko

Siberian Federal University

Email: cat.kras@gmail.com
postgraduate student

Natal’ya Anatol’evna Kolbasina

Siberian Federal University

Email: natalya-kolbasina@yandex.ru
Candidate of Engineering Sciences, associate professor

Denis Ivanovich Morozov

Siberian Federal University

Email: danner2@mail.ru
senior teacher

References

Supplementary files