Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

CALCULATION THE MAGNETIC FIELD OF THE ACTUATOR WITH SHAPE MEMORY EFFECT BY COMBINED METHOD OF FUNDAMENTAL SOLUTIONS AND FINITE ELEMENTS

Bakhvalov Y.A. 1 Grechikhin V.V. 1 Yufanova A.L. 1
1 Federal State Budget Educational Institution of Higher Professional University «Platov South-Russian State Polytechnic University»
1156 KB
The article is devoted to the calculation of the magnetic field of an actuator with shape memory effect. The mathematical model of magnetic system of the actuator is constructed based on decomposition of a magnetic field in the linear environment for the sum of fields: a field created by currents of coils in the absence of an active element from the ferromagnetic with shape memory and at the field formed by magnetization of an element in the absence of current in coils. For the calculation of the magnetic field in the surrounding space and in the active element is used method of fundamental solutions and finite element method respectively. The article describes the computational algorithm developed on the basis of mathematical model and shows the example of its application. Accuracy of calculation of a magnetic field, sufficient for engineering tasks can be received at rather small number of fictitious magnetic charges. It allows to realize the «fast» computer program for calculation of magnetic fields of actuators with effect of shape memory that is important in problems of functional diagnostics.
actuator with shape memory effect
magnetic field
finite element method
method of fundamental solutions
computational algorithm
1. Bakhvalov Yu.A., Knyazev S.U., Scherbakova E.E., Scherbakov A.A. Modelirovanie potentsialnykh poley s primeneniem metoda tochechnykh istochnikov [Modeling of potential fields using the method of point sources]. Novocherkassk, YURGPU (NPI), 2012. 158 p.
2. Vasilev A.N., Buchelnikov V.D., Takagi T., Khovailo V.V., Estrin E.I. Shape memory ferromagnets, PHYS-USP, 2003, vol. 46, no. 6, pp. 559–588.
3. Tozoni O.V., Maergoyz I.D. Raschet trekhmernykh elektromagnitnykh poley [The calculation of three dimensional electromagnetic fields]. Kiew, Tekhika, 1974. 352 p.
4. Yavorskiy B.M. Detlaf A.A. Spravochnik po fizike dlya inzhenerov i studentov vuzov [Handbook of physics for engineers and University students]. Moscow, Nayka, 1974. 943 p.
5. Asua E., Feuchtwanger J., Garcнa-Arribas A., Etxebarria V., Orue I., Barandiaran J.M. Ferromagnetic shape memory alloy actuator for micro- and nanopositioning. Sensor Letters, 2009, no. 7, pp. 348–350.
6. Gorbatenko N.I., Grechikhin V.V., Shaikhutdinov D.V. Measuring and Actuating Devices Based on Shape Memory Ferromagnets. Metal Science and Heat Treatment, 2015, vol. 56, pp. 609–613.
7. Suorsa I., Pagounis E., Ullakko K. Magnetic shape memory actuator performance. Journal of Magnetism and Magnetic Materials, 2004, vol. 272, pp. 2029–2030.

Одним из современных направлений развития промышленных систем управления является применение интеллектуальных материалов в измерительных и исполнительных элементах. К таким материалам относятся ферромагнетики с памятью формы (ФМПФ), способные изменять геометрические размеры под воздействием электромагнитных и механических сил и запоминать состояние [2, 5, 7]. Это позволяет повысить точность преобразования одного вида энергии в другой, упростить конструкцию устройств, сократить количество составных частей, подверженных поломке и износу. Как следствие, повысить эффективность систем управления. Поэтому разработка теории, принципов построения и средств исследования устройств на основе ФМПФ является актуальной задачей.

Постановка задачи

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

Общий вид двухтактного актуатора показан на рис. 1 [6]. В нем используется распределенная намагничивающая система, содержащая n пар катушек без ферромагнитных сердечников, что позволяет повысить быстродействие, уменьшить массу и габариты устройства.

pic_9.tif

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

Математическая модель

Построим математическую модель магнитного поля одной пары катушек, между которыми находится часть активного элемента из ФМПФ, размеры которого совпадают с размерами окон катушек (рис. 2). Такую модель достаточно просто адаптировать на случай n пар катушек.

Стационарное трехмерное магнитное поле задачи описывается системой уравнений

bahvalov01.wmf (1)

