CC BY
18
УДК 681.31:620.179.14 Дедегкаев А.Г., Степанов А.Л.
ОРИЕНТИРОВАННАЯ ДЛЯ САПР УНИВЕРСАЛЬНАЯ МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ВИХРЕТОКОВОГО ПРЕОБРАЗОВАТЕЛЯ. ЧАСТЬ 1. РЕЗУЛЬТИРУЮЩЕЕ ПОЛЕ ОБМОТКИ ВОЗБУЖДЕНИЯ
Dedegkaev A.G., StepanovA.L.
ORIENTED FOR CAD UNIVERSAL MATHEMATICAL MODEL OF THE EDDY CURRENT PROBE. PART 1. THE RESULTING FIELD OF THE DRIVE WINDING
Работа относится к области САПР методов и средств неразрушаю-щего электромагнитного контроля проводящих плоских изделий. Рассмотрен вопрос построения математической модели обмотки возбуждения трансформаторных вихретоковых преобразователей линейно-протяженной формы с поперечным сечением в виде выпуклого многоугольника. Это позволяет создать универсальную модель, в которую включены
70
модели более простой формы в поперечном сечении, как частный случай.
Ключевые слова: математическая модель, вихретоковый преобразователь, обмотка возбуждения, вектор-потенциал электромагнитного поля.
The work relates to the field of CAD methods and means of nondestructive electromagnetic testing of conductive flat products. The problem of constructing mathematical model of the drive winding eddy current probe linearly extended form is polygon in cross section. This allows you to create a universal mathematical model, which includes a model of a simpler shape in cross-section, as a special case.
Key words: mathematical model, eddy current probe, the drive winding, the vector potential of the electromagnetic field.
Внедрение автоматизированного проектирования (САПР) в создание контрольно-диагностического оборудования позволяет сократить время разработки и повысить надежность проектных решений. Важной составной частью автоматизированного проектирования является использование готового математического обеспечения - математических моделей (ММ) [1].
В данной работе продолжены результаты исследований по созданию ориентированной для САПР универсальной ММ трансформаторного вихре-токового преобразователя (ВТП) линейно-протяженной формы [2]. В развитие указанных исследований рассмотрена базовая ММ ВТП указанной формы, включающая обмотку возбуждения (ОВ1) с поперечным сечением в виде плоского выпуклого многоугольника (рис.1), и измерительную обмотку (ИО1), расположенные (в пустоте 2) между проводящими полупространствами (1 и 3). Пренебрегая влиянием коротких сторон обмоток, ОВ1 и ИО1 представлены в виде двух бесконечно длинных половин: правой (B) и левой. (A). Упрощающие анализ допущения для такой формы обмоток изложены в работе [3].
Целью и задачей работы являлось построение математических выражений, описывающих в пустоте (2) начальное (первичное) электромагнитное поле ОВ1 указанной формы и вторичное поле вихревых токов, протекающих в проводящей среде. Эти выражения составляют основу ММ ОВ и в развернутом виде весьма громоздки. Поэтому приведем их в обобщенном виде. Так начальное электромагнитное поле описывает вектор-потенциал [3,4]
2 2 г IVcW, (x-x2lj ) + (У-У2Ь ) ^ An, = f ln--^^-dSi 1k; (1)
4*S1 S\ (x - x^ )2 + (y - y^ )2 1 ; ()
вторичное поле вносимого векторного потенциала определяет выражение [3,4]
4 ^ 0<?„[ Е ) еХх + Е2&Х) е—Хх /
4 вн к = —^ОГ Ш-ХЕ-с°8 х(у " у1-х
1
1
Ез(-х)ехХ + Е4(&) е—Хх X Е
е°8 Х(у —
(? — у2-у)
d X 1к
(2)
Рисунок 1 - Схема расположения ВТП линейно-протяженной формы
где х,у - координаты системы Оху, I - комплекс действующего значения тока в проводниках ОВ1, ^о - магнитная постоянная, ^ - число проводников, уложенных на каждой половине ОВ1, Sl - площадь поперечного
—
сечения каждой половины обмотки, 1к - орт по оси 2 , -, х - некоторые взаимно ортогональные переменные интегрирования, имеющие область определения в пределах площади поперечного сечения обмотки, у^ , х^ ;
х2^ , - координаты срединных точек на площади поперечного сечения
ОВ1, определенные в системе декартовых координат Оху, Е\(-х),Ег(-х),Еъ(-х), Еа(-х),ХЕ комплексные величины, включающие постоянные интегрирования (их развернутые выражения приведены в [3,4]). Вектор-потенциалы направлены вдоль оси 2 . Будем рассматривать их комплексные модули.
В процессе анализа и вычисления (1) и (2) установлено: •интегралы (1) и (2) можно представить, в общем случае, в виде повторных интегралов, поскольку:
" у^х , х^с; ^х, у^х можно выразить через -, х с помощью формул линейного преобразования координат [5];
- при замене переменных X и у на - их не происходит искажение площади поперечного сечения обмотки и функциональный определитель прямого отображения (якобиан) [6] в данном случае равен
^ (4; х) =
дх дх д4 дх
ду дУ дХ
= +1-
Получим первоначальную форму представления (1) и (2)
:Б еЛ „Л
А 0к =
/
V 0 4к^1
4Б хб (4) ^ 1
I 11п[оЬ(х)]¿х- | ^ | 1ПНХ)]
х л (4)
Б
х Б (4)
4
Л
х Л (4)
(3)
А
1 V 0 ^
Л .Л
хл (4)
4
Б „.Б
хБ (4)
вн к
2тс£
й4 |уп1(х)^Х-|d4 |уп2(хМх
10
Л Л Б Б
4л ХЛ (4) 4
х4 (4)
(4)
2 2 4
где ch(x),zп(x),vпl(x),vп2(x) - соответственно числитель и знаменатель в (1), а также 1-ое и 2-ое слагаемые в (2). Они соответствуют выражениям для левой (Л) и правой (Б) половины обмотки; 4 гЛ(Б) - координата пока неизвестной точки на площади поперечного сечения обмотки (на любой половине обмотки), как предел внешнего интеграла в (3) и (4); х Л(Б) (4) - уравне-
1
ние пока неизвестной прямой также как предел внутреннего интеграла в (3) и (4);
•каждое слагаемое в (3) и в (4) можно представить в виде алгебраической суммы простых интегралов (это заключение обосновано анализом (3) и (4) при произвольном числе сторон указанного многоугольника, используя приемы взятия внутреннего интеграла и свойство аддитивности для простых определенных интегралов [6]). Их число равно удвоенному числу сторон многоугольника (2п) на каждой половине обмотки. Так (3) примет вид
А 0к =
1 V 0 4тс^
21 I
22
В
[СН(х) ]Х12
(4)
х22 ¿4+ I
X
В
[СН(х) ]х 23
(4)
23
¿4 +...
...+
2п
I
X
В
[СН(х) ]х п-1, п
-1 п(4)
Х21 ¿4+ I
2, п-1
Чп
В
[СН(х) ]х п1
(4)
¿4
11
12
.Л
[Ш (х) ]Х12
42
+ I
[2м (х) ]х л3(4)
х13
х1п
¿4 +... + |
X
^ (х) ]х л-1, п (4)
1, п-1
11
¿4+ I
Чп
^ (х) ]х л1(4)
(4)
¿4
¿4 +
;(5)
приведем обобщенную форму интеграла (4)
w
11 V,
внк
2^ К д
41
I
.Л
[ум (х)]х12
(4)
х12
12
¿4+ I
Л
[ум (х)]х 23
(4)
х13
¿4 +.
<
>
4
4
2
ж
1
3
Ж
. +
J
x
[vn (x)F n (4)
1, n-1
41
d4 + J
1n
A
[VN (x)]x nl
(4)
d4
21
22
B
[VN (x)]Xl2
(4)
d4 +
x22
+ J
[VN (x)FB3(4)
23
x2n
d4 +... + J
x
[VN (x)]x«-l,n(4)
2,n-1
21
d4+ J
x2n
[VN (x)]X«1(4)
d4
>d^, (6)
,b(a)
b(a),
b(a),
b(a),
где [СИ(х); [Ш(х)]хм+1®; [^(х)рм+1®; [^(х)рм+1® - первообразные внутренних интегралов (3) и (4), определенные одним из пределов
х В(+Д(4) +(В)
интегрирования '''; Хц+{ - уравнение прямой проходящей через вершины многоугольника +/+/+1 (В/В/+д), как уравнение прямой с угловым коэффициентом [5];
•произвольный член в (5) и (6) представляет собой линейный табличный интеграл, взятый вдоль / - той стороны правого многоугольника и / -той стороны левого многоугольника соответственно. При этом
хшахД/ ^ хш1пД и хшах,2/ ; хтт,И- с00тветственн0, большая и меньшая по величине координата двух вершин, ограничивающих / -тую сторону левого и правого многоугольника соответственно;
•знаки слагаемых определяют соотношения между координатами вершин многоугольников и координатами основания перпендикуляра, опущенного из вершины на диагональ многоугольника [7].
Результирующее поле определяем наложением рассмотренных полей
Adk = (a0k + Авнk)1k •
(7)
В качестве примера рассмотрим порядок (алгоритм) вычисления вектор-потенциалов для ММ ОВ с поперечным сечением в виде ромба (рис.2):
Считают заданными координаты вершин в системе oxy' A^x| 1; y' 1);
A2 (x12 ; У12) ; A3 (x13; У13); A4 (x14; У14 ); B У21); B 2 fe ; У22 k B3 (x'3 ; У'3) B4 (x24;У24) и, следовательно . Также считают заданным число проводников в пределах каждой половины ОВ1;
• определяют координаты вершин ромба в системе Oxy (рис.2) по формулам перехода [5]
X = xQ + x - cos a - y' sin a; y = yQ + x - sin a + y' cos a,
где xo,Уо - координаты центра O' системе Oxy;
• рассчитывают параметры уравнения прямой, проходящей через соседние вершины ромба: в общем случае, для стороны A- A-+1
(i е {1,2,3,4}), как уравнение прямой с угловым коэффициентом
x = -(k 1i 4 + d1i) ,
= (*i, i+1 - xi i)sin а +U, i+1 - У1 i)cos а .
1i (xi, i+i - xi i)cos а - к i+i - У1, i)sin а '
dii =
(x o-xi, i К i+i- У, i )-(xi, i+i- i )y о - У, i)
(xi, i+i- xi, i)cos а-(У1, i+i- >i i)sin а '
Рисунок 2 - Схема ОВ с поперечным сечением в виде ромба Аналогично определяют коэффициенты к2г и ^ • уравнения прямой
B, Bi+i
х =
-(k 2i^ + d2i) •
k 2 i ~
(x2, i+1 - x2 i)sin а +(У2, i+1 - У2 i)
icos а
(x2, i+i - X2 i)cos а - (y2, i+1 - У2, i)sin а '
d 2 i =
(x 0 - x2, i)(y2, i+1 - y2, i )-(x2, i+1 - x2, i)(y 0 - y2, i) (x2.,+i-X2, ,)cosa-(y2.,+i - y2,)sin а ;
фиксируют количество слагаемых для данной формы поперечного сечения обмотки, определяющих модуль А ок и А к. Каждое слагаемое подобно
приведенным в (5) и (6), их число равно 8 (для каждой стороны ромба Л1Л 2 Л 3Л 4 и Б1Б 2Б 3Б 4 );
Рисунок 3 - Схемы ОВ с различным профилем в поперечном сечении
Используют полученное развернутое выражение первообразной в (1) и (2) для каждой стороны А/ А/+д и В/В/+д многоугольника [7], вычисляют величину модуля каждого слагаемого в формулах (5) и (6); определяют знак каждого слагаемого в (5) и (6). Фиксируют вершины ромба (на каждой половине обмотки), имеющие максимальную и минимальную координату х . Определяют уравнение прямой, проходящей через эти вершины и координаты основания перпендикуляра, опущенного из каждой вершины на эту прямую. Сопоставляя координату у вершины ромба и координату у основания перпендикуляра по известному правилу [7], определяют знаки этих слагаемых; алгебраически суммируя указанные слагаемые, определяют величи-
ну А0£, А ^, и вектор-потенциал результирующего поля, создаваемого ОВ1
(7). Аналогичным образом рассчитывают поля, создаваемые линейно-протяженными ОВ с поперечным сечением в виде: прямоугольника (рис.3а), идеально тонкостенной ОВ (рис. 3б) и идеально короткой (спиральной) ОВ (рис. 3в), а также ММ ОВ многопроводных со средоточенными размерами в поперечном сечении (рис. 3г) и в виде нитевидных двухпроводных линий (рис.3д).
Алгоритм расчета таких форм обмоток проверен в [3]. Указанная терминология ОВ линейно-протяженной формы собрана и предложена авторами в результате литературного обзора и собственного опыта (Таблица 1).
Таблица 1 - Варианты ММ линейно-протяженных ОВ1 и их аббревиатура
Русская аббревиатура обозначения Abbreviation marks in English
Обмотка возбуждения вихретоково-го преобразователя [8] (ОВ) Drive winding of eddy current probe [8] (DW)
Нитевидная двухпроводная линия (НДЛ) Filiform two-wire line (TWL)
Обмотка многопроводная сосредоточенная в поперечном сечении, обмотка в виде двухпроводной многожильной линии (ОМС) Multi-wire winding centered in the cross section (MWC)
Линейно-протяженная идеально тонкостенная обмотка (ТО) Linearly extended perfectly thin walled winding (TW)
Линейно-протяженная идеально короткая обмотка (спиральная обмотка) (СО) Spiral winding (SW)
Линейно-протяженная обмотка прямоугольная в поперечном сечении (ОПС) Linearly extended winding is rectangular in cross section (WRS)
Линейно-протяженная обмотка ромбовидная в поперечном сечении (ОРС) Linearly extended winding diamond-shaped in cross section (WDS)
Линейно-протяженная обмотка эллиптическая в поперечном сечении (ОЭС) Linearly extended winding is elliptical in cross section (WES)
Линейно-протяженная обмотка круглая в поперечном сечении (ОКС) Linearly extended winding round in cross section (WRS)
Линейно-протяженная обмотка с формой поперечного сечения в виде многоугольника (ОМС) Linearly extended winding is polygon in cross section (WPS)
Вывод.
Таким образом в данной работе построена универсальная ММ ОВ1 трансформаторного ВТП линейно-протяженной формы с поперечным сечением в виде выпуклого многоугольника, которая включает ОВ разнообразные по профилю поперечного включения, как частные случаи. При этом:
• получены расчетные выражения, описывающие поле этой обмотки в отсутствие проводящей среды (начальное или первичное поле);
• построены расчетные выражения, описывающие поле вихревых токов в проводящей среде в области пространства, где расположена ОВ (вторичное поле);
• предложен алгоритм расчета этих полей (в словесной форме) и результирующего поля обмотки возбуждения, который конкретизирован на примере расчета поля обмотки с поперечным сечением в виде ромба;
• предложена терминология обозначения подобных по форме обмоток возбуждения.
Библиографический список:
1. Норенков И.П. Основы автоматизированного проектирования. - М.: Изд-во МГТУ имени Н.Э. Баумана, 2002. - 336 с
2. Степанов А.Л., Дедегкаев А.Г. Ориентированная для САПР универсальная ММ ВТК плоских проводящих сред. Часть 2 Решение задач для построения ММ. Труды Северо-Кавказского горно-металлургического института (ГТУ). Выпуск 21. Владикавказ, 2014. С. 112-121.
3. Степанов А.Л., Воронин П.А. Математическая модель линейно - протяженного вихретокового преобразователя для контроля плоских многослойных проводящих сред. Деп. в ВИНИТИ 06.08.97. № 2619 - В97.
4. Воронин П.А., Степанов А.Л. Математические основы автоматизированного проектирования режимов контроля и параметров вихретоковых преобразователей в однородных и неоднородных электромагнитных полях. Известия вузов. Сев.- Кав. регион. Техн. науки. Специальный выпуск «Математическое моделирование и компьютерные технологии». 2004. С. 89-93.
5. Корн Г., Корн Т. Справочник по математике для научных работников и инженеров. - М.: Наука, 1974. - 832 с.
6. Будак Б.М., Фомин С.В. Кратные интегралы и ряды. М.: Наука, 1965. 608 с.
7. Степанов А.Л., Воронин П.А. Математическая модель обмотки возбуждения вихретокового преобразователя линейно-протяженной формы с поперечным сечением в виде многоугольника. Деп. в ВИНИТИ 26.03.99 № 925 -В99.
8. ГОСТ 24289-80. Контроль неразрушающий вихретоковый. Термины и определения.
