Scientific journal
Modern problems of science and education
ISSN 2070-7428
"Перечень" ВАК
ИФ РИНЦ = 1,006

TOPOGRAPHY OF THE MAGNETIC FIELD AROUND SAMPLE SOFT MAGNETIC MATERIAL

Lankin M.V. 1 Narakidze N.D. 1 Lankin A.M. 1
1 Federal State Budget Educational Institution of Higher Professional Educational “Platov South-Russian State Poly-technic University (NPI)”
The article describes a method of obtaining the topography of the magnetic field in the vicinity of the sample magnetic material using multiple regression analysis (MRA) when tested in a semi-open magnetic system (MS). To build the model was used multiple regression analysis, the task of which is to develop an equation of such a surface in three-dimensional space, which is located at a minimum distance from the observation results. The result was obtained by the equation model of the surface in three dimensions of the sample became meaningless E330 all insignificant coefficients and the residual dispersion. Thus, the experiment confirms that the change in the magnetic field along the x axis from the sample surface of the magnetic material when tested in a semi-open magnetic system is nonlinear.
topography of the magnetic field
the magnetic system
magnetic materials
regression analysis

В [2] с помощью нелинейного регрессионного анализа была получена модель топографии магнитного поля (МП) постоянного магнита при его испытании в полуразомкнутой магнитной системе (МС). Для разработки метода, позволяющего измерять напряженность МП поля на поверхности изделий из магнитомягких материалов (МММ) при их испытаниях в полуразомкнутых МС, необходимо решить задачу построения математической модели изменения напряженности МП.

Такая математическая модель в случае неизменности зазора δ между полюсами намагничивающей системы (НС) и изделия из ферромагнитного материала (ФММ) и его геометрических размеров λ (λ – параметр, характеризующий геометрические размеры изделия из ФММ; λ = l/dmax, где l – длина изделия из ФММ; dmax – максимальный размер поперечного сечения изделия из ФММ) имеет вид , где х – расстояние от испытуемого образца, I – ток НС.

Для построения модели топографии МП изделия из МММ при помощи пакета Femm 4.0 [1] была рассчитана зависимость напряженности МП от тока НС и расстояния x от центра на поверхности изделия из МММ (сталь Э330) (рис. 1 и рис. 2) с размерами образцов 29×15×15 мм и зазором между полюсами НС и изделия из МММ 3 мм (λ=2; δ=3; λ/2δ=0,33).

Рис. 1. Семейство КР образца из стали Э330 при различных значениях расстояния х (мм)

Рис. 2. Зависимость тангенциальной составляющей напряженности Hτ образца из стали Э330 от расстояния х при различных значениях тока I

Для построения модели был использован множественный регрессионный анализ (МРА). Задача МРА состоит в построении уравнения такой поверхности в трехмерном пространстве, которая располагается на минимальном расстоянии от результатов наблюдений Нi. Для определения этого расстояния воспользуемся выражением остаточной дисперсии S2ост. [4]:

,

где n – количество значений тока I, при которых измерялась напряженность МП Нi; m – количество точек пространства, в которых измерялась напряженность МП; р – количество слагаемых в уравнении поверхности; – расчетное значение напряженности МП, полученное из уравнения поверхности .

Для определения коэффициентов уравнения регрессии bi использовалось выражение в матричном виде [3]:

, (1)

где В – вектор-столбец коэффициентов уравнения регрессии bi; Y – вектор-столбец результатов наблюдений Hi; Х – матрица факторов, причем для определения нулевого коэффициента регрессии в нее добавляется столбец единичных значений, т.е. ,

ХТ – транспонированная матрица Х, (ХТХ)-1 – матрица, обратная ХТХ.

В общем виде нелинейное уравнение регрессии можно представить как:

Обычно задачу решают за несколько шагов путем повышения степени уравнения до тех пор, пока уменьшение остаточной дисперсии остается значимым. Для расчета коэффициентов bi воспользуемся выражением (1), но, начиная со второго шага перед повышением степени полинома, будем производить замену переменных, линеаризующих функций: х12=х2, х13=х3, I12=I2, I13=I3 и так далее [4].

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