с граничными условиями на границе раздела сред: bahvalov02.wmf bahvalov03.wmf Здесь µ0 – магнитная постоянная; bahvalov04.wmf – плотность тока в катушках; V+ – подобласть, занятая активным элементом (ферромагнетиком); V– – подобласть, окружающего элемент пространства, заполненная линейной средой с µ0; bahvalov05.wmf bahvalov06.wmf –магнитная индукция и напряженность магнитного поля соответственно; bahvalov07.wmf bahvalov08.wmf – нормальные составляющие магнитной индукции; bahvalov09.wmf bahvalov10.wmf – тангенциальные составляющие напряженности магнитного поля на границе раздела сред V+ и V–.

Для уменьшения размерности задачи будем далее использовать скалярный магнитный потенциал и фиктивные магнитные заряды. Это позволит применить для расчетов параметров поля в V бессеточный численный метод фундаментальных решений (МФР).

pic_10.tif

Рис. 2. Сечение исследуемой системы плоскостью y0z

Представим в V напряженность магнитного поля в виде суммы двух полей [3]:

bahvalov11.wmf (2)

где bahvalov12.wmf – напряженность магнитного поля, созданного токами катушек во всем пространстве V (V = V+ + V) при удалении из V+ ферромагнетика; bahvalov13.wmf – напряженность магнитного поля в V, созданного намагниченностью ферромагнетика при отсутствии тока в катушках.

При использовании разложения (2) в V+ результирующее поле будет определяться как разность близких величин, что приведет к большим погрешностям. Поэтому, в отличие от [5], в V+ будем определять результирующее поле bahvalov14.wmf.

Исходную систему уравнений (1) после введения магнитного потенциала j заменим тремя системами:

bahvalov15.wmf bahvalov16.wmf

bahvalov17.wmf в V; (I)

bahvalov18.wmf

bahvalov19.wmf bahvalov20.wmf в V–; (II)

bahvalov21.wmf

bahvalov22.wmf B = B(H) в V+. (III)

Решение системы (I) заменим вычислением bahvalov23.wmf путем интегрирования по объемам катушек Vc (Vc = Vc1 + Vc2), используя формулу, полученную из закона Био – Савара – Лапласа [4]

bahvalov24.wmf (3)

где bahvalov25.wmf iwc – магнитодвижущая сила, создаваемая током i в катушке с числом витков wc; Sc – площадь сечения катушки; bahvalov26.wmf – единичный вектор, направленный из точки N к точке M.

Решение системы (II) в V выполним МФР, располагая в V+ фиктивные магнитные заряды.

Решение системы (III) в V+ с нелинейной средой из ФМПФ выполним методом конечных элементов (МКЭ).

Начало декартовой системы координат поместим в геометрическом центре активного элемента (рис. 2). Учитывая симметрию магнитной системы относительно плоскости x0y, будем решать задачу в области z ≥ 0. Принимаем, что на плоскости x0y все потенциалы bahvalov27.wmf равны нулю. На боковых гранях параллелепипеда активного элемента положим bahvalov28.wmf. На грани, лежащей в плоскости x0y, φ+ = 0. На верхней грани Sh имеем следующие соотношения, вытекающие из исходных граничных условий (1),

bahvalov29.wmf (4)

где bahvalov30.wmf

bahvalov31.wmf

bahvalov32.wmf

точка Mi ∈ Sh; i = 1, 2, ..., N; N – число точек коллокации на Sh.

Вычислительный алгоритм задачи

Считаем известными геометрические размеры активного элемента: ширину Ae и длину Be (размеры по осям 0x и 0y); толщину He (размер по оси 0z), геометрические размеры катушки (ширину Ac, длину Bc, толщину Hc, толщину намотки провода Cc) и магнитодвижущую силу iwc, создаваемую токами в катушках.

На верхней грани активного элемента разместим точки коллокации Mi (i = 1, 2, ..., N). Точки коллокации Mi имеют одинаковую z-координату bahvalov33.wmf. Фиктивные магнитные заряды (количество их 2N) разместим в точках, x-координаты и y-координаты которых совпадают с соответствующими координатами точек коллокации Mi, причем bahvalov34.wmf, что следует из равенства нулю потенциалов всех точек плоскости x0y.

На рис. 3 показано положение точек коллокации M, при их числе N = 5, и фиктивных магнитных зарядов q, число зарядов 2N. Размеры a = b = 2∙10–3 м.

pic_11.wmf

а

pic_12.wmf

б

Рис. 3. Верхняя грань активного элемента Sh с точками коллокаций M(a) и расположением фиктивных магнитных зарядов q(b)

