
Расчет напряженности магнитного поля плоского витка с током (80
..pdf4 4 |
Загрядцкий В.И., Кобяков Е. Т. |
«ЭЛЕКТРИЧЕСТВО» № 8/2006 |
Расчет напряженности магнитного поля плоского витка с током
ЗАГРЯДЦКИЙ В.И., КОБЯКОВ Е.Т.
Предложен универсальный алгоритм для прибли A versatile approximate algorithm for calculating the женного расчета напряженности магнитного поля magnetic field developed by a planar turn with a current плоского витка с током. Напряженность определя through it is proposed. The algorithm allows thefieldde
ется в плоскости витка произвольной геометриче |
veloped by turn of an arbitrary shape in its plane to be de |
ской формы. |
termined. |
К л ю ч е в ы е с л о в а : магнитное поле, напря |
Key w o r d s : magnetic field, intensity, current, |
женность, ток, плоский виток, изотропная среда |
planar turn, isotropic medium |
В статье [1] получены аналитические выраже ния для напряженности магнитного поля плоско го витка с током, состоящего из двух радиальных и двух дуговых участков. Выражения позволяют определить напряженность магнитного поля в точках, лежащих в плоскости витка, расположен ного в однородной изотропной среде. Задача ре шалась на основе закона Био-Савара [2]. При этом исходная зависимость напряженности с/Я в точке М плоскости от тока / в элементе dl контура витка была преобразована к виду, удобному для решения поставленной задачи, что позволило пу тем интегрирования найти аналитические выраже ния напряженности магнитного поля в точке М от тока / радиальных и дуговых участков контура.
A R |
|
Для напряженности поля Нм |
в точке М |
(рис. 1) от тока / радиального (прямолинейного)
участка АВ аналитическое выражение известно [3, 4]. В []] оно представлено в виде:
^ = 4 ^ ( s i n y l + S i n y 2 ) > |
(D |
где до ~ длина перпендикуляра МК, опущенного из точки М на направление АВ; у] и у2 — углы, отсчитываемые от перпендикуляра МК к лучам МА и MB соответственно. Аналогично и выраже ние для точки М\, расположенной по другую сто рону от направления АВ (рис. 1):
нщ--^Шг',+Л,г'1У |
(2) |
Зависимости (1) и (2) характеризуются одина ковой структурой и могут быть представлены в виде единого аналитического выражения:
«ЭЛЕКТРИЧЕСТВО» № 8/2006 Расчет напряженности магнитного поля плоского витка с током 45
<~t \ла^ |
(sin у , + s i n y 2 ) , |
(3) |
|
где £=1, если точка М, в которой определяется напряженность магнитного поля, находится слева от отрезка АВ по направлению тока, и £=—1 при расположении точки М справа от направления АВ. При этом в формуле (3) целесообразно дать еди ное обозначение углов, поскольку они в обоих случаях образуются единообразно.
А
Mi
Рис. 1
Выражение для напряженности в точках плос кости от тока, протекающего по дуговым участкам контура, содержит эллиптические интегралы [1], что значительно осложняет расчет.
В настоящей статье предлагается алгоритм приближенного расчета напряженности магнитно го поля витка, который позволяет значительно уп ростить расчет и сократить его время. Алгоритм разработан применительно к плоскому витку про извольной геометрической формы. Некоторые виды витков, встречающиеся в задачах электро техники, показаны на рис. 2. Для решения постав ленной задачи криволинейные участки витка (рис. 2,а-в) аппроксимируются набором прямоли нейных отрезков.
На рис. 3 приведены обозначения геометриче ских параметров, соответствующих выделенному элементу AjBj. Каждому элементу присвоен но меру (у= 1, 2 п) при их общем числе п. Начало
а) |
б) |
Рис.3 |
|
элемента обозначено Aj, а его конец |
Bj. Все па- |
раметры, относящиеся к элементу у, имеют соот ветствующий индекс у (на рис. 3 в целях упроще ния эти индексы отсутствуют). Координаты гра ничных точек контура и точек, в которых опреде ляется напряженность, указаны в полярной систе ме координат с началом в точке О. Выбор этой точки достаточно произволен и подчинен сообра жениям удобства расчета для каждой конкретной задачи. Полярные координаты точки Kj (основа ние перпендикуляра МК:=а() обозначены как
pi и (pi. Длина лучей MAj и MBj обозначена соот ветственно через ау и bj. Координату точки Щ на направлении AjBj относительно точки Aj обозна чим через S- = А-К •, а координату точки Bj отно сительно точки К: — чt;=KKJBJ: :. . Заметим, что
значения SJ и tj могут быть положительными и от рицательными. Если точка Kj находится между точками Aj и Bj, то Sj и tj положительны. Если точ ка Kj находится перед точкой Aj (по направлению
тока), то S: <0, a t • >0. Если же точка К,- находит- |
|||
j |
j |
J |
|
ся перед точкой Bj, то |
s- >0, а / < 0 . Знаки этих |
||
параметров определяют |
и знаки углов у. и |
yl, |
|
что должно быть учтено в расчете. |
|
||
Пользуясь |
обозначениями, приведенными |
на |
рис. 3, получаем следующие геометрические соот ношения, необходимые для построения расчетно го алгоритма:
|
|
°г\рМ |
+PAJ ~2PMPAJ |
C0S(<PM |
~<PAJ >•" <4> |
в) |
г) |
bj=\p2M |
+ Р \ . ~2PMPBJ |
COS(PBJ |
-Ч>МУ> (5) |
|
Рис. 2 |
||||
|
|
|
|
|
46 Загрядцкий В.И., Кобяков Е. Т. «ЭЛЕКТРИЧЕСТВО» № 8/2006
/ 2 + д 2 _ л 2 |
|
Для определения |
величины р ^., |
пользуясь |
|
J . |
(6) |
обозначениями рис. 3, можно |
получить |
формулу |
|
21 : |
|
рв |
s\n(<pB |
-<pA ) |
|
|
|
|
|
|
|
|
|
(7) |
Siny{ = Sj/Oj-, |
(8) |
||||
sin- |
= |
0 " |
/ |
' * / ' |
(9) |
1^2 |
|
|
|
||
CQ = fl • COS у j7 = |
|
6 • COS у 2 , |
(10) |
где / • = /1 • i? • — длина участка А • 5 • контура витка. Формулы (4)—(10) справедливы как для точек
Л/, расположенных внутри контура, так и для то чек М\ снаружи контура.
Напряженность магнитного поля определяется по формуле (3) для каждого участка контура, а ре зультирующая напряженность должна быть найде на в соответствии с принципом суперпозиции пу тем алгебраического суммирования результатов, полученных для отдельных участков.
Для организации вычислительного процесса алгоритм расчета может быть представлен в таб личной форме. Определение напряженности маг нитного поля витка упрощается, если стороны витка состоят из прямолинейных участков (рис. 2,г), В этом случае отпадает процесс аппрок симирования.
Остановимся далее на определении коэффици ента £ в формуле (3). Как уже отмечалось в пояс нении к формуле (3), он может принимать два значения: 1 или -1 — в зависимости от направле ния тока / и расположения точки М, в которой определяется напряженность. Указанный способ определения коэффициента £ удобен при «руч ном» счете. Однако при применении ЭВМ целесо образно указать алгоритм его определения, необ ходимый для программирования вычислительного процесса.
Будем различать участки / и // аппроксимиро ванного контура (рис. 3), граничные точки кото рых обозначены через С и D. На участке / ток i направлен от точки С к точке D, а на участке // — от точки D к точке С.
Введем в рассмотрение точку Lj (рис. 3), в кото рой направление элемента,/ контура пересекает луч ON, проведенный к точке N, в которой определяет ся напряженность поля. Соответствующие ради альные расстояния от начала координат О до точек Lj w К та рис. Ъ обозначены, через р L и p N .
|
|
|
|
,(П) |
LJ |
Aj |
PA J.S'm{-VN'VА |
J)~РВ "J.^(tN~PВ |
~J ) |
WPAJ'PAJ |
"PBJ'VBJ |
полярные координа |
||
ты точек начала Aj "и конца Bj элемента j |
контура |
витка с током; tpN — угловая координата точки N.
Вычислив значение р, |
по формуле (11), его |
|||
следует сопоставить со значением р N. |
||||
Для |
участка / контура |
(рис, 3) |
имеем: |
|
|
£ у = - 1 при |
pLj<pN; |
( ] 2 ) |
|
|
£ := + 1 ПРИ PL |
>pN. |
||
|
J |
J |
|
|
Для |
участка //: |
|
|
|
|
£ у = - 1 при |
pL.>pN; |
( ] 3 ) |
|
|
£у = + 1 при |
pL |
<pN. |
|
При |
совпадении точек Lj и |
N, |
т.е. в случае |
|
|
|
А.-В: |
|
|
PL=PN> |
слеДУет принять Яд/ |
=0. |
В качестве примера рассмотрим определение напряженности магнитного поля в точках М и N плоского витка, образованного двумя радиальны ми и двумя дуговыми участками (элемент обмотки торцевого электродвигателя). Аппроксимируя этот виток прямолинейными отрезками, получаем многоугольник, представленный на рис. 4.
Соответствующая длина элементов контура: /, =/2 = /,0 = /9 = 2 • 13sin 7,5°= 3,3937 см;
/3 =/8 =4см;
/4 = /5 = /6 = /? = 2 • 9 sin 7,5°= 2,3495 см.
Точка М расположена в плоскости симметрии магнитного поля, поэтому для получения оконча-
/^jJ-iUUio
Рис.4
«ЭЛЕКТРИЧЕСТВО» № 8/2006 |
Расчет напряженности магнитного поля плоского витка с током |
47 |
Таблица I
Номер
контура j полярные координаты
РАГ |
Ш |
<РА> г Р а д . |
1 |
13 |
0 |
2 |
13 |
15 |
3 |
13 |
30 |
4 |
9 |
30 |
5 |
9 |
15 |
|
Исходные геометрические |
параметры |
|
|
|
|
граничных точек контура |
длина |
полярные |
координаты |
длина лучей |
МА: и MB;, |
|
|
|
элементов |
точки М |
см |
|
|
рв., см |
<рв , град. |
контура lj, |
|
|
|
|
см |
РМ' с м |
<РМ> ГРад- |
aJ |
bj |
||
13 |
15 |
3,3937 |
11 |
0 |
2 |
3,7075 |
13 |
30 |
3,3937 |
— |
— |
3,7075 |
6,5051 |
9 |
30 |
4 |
— |
— |
6,5051 |
5,5251 |
9 |
15 |
2,3495 |
— |
— |
5,5251 |
3,2782 |
9 |
0 |
2,3495 |
- |
- |
3,2782 |
2 |
Таблица 2
Номер |
|
|
Дополнительные |
геометрические параметры |
|
|
контура j |
линейные |
координаты |
sin y{ |
sin yJ2 |
длина луча MKj |
?у |
|
AKj и KBj, см |
|
|
Q Q , CM |
|
|
|
SJ |
Ч |
|
|
|
|
1 |
0,26107 |
3,13263 |
0,130537 |
0,84495 |
1,382887 |
1 |
2 |
-2,51265 |
5,90635 |
-0,677356 |
0,907953 |
2,727406 |
1 |
3 |
3,47372 |
0,52628 |
0,533997 |
0,0952522 |
5,499990 |
1 |
4 |
5,38426 |
-3,03479 |
0,974506 |
-0,925746 |
1,239625 |
J |
5 |
2,61052 |
-0,26105 |
0,796325 |
-0,130525 |
1,982888 |
1 |
тельного результата достаточно рассмотреть влия ние участков / - 5 с последующим удвоением сум мы значений напряженностей, найденных от тока этих элементов.
Граничную точку С совместим с точкой Ay, a точку D — с точкой ВТ, (рис. 4).
Составим расчетные таблицы для определения напряженности магнитного поля в точке М. В табл. 1 указаны значения известных геометриче ских параметров и параметров, вычисленных по формулам (4) и (5).
В табл. 2 приведены параметры, рассчитанные по формулам (6)—(13).
Напряженность магнитного поля в точке М на
ходим, пользуясь |
принципом |
суперпозиции: |
|||
' |
^,5 |
sin• |
yj"j + sin• |
УЛ |
2 = --2,13204. |
yJ2 |
|||||
Я ^ ~ 4 л " |
У=1 |
|
< |
|
4л |
Соответствующее значение, найденное по фор мулам [1] с использованием эллиптических инте
гралов: Н ^ = — 2,11776. Сравнивая результаты,
И—Я
получаем: |
м - ^ 1 0 0 % 2=0,674%. |
|
Н М |
Определенный интерес представляет определе ние напряженности магнитного поля в точке /V
(рис. 4), расположенной за пределами контура витка с током. Результаты вычислений представ лены в табл. 3 и 4.
Напряженности магнитного поля в точке JV ОТ тока / каждого элемента./ контура витка находим по формуле (3) с учетом значения коэффициента £ •:
HANJB{ |
= /-0,0399027; |
НNА2В2 |
_ |
4лi 0,143534; |
|
|
4л |
|
|
|
|
# £ * з =--/-0,755212; |
Я АУА |
=/-0,0881304; |
|||
/V |
4л |
N |
|
4л |
|
|
|
|
|||
Н*&=±- |
0,00434386; |
Н ^6 *6 |
= - -/- 0,00649975; |
||
N |
4л |
' |
N |
|
4л |
|
N |
|
|||
HpBl |
|
|
• AQBo |
|
|
= --^0,00799604; H'tf"* = -^0,02073792; |
|||||
HAJB<) =-/-0,0138341; |
И ^о*н> =-/-0,0205659. |
||||
/v |
4л |
|
" |
|
4л |
Результирующая напряженность поля в точке N |
|||||
|
|
10 |
|
|
|
|
Н N |
|
|
|
|
|
|
7=1 |
|
|
|
Находим |
отношение: НN |
/ Нм =-0,20575. |
|||
|
|
|
|
|
М |
Значение того же отношения, найденное в [1], равно —0,208. Расхождение составляет 0,96%.
Предложенный приближенный алгоритм рас чета может быть рекомендован в тех случаях, ко гда точные методы связаны со сложными и дли тельными вычислениями, а точность результата
48 |
|
|
|
|
Загрядцкий В,И., Кобяков Е.Т. |
|
«ЭЛЕКТРИЧЕСТВО» № 8/2006 |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 3 |
|
Номер |
|
|
|
|
Исходные геометрические параметры |
|
|
|
|
|||||
элемента |
полярные координаты граничных точек контура |
длина |
полярные координаты |
дайна лучей NA; и NB/, |
||||||||||
контура j |
||||||||||||||
|
|
|
|
|
|
|
элементов |
|
точки N |
|
см |
|
||
|
|
|
|
|
|
|
контура lj, |
pN, |
см |
<pN, град. |
|
|
|
|
|
|
РА-< с м |
<РЛ-' гРад- |
Рд.. см |
<Рв.< г Раа. |
см |
aJ |
|
bJ |
|||||
1 |
|
13 |
0 |
13 |
15 |
3,3937 |
|
11 |
40 |
8,42088 |
|
5,549411 |
||
2 |
|
13 |
15 |
13 |
30 |
3,3937 |
|
— |
— |
5,549411 |
|
2,888768 |
||
3 |
|
13 |
30 |
9 |
30 |
4 |
|
— |
— |
2,888768 |
|
2,647276 |
||
4 |
|
|
9 |
30 |
9 |
15 |
2,3495 |
|
— |
— |
2,647276 |
|
4,748795 |
|
5 |
|
9 |
15 |
9 |
0 |
2,3495 |
|
— |
— |
4,748795 |
|
7,093883 |
||
6 |
|
|
9 |
0 |
9 |
-15 |
2,3495 |
|
— |
— |
7,093883 |
|
9,403822 |
|
7 |
|
|
9 |
-15 |
9 |
-30 |
2,3495 |
|
— |
— |
9,403822 |
11,587924 |
||
8 |
|
|
9 |
-30 |
13 |
-30 |
4 |
|
— |
— |
11,587924 |
|
13,86298 |
|
9 |
|
|
13 |
-30 |
13 |
-15 |
3,3937 |
|
— |
— |
13,862298 |
11,223062 |
||
10 |
|
13 |
-15 |
13 |
0 |
3,3937 |
|
- |
- |
11,223062 |
|
8,420088 |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 4 |
|
Номер |
|
|
|
Дополнительные геометрические параметры |
|
|
ч |
|||||||
контура j |
|
линейные |
координаты |
sin y\ |
sin yJ2 |
длина луча NKj |
|
|||||||
|
|
|
|
AKj и KBj, см |
|
|
|
|
fljj, |
CM |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
SJ |
Ч |
|
|
|
|
|
|
|
|
|
|
1 |
|
|
7,607104 |
-4,213404 |
0,903362 |
-0,759252 |
3,611538 |
|
1 |
||||
|
2 |
|
|
5,004594 |
-1,610894 |
0,901824 |
-0,55764 |
2,397921 |
|
I |
||||
|
3 |
|
|
2,167114 |
1,832886 |
0,750186 |
0,692367 |
1,910131 |
|
-1 |
||||
|
4 |
|
|
-2,133029 |
4,482499 |
-0,805745 |
0,943923 |
1,567883 |
|
1 |
||||
|
5 |
|
|
-4,735560 |
7,08503 |
-0,997213 |
0,998752 |
0,354293 |
|
1 |
||||
|
6 |
|
|
-6,935326 |
9,284796 |
-0,977649 |
0,987343 |
1,491443 |
|
-1 |
||||
|
7 |
|
|
-8,582382 |
10,93185 |
-0,912648 |
0,943383 |
3,843778 |
|
- 1 |
||||
|
8 |
|
|
-5,237779 |
9,237779 |
-0,452003 |
0,666363 |
10,33662 |
|
1 |
||||
|
9 |
|
|
11,45391 |
-8,060213 |
0,826222 |
-0,718183 |
7,80964 |
|
1 |
||||
|
10 |
|
9,806863 |
-6,413163 |
0,873814 |
-0,761579 |
5,457331 |
|
1 |
|||||
расчета по приближенному алгоритму вполне дос |
защитил докторскую диссертацию «Исследование |
|||||||||||||
таточна |
для |
поставленных задач. |
|
совмещенных электрических машин (Основы теории |
||||||||||
|
|
|
|
|
|
|
машин и устройств с вращающимся магнитным по |
|||||||
|
|
|
СПИСОК ЛИТЕРАТУРЫ |
|
лем и разнополюсными обмотками)» в Харьковском |
|||||||||
|
|
|
|
политехническом институте. Профессор |
кафедры |
|||||||||
|
|
|
|
|
|
|
||||||||
1. |
Загрядцкий В.И., Кобяков Е.Т. К анализу напряженности |
«Электрооборудование и энергосбережение» Орлов |
||||||||||||
магнитного поля некругового витка с током в однородной изо |
ского государственного технического университета |
|||||||||||||
тропной среде. — Электричество, 2002, № 3. |
|
(ОрелГТУ). |
|
|
|
|
|
|||||||
2. Нейман Л.Р., Калантаров П.Л. Теоретические основы |
|
|
|
|
|
|||||||||
Кобяков Евгений Тихонович окончил механический |
||||||||||||||
электротехники. Ч. 3. — М.;Л.: Госэнергоиздат, 1954. |
||||||||||||||
3. |
Говорков В.А. Электрические и магнитные поля. — М.: |
факультет Всесоюзного заочного института тек |
||||||||||||
Госэнергоиздат, |
1960. |
|
|
|
стильной и |
легкой |
промышленности |
в |
1965 г. В |
|||||
4. |
Зайчик М.Ю. Сборник задач и упражнений по теорети |
1995 г. защитил кандидатскую диссертацию «Анализ |
||||||||||||
ческой |
электротехнике. — М.: Энергоатомиздат, 1988. |
|||||||||||||
и синтез динамических систем в задачах проектиро |
||||||||||||||
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
[19.05.05] |
вания испытательных машин осевого |
циклического |
||||||
Авторы: |
Загрядцкий |
Владимир Иванович |
нагружения и роторов» в ОрелГТУ. Профессор ка |
окончил электротехнический факультет Горьковско-федры «Динамика и прочность машин» этого универ го политехнического института в 1954 г. В 1973 г. ситета.