,

где – критическое значение критерия Фишера, для числа степеней свободы , и выбранного уровня значимости % [4].

После проведения нескольких итераций (j=5) была построена модель напряженности МП вплоть до пятого порядка и рассчитаны остаточные дисперсии для образца из МММ стали Э330. Результаты вычислений сведены в табл. 1.

Таблица 1 Результаты вычислений остаточной дисперсии

Образец из стали Э330

j

p

f1

f2

1

276.93

2

119

117

52.85

1.356

2

5.24

5

119

114

116.44

1.359

3

0.045

15

119

104

3.66

1.371

4

0.0123

24

119

95

1.09

1.384

5

0.0113

35

119

84

График зависимости остаточных дисперсий S2ост от количества слагаемых в уравнении поверхности p для образца из стали Э330 показан на рис.3.

Рис. 3. Зависимость S2ост от количества слагаемых в уравнении поверхности p для образца из стали Э330

Таким образом, в качестве модели изменения напряженности МП при перемагничивании образцов из стали Э330 по кривой размагничивания (КР) может быть использовано уравнение регрессии четвертой степени:

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

Таблица 2. Результаты вычислений коэффициентов bi (сталь Э330, tкр =1,661)

i

0

1

2

3

4

bi

-0.472

0.050

0.827

10.141

0.00003

tbi

-7.550

0.517

4.644

37.022

0.001

i

5

6

7

8

9

bi

0.668

0.013

-0.409

0.131

-0.0002

tbi

1.045

0.105

-0.416

0.301

-0.033

i

10

11

12

13

14

bi

-3.043

-0.048

-0.017

0.016

-0.139

tbi

-3.609

-2.530

-0.257

0.174

-0.242

i

15

16

17

18

19

bi

1.094

0.00001

2.280

0.002

0.001

tbi

0.844

0.026

8.036

2.032

0.231

i

20

21

22

23

24

bi

-0.001

0.0001

-0.002

0.026

-0.457

tbi

-0.160

0.057

-0.058

0.136

-1.047

Проверка значимости полученных коэффициентов уравнения регрессии bi произведена с использованием t-критерий Стьюдента [5]:

; ,

где cii – диагональный элемент матрицы обратной к матрице нормальных уравнений (ХТХ)-1; – погрешность коэффициента регрессии.

Результаты вычислений – коэффициентов сведены в табл. 2. Из данных табл. 2 можно сделать вывод о незначимости восемнадцати коэффициентов полученного уравнения регрессии. Поэтому согласно [5], был произведен постепенный отсев незначащих коэффициентов, у которых tbi<0,5, а затем tbi<1 и, кроме того, по критерию Фишера проводилась оценка значимости отношений остаточных дисперсий. В результате такого отсева было получено уравнение модели поверхности в трехмерном пространстве для образца из стали Э330 со всеми значащими коэффициентами и незначащим отношением остаточных дисперсий (табл. 3).

Таблица 3. Результаты вычислений коэффициентов bi (сталь Э330, tкр =1,659)

i

0

1

2

3

4

S2ост

bi

-0.422

0.034

0.791

10.196

-1.961

0.014

tbi

-17.098

7.226

8.045

411.583

-11.753

i

5

6

7

8

9

bi

-0.047

0.483

1.893

0.002

-0.268

tbi

-64.270

15.050

24.561

30.452

-18.079

После исключения незначащих коэффициентов уравнение модели для образца из стали Э330 выглядит следующим образом:

.

При подстановке bi :

На рис. 4 показано распределение напряженности МП для образца из стали Э330, построенное с помощью полученного уравнений регрессии.

Рис. 4. Регрессионная модель распределения напряженности МП стали Э330

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

Статья подготовлена с использованием оборудования ЦКП "Диагностика и энергоэффективное электрооборудование" ЮРГПУ(НПИ).

Рецензенты:

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

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