Аппликаты всех зарядов равны h для qi и –h для bahvalov35.wmf. Принимаем h = He/20.

Потенциал в точке коллокации Mi определим по формуле [1]

bahvalov36.wmf (5)

где bahvalov37.wmf

bahvalov38.wmf

(xi, yi, zi) – координаты точки коллокации Mi; (xj, yj, zj) – координаты точки, в которой находится заряд qj; bahvalov39.wmf – координаты точки, в которой находится заряд bahvalov40.wmf.

Z-составляющую напряженности магнитного поля в точке Mi, созданного зарядами, будем определять на основании (5) по формуле

bahvalov41.wmf (6)

где zj = h; bahvalov42.wmf

Система (4) с учетом соотношения (5) и (6) принимает следующий вид:

bahvalov43.wmf (7)

Величина bahvalov44.wmf системы (7) определяется по формуле (3), bahvalov45.wmf находиться из соотношения

bahvalov46.wmf (8)

Таким образом, система (7), состоящая из 2N уравнений, содержит 3N неизвестных bahvalov47.wmf. Дополнительно N уравнений, связывающих φ+(Mi) и bahvalov48.wmf, получим, решая систему (III) МКЭ.

Алгоритм решения исходной задачи состоит из следующих этапов:

1. Выбираем расположения N точек коллокаций. Формируем массивы координат точек коллокации и точек, в которых расположены фиктивные магнитные заряды. Например, для точки Mi имеем Mi(xi, yi, H/2), qi(xi, yi, h), bahvalov49.wmf.

2. Вычисляем по формулам (3) и (8) в точках коллокации значения bahvalov50.wmf и bahvalov51.wmf, i = 1, 2, ..., N.

3. Учитывая, что bahvalov52.wmf и bahvalov53.wmf, полагаем φ+(Mi) = 0 и, решая первые N уравнений системы (7), найдем qi j = 1, 2, ..., N.

4. Решая вторые N уравнений системы (7), находим bahvalov54.wmf i = 1, 2, ..., N.

5. Решая краевую задачу МКЭ – систему (III) с граничными условиями: на Sh значения bahvalov55.wmf, где M ∈ Sh, в узлах конечной элементной сетки, которые определяются по второй формуле системы (7); на грани V+, лежащей в плоскости x0y, полагаем φ+ = 0, на боковых гранях bahvalov56.wmf.

6. В тех случаях, когда магнитная проницаемость среды V+ зависит от H+, производится уточнение значений m и H+ с использованием зависимости B(H) ФМПФ.

7. Определив φ+(Mi), можно уточнить значения зарядов qi, решая систему (7).

Результаты экспериментальных исследований

Выберем для исследования систему (рис. 2), имеющую следующие параметры:

Ae = Be = Ac = Bc = 2∙10–3 м; He = 1∙10–3 м;

Hc = 2∙10–3 м; Cc = 0,5∙10–3 м.

Для оценки погрешности моделирования комбинированным методом (МФР и МКЭ) решим задачу МКЭ, при этом катушки и активный элемент поместим в центре параллелепипеда со сторонами 0,5×0,25×0,25 м, который заполним 87434 тетраэдрами. На гранях параллелепипеда задаем векторный магнитный потенциал bahvalov57.wmf. На первом этапе расчетов определяем при iwc = 103. А каждой катушки значения bahvalov58.wmf в точках коллокации Mi, а также в ряде точек на отрезке [0, zi], где zi = He/2. Далее численным интегрированием вычисляем bahvalov59.wmf с использованием соотношения (8). Затем определяем величины bahvalov60.wmf и φ+ в точках коллокаций.

Перейдем к решению задачи комбинированным методом (МФР и МКЭ).

Вариант I. Выбираем N = 5 точек коллокаций (рис. 3). В силу симметрии их положения имеем q1 = q2 = q3 = q4. Таким образом, в этом варианте имеем два неизвестных заряда q0 и q1, а также bahvalov61.wmf, bahvalov62.wmf, bahvalov63.wmf и bahvalov64.wmf.

Величины bahvalov65.wmf и bahvalov66.wmf в точках коллокации M0 и M1 определяем по формулам (3), (8). В результате имеем систему уравнений

bahvalov67.wmf

bahvalov68.wmf

