Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Расчет напряженности магнитного поля плоского витка с током (80

..pdf
Скачиваний:
3
Добавлен:
15.11.2022
Размер:
174.28 Кб
Скачать

4 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):

нщ--^Шг',+Л,г'

(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

 

<

 

Соответствующее значение, найденное по фор­ мулам [1] с использованием эллиптических инте­

гралов: Н ^ = — 2,11776. Сравнивая результаты,

И—Я

получаем:

м - ^ 1 0 0 % 2=0,674%.

 

Н М

Определенный интерес представляет определе­ ние напряженности магнитного поля в точке /V

(рис. 4), расположенной за пределами контура витка с током. Результаты вычислений представ­ лены в табл. 3 и 4.

Напряженности магнитного поля в точке JV ОТ тока / каждого элемента./ контура витка находим по формуле (3) с учетом значения коэффициента £ •:

HANJB{

= /-0,0399027;

НNА2В2

_

i 0,143534;

 

 

 

 

 

# £ * з =--/-0,755212;

Я АУА

=/-0,0881304;

/V

N

 

 

 

 

Н*&=±-

0,00434386;

Н ^6 *6

= - -/- 0,00649975;

N

'

N

 

 

N

 

HpBl

 

 

• AQBo

 

 

= --^0,00799604; H'tf"* = -^0,02073792;

HAJB<) =-/-0,0138341;

И ^о*н> =-/-0,0205659.

/v

 

"

 

Результирующая напряженность поля в точке 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 г. ситета.

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]