Полагая φ+(M0) = φ+(M1) = 0, найдем q0 и q1. Далее определяем H+(M0) и H+(M1) при начальном значении μ(0) = 103μ0. В случае μ(H+) решаем краевую задачу в V+ МКЭ, уточняя bahvalov69.wmf φ+, q0, q1. Результаты сведем в табл. 1.

Среднеквадратичное отклонение (RMSD) вычисляем по формуле

bahvalov70.wmf

Таблица 1

Результаты расчетов при N = 5

Номер точки (i)

Расчет комбинированным методом (МФР и МКЭ)

Расчет МКЭ

Относительная погрешность

bahvalov71.wmf, А/м

bahvalov72.wmf, А

bahvalov73.wmf, А/м

φ+, А

bahvalov74.wmf, А/м

bahvalov75.wmf, %

0

176613,4

83,8

445,8

0,1

266,2

67,4

1

125161,4

57,7

311,6

0,082

212,5

46,6

Относительное RMSD для N = 5 составило 53,3 %.

Вариант II. Разместим N = 13 точек коллокаций (рис. 4) на Sh. В силу симметрии положения точек имеем четыре неизвестных заряда q0, q1, q2, q3, а также неизвестные bahvalov76.wmf и φ+(Mi), i = 0, 1, 2, 3.

Система уравнений составляется так же, как в варианте I. Результаты сведены в табл. 2. Относительное RMSD для N = 13 составило 9,1 %.

Замечание. Результаты, приведенные в табл. 1, 2, получены при μ = 103μ0.

pic_13.wmf

Рис. 4. Расположение точек коллокации N = 13

Таблица 2

Результаты расчетов при N = 13

Номер точки (i)

Расчет комбинированным методом (МФР и МКЭ)

Расчет МКЭ

Относительная погрешность

bahvalov77.wmf, А/м

bahvalov78.wmf, А

bahvalov79.wmf, А/м

φ+, А

bahvalov80.wmf, А/м

bahvalov81.wmf, %

0

176613,4

83,8

306,5

0,089

266,2

15,1

1

125161,4

57,7

238,9

0,070

212,5

12,4

2

161570,5

69,5

233,0

0,099

238,0

–2,1

3

157664,0

74,0

265,1

0,055

246,4

7,6

Выводы

Предложена математическая модель минимальной размерности и алгоритм для расчета магнитного поля актуатора на основе ФМПФ. Модель и алгоритм основаны на разложении магнитного поля в линейной среде, окружающей устройство, на два поля: bahvalov82.wmf – напряженность магнитного поля, созданная токами катушек при отсутствии активного элемента из ФМПФ и поля bahvalov83.wmf, создаваемого намагниченностью ФМПФ при отсутствии катушек с током. Для определения магнитного поля в окружающем пространстве bahvalov84.wmf и активном элементе bahvalov85.wmf используются МФР и МКЭ соответственно. При этом магнитное поле bahvalov86.wmf вычисляется по формуле, вытекающей из закона Био – Савара – Лапласа.

Достаточная для инженерных задач точность расчета параметров магнитного поля может быть получена при относительно небольшом числе фиктивных магнитных зарядов, что характерно для МФР [1]. Так при 26 зарядах погрешность расчета не превышает 9 %. Как следует из рис. 3, дальнейшее увеличение зарядов приведет к уменьшению погрешности. При 34 зарядах следует ожидать менее 5 %.

При учете симметрии магнитной системы размерность задачи существенно сокращается: при 10 зарядах число неизвестных равно двум (q0, q1), при 18 зарядах – трем (q0, q1, q2), при 26 зарядах – 4 (q0, q1, q2, q3). Это позволяет реализовать «быструю» компьютерную программу для расчета магнитных полей актуаторов на основе ФМПФ, что важно в задачах диагностики при прогнозировании состояния объекта и при анализе неисправностей в реальном времени.

Результаты работы получены при поддержке гранта РФФИ № 14-08-01288 «Разработка теории натурно-модельных испытаний измерительных и исполнительных систем, построенных на основе ферромагнитных материалов с эффектом памяти формы».

Рецензенты:

Савелов Н.С., д.т.н., доцент, профессор кафедры «Автоматика и телемеханика», ФГБОУ ВПО «Южно-Российский государственный политехнический университет имени М.И. Платова», г. Новочеркасск;

Седов А.В., д.т.н., доцент, профессор кафедры «Автоматика и телемеханика», ФГБОУ ВПО «Южно-Российский государственный политехнический университет имени М.И. Платова», г. Новочеркасск.