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

Шаровая молния Баранов М.И

.pdf
Скачиваний:
28
Добавлен:
02.02.2015
Размер:
4.43 Mб
Скачать

Линейная ЭДС (рис. 2,а) снималась на зажимах

обмотки статора при вращении ротора со скоростью ω =200 с1 , а фазная ЭДС (рис. 2,б) определялась с

помощью пробной обмотки с числом витков N = 11 при вращении ротора со скоростью ω =110 с1 .

Из представленных осциллограмм следует, что, вследствие синусоидальности ЭДС, ВДПМ данной конструкции имеет широкие возможности примене- ния как в простых электроприводах со 120-градусной коммутацией транзисторных ключей, так и в преци- зионных электроприводах с векторным управлением.

Для анализа динамических и квазиустановивши- ся режимов работы такого двигателя используется система уравнений, записанная в координатных осях

(d;q) [2]:

 

u

 

= r i

 

+ L

did

− ωL i

;

 

 

 

 

 

 

 

d

S d

 

d

dt

q q

 

u

 

= r i

+ L

 

diq

+ ωL i

+ ωΨ ;

 

 

 

 

 

 

q

S q

 

q

dt

d d

m

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

(1)

 

MЭ =

p[i

Ψ

 

+ i i

(L

L )];

 

 

 

 

 

 

2

 

 

q

 

m

 

 

d q

d

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

JΣ

dω

= MЭ МС,

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

где rS -

активное сопротивление обмотки статора;

Ld , Lq -

продольная и

поперечная

индуктивности;

Ψm - максимальное потокосцепление обмотки статора

с потоком ротора; ω - скорость

вращения ротора;

MЭ - электромагнитный момент;

JΣ - суммарный мо-

мент инерции; MС - момент сопротивления.

Для решения системы уравнений (1) необходимо

знать индуктивности Ld , Lq и потокосцепление Ψm ,

которые могут быть получены в результате решения полевой задачи. Полевая модель является достаточно точной для этих целей, т.к. позволяет учесть реаль- ную геометрию ВДПМ, насыщение стали, характери- стики постоянных магнитов. Для нахождения индук- тивностей Ld , Lq были проведены полевые расчеты для двух положений ротора при МДС катушки обмот- ки статора Fs=0 A и Fs=66 A. При этом две фазы об- мотки статора включались последовательно, т.е. из шести катушек задействованы четыре.

На рис. 2 показано распределение магнитного поля при ориентации ротора по оси d. Сопоставление распре- деления векторного магнитного потенциала на рис. 2,а (Fs=0 A) и рис. 2,б (Fs=66 A) показывает, что реакция якоря является намагничивающей, о чем свидетельству-

ет увеличение числа магнитных силовых линий. Максимальное потокосцепление Ψm обмотки

статора с потоком ротора:

 

 

Ψm = 2 nk ( 1Fs =0 +

2Fs =0 ) =

(2)

 

,

= 2 44 (7.048 + 7.047) 105 = 0.0124 Вб

 

где nk - число витков

катушки обмотки статора,

Ф1Fs =0 , Ф2Fs =0 - магнитные потоки через зубцы статора (см. рис. 2) при МДС катушки обмотки статора Fs=0 A.

а)

б)

Рис. 2. Магнитное поле ВДПМ по оси d

при Fs=0 (а) и Fs=66 А (б)

Индуктивность Ld по продольной оси по дан- ным полевого расчета определяется следующим обра- зом:

2

nk2

(

1

1

+ 2

Fs =66

2

Fs =0

)

 

L =

 

 

 

Fs =66

Fs =0

 

 

 

, (3)

 

 

 

 

 

 

 

 

 

 

d

 

 

 

FS

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где Ф1Fs =66 , Ф2Fs =66 - магнитные потоки при Fs=66 A. На рис. 3 показано распределение магнитного поля при ориентации ротора по оси q. По отношению к оси d ротор смещен против часовой стрелки на ме-

ханический угол 22,5° (90 эл. град.).

Електротехніка і Електромеханіка. 2008. 6

41

 

 

 

 

 

Таблица 2

 

Ф1Fs =0

Ф2Fs =0

Ф1Fs =66

Ф2Fs =66

 

L

 

Вб*10-5

Вб*10-5

Вб*10-5

Вб*10-5

 

мГн

Ось d

7.048

7.047

8.459

8.458

 

1.656

Ось q

-4.084

4.092

-2.736

5.433

 

1.578

Данные табл. 2 подтверждают известный факт о том, что у двигателей с постоянными магнитами ин- дуктивности по осям d и q мало отличаются [3]. В данном случае их отличие не превышает 5%.

Для идентификации значения индуктивности обмотки статора LS Ld Lq был проведен экспери-

мент на установке, схема которой представлена на рис. 4.

а)

б)

Рис. 3. Магнитное поле ВДПМ по оси q

при Fs=0 (а) и Fs=66 А (б)

Индуктивность Lq по поперечной оси по дан- ным полевого расчета определяется следующим обра- зом:

 

 

2

nk2

( 1

1

+ 2

Fs =66

2

Fs =0

)

 

 

L

q

=

 

 

Fs =66

Fs =0

 

 

 

,

(4)

 

 

 

FS

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

Ф1Fs =0 ,

Ф2Fs =0 -

магнитные потоки (см.

рис.

3)

при МДС катушки обмотки статора Fs=0 A, Ф1Fs =66 , Ф2Fs =66 - магнитные потоки при Fs=66 A.

Результаты определения магнитных потоков и результаты расчета индуктивностей по (3) и (4) све- дены в табл. 2.

Рис. 4.

В данной схеме при подаче однополярного пита- ния +U осуществляется импульсное регулирование тока, проходящего через две последовательно вклю- ченные фазы статора (Аx, By). Ротор ВДПМ при этом автоматически ориентируется по оси d. Амплитуда изменения тока составляет i, измерение тока выпол- няется резистором Rизм.

При включении транзистора VT уравнение элек- трического равновесия для цепи [+U – Ax – yB – Rизм

– VT – общий провод]:

U = R

i

+ L

di1

,

(5)

 

1

1

S

dt

 

где R1 - активное сопротивление контура прохожде- ния тока i1 , t - время. В данном режиме работы схе-

мы di1 > 0 . dt

При выключении транзистора VT уравнение электрического равновесия для цепи [ Ax – yB – Rизм

диод VD]:

0 = R

i

+ L

di2

,

(6)

 

2

2

S

dt

 

где R2 - активное сопротивление контура прохожде-

42

Електротехніка і Електромеханіка. 2008. 6

ния тока i2 . В этом режиме di2 < 0 . dt

График изменения тока двигателя в функции времени, построенный по результатам эксперимента, приведен на рис. 5.

Рис. 5. График изменения тока ВДПМ

составляет 1,57 Ом, а сопротивление токоизмеритель- ного резистора в схеме по рис. 4 Rизм = 0,2 Ом. Учи- тывая сопротивление полупроводниковых приборов VT и VD, можно сделать вывод о достаточно хоро- шем соответствии расчетных значений сопротивлений R1 и R2 реальным значениям сопротивлений этих контуров.

На рис. 6 представлена расчетная механическая характеристика исследуемого в данной работе ВДПМ, полученная на основе решения системы уравнений (1) с учетом изложенной выше методики определения индуктивностей обмотки статора по осям Ld , Lq и потокосцепления Ψm . Также на рис. 6 представлена экспериментальная механическая характеристика ВДПМ. Сопоставление расчетной и эксперименталь- ной механических характеристик позволяет сделать заключение о том, что решение полевой задачи в двухмерной постановке дает возможность с достаточ- ной точностью определить все необходимые данные для расчета динамических режимов электропривода с явнополюсным ВДПМ.

Из уравнений (5), (6) могут быть получены сред- ние значения токов I1 и I2 за интервалы времени t1 и t2 в квазиустановившемся режиме регулиро-

вания тока:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U LS

 

i1

 

 

 

 

 

LS

 

i2

 

 

 

 

t

 

 

 

 

 

t

2

 

 

I

=

 

 

1

,

 

I

2

=

 

 

 

 

.

(7)

 

 

 

 

 

 

 

 

 

 

1

 

R1

 

 

 

 

 

 

 

 

 

R2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Исходя из того, что

 

i1 = − i2 =

 

i ,

I1 = I 2 и

считая, что R1 R2

индуктивность статора определя-

ется так:

 

 

 

 

 

 

 

 

t1

 

 

 

 

 

 

 

 

 

 

L

 

=

U

 

 

t2

.

 

 

 

(8)

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

 

i

 

t1 +

t2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω,c

 

500-1

 

450

Расчет

400

 

350

 

300

Эксперимент

250

 

200

 

150

 

100

 

50

 

0

0

0.02

0.04

0.06

0.08

M,Нм0.1

По

экспериментальным данным

U = 24,1 В,

Рис. 6. Механические характеристики ВДПМ

i = 0,43 А,

t1 = 32 мкс,

t2 = 228 мкс. Тогда индук-

 

 

тивность по (8) LS = 1,572 мГн, что близко к расчет-

ЛИТЕРАТУРА

ному значению, полученному исходя из решения по-

[1] Калюжный А.Н., Рымша В.В. Вентильные электродви-

левой задачи. Здесь же отметим, что полевая задача

гатели и электроприводы производства николаевского

завода "Электротехника"

// Рынок электротехники. –

решалась в плоской постановке, при которой вклад

2007. 4(8). – С. 90.

 

лобовых частей обмотки статора в общую индуктив-

 

[2] Столов Л.А., Афанасьев

А.Ю. Моментные двигатели

ность не учитывается. Однако, ввиду того, что в рас-

постоянного тока. – М.: Энергоатомиздат, 1989. – 224 с.

сматриваемой

конструкции ВДПМ с

явновыражен-

[3] Косулин В.Д., Михайлов Г.Б., Омельченко В.В., Путни-

ными полюсами вылет лобовых частей мал, результат

ков В.В. Вентильные электродвигатели малой мощности

расчета индуктивности в плоской постановке полевой

для промышленных роботов. – Л.: Энергоатомиздат,

задачи получается достаточно корректным.

Ленингр. Отд-ние, 1988. – 184 с.

Для

дополнительной

проверки

правильности

 

 

идентификации индуктивности обмотки статора и

 

Поступила 03.04.2008

ввода в расчет условия R1 R2 по уравнениям (7)

 

 

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

 

 

ров прохождения токов i1

и i2 . В эксперименте при

 

 

регулировании ток обмотки статора изменялся в пре-

 

 

делах от 1,14 А до 1,57 А, т.е. его среднее значение

 

 

составило I = 1,355 А.

 

 

 

 

Тогда по (7) сопротивление R1 = 2,19 Ом, сопро-

 

 

тивление

R2 =2,2 Ом. В то же время собственное со-

 

 

противление обмотки статора при температуре 20 °С

 

 

 

 

 

Електротехніка і Електромеханіка. 2008. 6

 

43

Теоретична електротехніка

УДК 621.317.4

УНИВЕРСАЛЬНЫЕ СЕЛЕКТИРУЮЩИЕ КОНТУРЫ ИЗМЕРИТЕЛЬНОЙ СИСТЕМЫ РАЗНОТИПНЫХ ПРОСТРАНСТВЕННЫХ ГАРМОНИК МАГНИТНОГО ПОЛЯ

Гетьман А.В., к.т.н.

Научно-технический центр магнетизма технических объектов Национальной академии наук Украины Украина, 61106, Харьков, ул. Индустриальная, 19

тел. (0572) 99-11-75

Досліджуються практичні аспекти просторового гармонічного аналізу з метою наближення до поверхні технічних об'єктів області, доступної для аналітичного опису зовнішнього магнітного поля. Теоретично обґрунтована просто- рова конфігурація селектуючих контурів систем виміру сферичних і сфероїдальних гармонік магнітного поля. Уза- гальнено методику застосування контурних вимірювальних систем для визначення амплітудних коефіцієнтів різно- типних гармонік скалярного магнітного потенціалу.

Исследуются практические аспекты пространственного гармонического анализа с целью приближения к поверхно- сти технических объектов области, доступной для аналитического описания внешнего магнитного поля. Теоретиче- ски обоснована пространственная конфигурация селектирующих контуров систем измерения сферических и сферои- дальных гармоник магнитного поля. Обобщена методика применения контурных измерительных систем для опреде- ления амплитудных коэффициентов разнотипных гармоник скалярного магнитного потенциала.

ВВЕДЕНИЕ.

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

Однако практическое проведение ПГА МП вбли- зи поверхности ТО возможно лишь при знании ам- плитуд ряда пространственных гармоник, позволяю- щего с необходимой точностью восстанавливать маг- нитное поле во всей области практического использо- вания модели. В свою очередь измерение достаточно- го количества и с необходимой точностью амплитуд пространственных гармоник связано с затруднениями технического и методического характера [2]. Тем са- мым методическая погрешность выделения вклада каждой гармоники в суммарный скалярный потенци- ал, ограниченность количества измеряемых гармоник и технологические сложности эксплуатации сущест- вующих измерительных систем являются сдержи- вающим фактором на пути применения ПГА МП вблизи поверхности ТО.

Как известно [3], разработанные ранее динами- ческие контурные измерительные системы, для прак- тического определения полного набора мультиполь- ных коэффициентов построены на интегральном ме- тоде измерения. Измеряемой величиной в них являет- ся магнитный поток, пронизывающий контурные об- мотки специальной конфигурации [4]. Благодаря спе- циальной форме контуров, поток в них оказывается состоящим из суммы вкладов мультиполей только одного порядка.

Однако, как поиск пространственной конфигура- ции селектирующих контуров, так и разработка мето- дики обработки сигнатур изначально проводилось для

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

С другой стороны обобщенный метод ПГА [5] изначально предопределяет практическую возмож- ность его применения к ТО с различными соотноше- ниями их габаритных размеров на единой метрологи- ческой базе. В связи с чем, очевидна необходимость использования теоретических основ обобщенного метода ПГА для построения аналогичной динамиче- ской контурной системы, позволяющей проводить селекцию и измерение как сферических, так и сфе- роидальных гармоник магнитного поля. Для чего не- обходимо провести поиск пространственной конфи- гурации контуров, селектирующих разнотипные про- странственные гармоники заданного порядка.

Целью работы является синтез универсальных селектирующих контуров измерительной системы сферических и сфероидальных гармоник МП ТО.

ПОСТАНОВКА ЗАДАЧИ Обобщенный метод ПГА использует общее

представление базисных решений уравнения Лапласа - для вытянуто-сфероидальных, сферических и сплюснуто-сфероидальных гармоник [6] в виде:

n

σ Fnm (x1 )

σPnm (x2 )

σ

m

cos mϕ

 

 

U = ∑ ∑

 

cn

 

, (1)

n=1m=0

 

 

σ sm sin mϕ

 

 

 

 

 

n

 

 

 

где x1 и x2 соответствующие "радиальная" и "угло- вая" координаты;

n и m – соответственно степень и порядок простран-

ственной гармоники;

σFnm(x1) и σGnm(x2) – "радиальная" и "угловая" функ- циональные зависимости пространственных гармоник

[6], тип которых определяется параметром σ, так что при σ=1 они вытянутосфероидальные, при σ=2 – сфе- рические, а при σ=3 – сплюснутосфероидальные;

44

Електротехніка і Електромеханіка. 2008. 6

cnm , snm

амплитуды соответствующих пространст-

Воспользовавшись

геометрическим

 

условием

венных гармоник скалярного потенциала магнитного

для представления элемента длины контура и ортого-

нальностью тригонометрических

функций,

 

получим

поля.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

окончательное выражение для потока в виде

 

 

Будим исходить из того, что величина Ф являет-

 

 

 

 

h

 

U 0(x , x )

 

 

 

 

тирующим контуром, находящимся в воздухе с отно-

 

 

dx1 .

(6)

Ф = 2πμ0R0 h

n

x

1

 

2

ся результатом измерения МП, сцепленного с селек-

 

 

1

 

 

 

 

 

 

 

сительной проницаемостью μ=1, и подлежит даль-

 

 

2

 

 

 

2

 

 

 

 

 

 

Таким образом, мы показали в (6), что круговой

нейшей обработке для определения амплитуд про-

контур аксиально-симметричный аппликате пред-

странственных гармоник. При этом ключевым момен-

ставляет собой аппаратный фильтр, выделяющий

том разработки селектирущих контуров динамиче-

суммарный вклад зональных (m=0) гармоник в прони-

ских измерительных систем является использование

зывающем его магнитном потоке. Причем такая се-

представления пространственной гармоники скаляр-

лекция одинаково характерна для вытянуто-

ного потенциала магнитного поля через его вектор-

сфероидального,

сферического

 

 

и

 

 

сплюснуто-

ный потенциал.

 

 

 

 

 

 

 

 

 

 

сфероидального базисных решений.

 

 

 

 

 

 

 

 

 

 

 

 

 

h

h1h2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A = −μ ×

 

ϕ

Udx e

,

(2)

СЕЛЕКЦИЯ ТЕССЕРАЛЬНЫХ ГАРМОНИК

 

 

 

 

 

 

 

 

0

 

 

h

h

1

x1

 

Проведем поиск возможных вариантов специ-

 

 

 

 

 

 

 

2

 

ϕ

 

 

 

где h1,

h2 ,

hϕ коэффициенты Ламе для соответ-

альной конфигурации селектирующих контуров для

использования их

в качестве аппаратного

 

фильтра

ственно "радиальной",

 

"угловой" и циклической ко-

 

 

тессеральных гармоник порядка m’, основываясь на

ординат;

erx1

единичный вектор в направлении ко-

следующих положениях.

 

 

 

 

 

 

 

 

 

 

 

ординаты x1.

 

 

 

 

 

 

 

 

 

 

 

 

 

Поскольку определение зональных амплитудных

Такое представление (2) позволяет проводить

коэффициентов проводится при прохождении источ-

синтез пространственной конфигурации селектирую-

ника поля через круговой контур по его оси симмет-

щих контуров, если воспользоваться выражением

рии, то очевидна необходимость использования того

магнитного потока через векторный потенциал и эле-

же перемещения и при определении остальных (тес-

 

 

 

 

r

 

 

 

 

 

 

 

 

 

серальных) амплитудных

 

коэффициентов.

 

Поэтому

мент длины контура dl

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

селектирующий контур тессеральной гармоники дол-

 

 

 

 

 

 

 

 

= Adl .

 

 

(3)

 

 

 

 

 

 

 

 

 

 

жен располагаться вместе с круговым контуром на

 

 

 

 

 

 

 

 

 

 

L

 

 

 

цилиндре, ось которого совпадает с трассой переме-

При этом именно специальная пространственная

щения. Для выделения гармоник, имеющих порядок

конфигурация контура L интегрирования в (3) "пре-

m’, из магнитного потока Ф, представленного через

вращает" его в аппаратный фильтр пространственных

векторный потенциал в выражении (3) воспользуемся

гармоник заданного порядка.

 

 

 

 

 

ортогональностью

тригонометрических

 

функций,

СЕЛЕКЦИЯ ЗОНАЛЬНЫХ ГАРМОНИК

 

входящих в векторный потенциал в (2). Поскольку

 

функциональные зависимости векторного потенциала

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Очевидно, что наиболее простую форму инте-

от координаты φ для сферического и сфероидальных

грал (3) примет, если проводить интегрирование

случаев идентичны, то такой контур должен выделять

только по циклической координате. Для всех трех

как сферические, так и сфероидальные амплитудные

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

коэффициенты порядка m’.

 

 

 

 

 

 

 

 

 

 

координата φ. В самом простом случае геометриче-

 

 

 

 

 

 

 

 

 

 

Однако применение

свойства

 

ортогональности

ской конфигурации мы имеем круговой контур ра-

для тригонометрических функций, входящих в век-

диуса R0 , параллельный плоскости XOY, с центром

торный потенциал (3) возможно, если линия контура,

на оси аппликат. При этом элемент длины контура dl

во-первых, ортогональна ортам радиальной и угловой

является функцией только координаты φ. Представим

координат, а во-вторых,

 

имеет

косинусоидальную

φ-тую составляющую векторного потенциала (един-

(синусоидальную) функциональную зависимость от

ственно дающую вклад в магнитный поток) в виде

 

циклической координаты. Кроме того, такой контур

 

 

Aϕ = μ0

h1

U (x1, x2 ,ϕ)dx1 =

 

не должен содержать трассу перемещения - ось ап-

 

 

 

 

пликат. Иначе в сцепленном с ним магнитном потоке

 

 

 

 

 

h2

 

 

 

x2

 

 

(4)

 

 

 

h1

 

U(x1, x2 )

U (ϕ)

появится вклад гармоник с осевой симметрией зо-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= μ0

 

 

 

x

dx1

 

нальных. Трудность совмещения предыдущих и по-

 

 

h

 

 

 

 

 

2

 

 

 

 

 

 

2

 

 

 

 

следнего условия преодолевается, если использовать

Из вида скалярного магнитного потенциала,

вместо одного два осе симметричных, но противопо-

представленного базисными решениями (1), следует

ложно коммутированных контура. Потребовав, чтобы

одинаковость их функциональной зависимости от

плоскость суммарного контура, образованного двумя

циклической координаты φ. Поэтому (3) можно пере-

функциональными от φ круговыми линиями, была

писать как

 

 

 

U (x1, x2 )

 

 

 

параллельна ап-пликате. Другими словами чтобы ра-

 

 

h1

 

 

 

диусы кругов были одинаковыми, а образовавшейся

 

Ф = μ0 h2

 

 

 

x2

dx1 U(ϕ)dlϕ .

(5)

суммарный контур лежал на цилиндре, осью которого

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

является трасса перемещения ТО. Тогда суммарный

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

контур будет описываться

разностью функциональ-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Електротехніка і Електромеханіка. 2008. 6

 

 

 

 

 

 

 

 

 

 

 

 

 

45

ных зависимостей линий контуров, лежащих на ци-

го контура, так и достаточного количества предлагае-

линдре. Выберем эти функциональные зависимости

мых контуров для практического определения необ-

L1(φ) и L2(φ) такими, чтобы выделялся бы вклад гар-

ходимого количества тессеральных пространственных

моник порядка m’, например, как на рис. 1. Т.е. чтобы

гармоник. Для этого необходимо расположить конту-

L1(φ) и L2(φ) задавали отклонение вдоль цилиндра от

ра соосно, как показано на рис. 1, и учесть смещение

круговой линии симметрии двух контуров в зависи-

между контурами при измерении потока. Преимуще-

мости от циклической координаты и находились бы в

ство использования таких контуров состоит в общно-

противофазе друг к другу.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

сти всей метрологической базы необходимой для

 

 

 

 

 

 

m=0

 

 

 

 

 

 

m=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m=2

 

 

 

 

 

 

практического определения коэффициентов как сфе-

R0

 

 

 

 

 

 

 

 

L2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

рических, так и сфероидальных гармоник.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L1

 

 

 

 

L2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L1

 

 

 

 

 

 

 

ВЫВОДЫ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

Таким образом, в работе теоретически обоснова-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

на

пространственная конфигурация

универсальных

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

контуров измерительной системы, селектирующих

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

сферические и сфероидальные гармоники заданного

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

порядка. Создание таких систем необходимо для

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

практического проведения пространственного гармо-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L1

 

 

нического анализа магнитного поля вблизи техниче-

 

 

 

 

 

 

 

 

 

 

 

 

L1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ских объектов инвариантно к особенностям их гео-

 

Рис. 1. Расположение селектирующих контуров на

 

 

 

 

общей цилиндрической поверхности

 

 

 

 

 

 

метрии.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда косинусоидальный (верхняя часть фигур-

 

Преимущество такой контурной системы заклю-

 

чается в возможности определения достаточного ко-

ных скобок) и синусоидальный (нижняя часть фигур-

личества гармоник, необходимого для проведения

ных скобок) контура порядка m’ запишем для кратко-

описания с заданной точностью магнитного поля тех-

сти в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

нического объекта.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

Принципиальным отличием геометрической кон-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cosm

ϕ

 

 

 

 

 

(7)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L1(ϕ) = − L2 (ϕ) = α

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

фигурации

предлагаемых

измерительных контуров

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin m,ϕ

 

 

 

 

 

 

является их более рациональное размещение на цилин-

где α - постоянный коэффициент, определяющий ам-

дрической поверхности по сравнению с контурами из-

плитуду отклонения

вдоль цилиндра,

 

 

выбираемый (в

мерительных устройств мультиполей [4], имеющих

отличие от известных из [4])

 

одного порядка с его

жёсткое ограничение на протяженность вдоль цилинд-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ра (α<0,1R0). Это существенно увеличивает полезную

радиусом R0. Тогда под в

 

 

dl

 

 

 

(3)

 

 

следует понимать

выражение:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

площадь контуров и как следствие улучшает соотно-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

шение измеренного сигнала к помехе.

 

 

 

 

 

 

 

 

 

 

 

dl = dlx

 

 

ex

 

+ dlx

 

 

 

ex

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ЛИТЕРАТУРА

 

 

 

 

 

 

 

 

 

 

 

= α

 

 

 

 

 

 

 

 

 

 

d[L1(ϕ) L2 (ϕ)] ex +

(8)

[1] Розов

В.Ю.

Системи

автоматичної компенсації

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h1

 

x1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

зовнішнього

 

магнітного поля

енергонасичених

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

обєктів: Автореф. дис... докт. техн. наук: 05.09.03. /

 

 

 

 

 

 

 

 

 

 

+ α

 

 

 

 

 

 

 

 

 

 

 

 

d[L (ϕ) L (ϕ)] e

x2

 

 

ІЕД НАНУ. – Київ, 2002. – 40 с.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

[2]

Гетьман А

В

.

Особенности исследования магнитного

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

Магнитный поток Фm’, сцепленный с контуром

 

поля технического объекта вблизи его поверхности//

порядка m’ согласно (3) запишется как

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Техническая электродинамика. – 2006.- Тематический

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

выпуск. Ч. 6. Проблеми сучасної електротехніки. -С.9-

 

 

 

 

 

 

 

 

m' = A dl = (A)x1 dlx1

+ (A)x2 dlx2

(9)

 

 

 

 

 

 

 

 

 

[3]

12.

 

 

 

 

 

 

и после упрощения выражения примет вид

 

 

 

 

 

 

Волохов С.А. Кильдишев А.В. Нормализованные маг-

 

 

 

 

 

 

 

нитные сигнатуры мультипольного источника, движу-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m'

(x1, x2 )

 

щегося по оси кругового контура // Электри-чество.

 

 

 

 

m'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U n

 

 

1997. - 3. - С. 65-67.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

− αμ0hϕ dx1 x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Φ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(10)

[4] Волохов С.А., Кильдишев А.В. Измерительные конту-

 

 

 

 

 

 

 

 

h1h2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U nm' (x1, x2 )

 

 

ры и селектирующие функции для определения тессе-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx 2 x1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ральных мультиполей интегральным преобразованием

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вид выражения

(10)

 

 

 

показывает

 

практическую

 

магнитных сигнатур // Космічна наука і технологія. –

 

 

 

 

 

 

1996. - Т.2, 5-6. - С.26-30.

 

возможность селекции из магнитного потока тессе-

 

 

[5]

Гетьман А.В. Структура контурных динамических

ральных гармоник порядка m’, с помощью контура

 

систем для практического гармонического анализа

специальной конфигурации. При этом пространст-

 

магнитного поля технических объектов // Техническая

венная конфигурация контура определяет порядок

 

электродинамика. – 2008.- Тематический выпуск. Ч. 3.

выделяемых гармоник. Тем самым контур превраща-

 

Проблеми сучасної електротехніки. -С. 97-100.

ется в аппаратный фильтр селекции гармоник задан-

[6]

Ерофенко В.Т. Теоремы сложения. Справочник. -

ного порядка. Конструктивно такой контур близок к

 

Минск: Наука и техника, 1962.-256 с

 

 

 

 

 

 

 

 

 

контуру для выделения зональных гармоник, что по-

 

 

 

 

 

 

Поступила 05.09.08

зволяет одновременное использование, как зонально-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

46

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Електротехніка і Електромеханіка. 2008. 6

УДК 621.3.011

ПОСТРОЕНИЕ И ИССЛЕДОВАНИЕ ХАРАКТЕРИСТИК ФЕРРОРЕЗОНАНСНОЙ ЦЕПИ НА ОСНОВЕ МЕТОДА СХЕМНОГО МОДЕЛИРОВАНИЯ

Канов Л.Н., к.т.н., доц.

Севастопольский национальный технический университет Украина, 99053, Севастополь, Стрелецкая бухта, СевНТУ, кафедра "Судовые и промышленные электромеханические системы"

тел. (0692) 235-272

Запропоновано схемні моделі основних лінійних і нелінійних електротехнічних елементів у стаціонарному режимі зминного струму, на підставі яких побудована схемна модель ферорезонансного кола. Запропоновано методику побу- дови вольтамперних і амплітудочастотних характеристик коли, а також визначення границь областей стійких режимів.

Предложены схемные модели основных линейных и нелинейных электротехнических элементов в стационарном ре- жиме переменного тока, на основании которых получена схемная модель феррорезонансной цепи. Предложена мето- дика построения вольтамперных и амплитудочастотных характеристик цепи и определения границ устойчивых режимов.

ВВЕДЕНИЕ

Вэлектроприводах переменного тока, где при- меняются магнитные усилители и все чаще частотно регулируемые асинхронные двигатели, важное значе- ние имеет анализ вольтамперных (ВАХ) и амплиту- дочастотных (АЧХ) характеристик феррорезонансных цепей [1]. Точное аналитическое исследование перио- дических режимов в таких цепях невозможно. Суще- ствует несколько методов приближенного аналитиче- ского определения параметров периодических режи- мов, постоянно появляются новые их модификации. Так в [2] на основе метода оптимальной линеаризации изложен итеративный подход к расчету установив- шихся режимов нелинейных цепей переменного тока, особенностью которого является возможность нахож- дения параметров нелинейных элементов по их ВАХ.

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

Вряду этих методов выделяется простой и эко- номный метод гармонической линеаризации [4, 5], который является основой частотного анализа нели- нейных цепей. Этот метод предполагает близкую к синусоидальной форму токов и напряжений в цепи и позволяет использовать комплексные амплитуды. Существует несколько модификаций этого метода, позволяющих повысить его точность, например, [6]. Обзор этих методов показывает, что все они связаны с большим количеством аналитических преобразований и не позволяют получить явных выражений для ха- рактеристик феррорезонансных цепей. Графические методы построения таких характеристик громоздки и имеют невысокую точность.

ПОСТАНОВКА ЗАДАЧИ Целью статьи является разработка методики

применения схемного моделирования для построения ВАХ и АЧХ нелинейных цепей в рамках условий ме- тода гармонической линеаризации, т.е. в предположе- нии малого отличия форм токов и напряжений от си- нусоидальной. Схемное моделирование представляет эффективный численный метод исследования режи- мов переменного тока в нелинейных цепях [7].

МАТЕРИАЛЫ ИССЛЕДОВАНИЯ Определим схемные модели основных элементов

в режиме синусоидального тока. Для линейного ком-

плексного

сопротивления

Z = R + jX

справедливо

(R + jX )(I′ + jI′′) = U′ + jU ′′ ,

где R, X

активное и

реактивное

сопротивление;

I , I ′′,U,U′′

веществен-

ные и мнимые части комплексных, действующих то-

ков и напряжений. Отделяя вещественные и мнимые части, получаем RI ′ − XI ′′ = U; RI ′′ + XI ′ = U′′ . Этим

уравнениям соответствует схемная модель в виде па- ры последовательных цепей, изображенная на рис. 1, на которой реактивное сопротивление входит коэф- фициентом управления в управляемые источники на- пряжения. Очевидно, в схемной модели резистора эти источники будут отсутствовать, а для реактивных элементов будет отсутствовать активное сопротивле- ние. Аналогично строятся схемные модели проводи- мостей, взаимных индуктивностей [7] и др.

Рис. 1. Схемная модель комплексного сопротивления

Нелинейные элементы с однозначной ВАХ по

действующим значениям токов и напряжений I = f (U ) представляются парой нелинейных управ-

ляемых источников тока, соответствующих вещест- венным и мнимым частям тока в элементе. Эти ис-

Електротехніка і Електромеханіка. 2008. 6

47

точники зависят от вещественных и мнимых частей

Магнитный поток кроме основной гармоники содер-

напряжения на них: I ′ = f1(U ,U ′′);

 

I ′′ =

f2 (U ,U ′′) .

жит постоянную составляющую Ф0 :

 

 

 

 

 

 

 

 

 

 

 

 

Для конкретизации этих функций обратимся к записи

 

 

 

 

 

=

0 +

 

Мcos ωt , Ф

 

 

 

 

 

=

U M

,

 

 

 

 

 

(6)

активной

и

реактивной

мощностей

элемента:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

М

 

 

 

 

 

 

 

 

 

 

 

 

P = U I ′ + U ′′I ′′ ; Q = U ′′I ′ − U I ′′ . В резистивном нели-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ωw

 

 

 

 

 

 

 

 

 

 

где U M амплитуда напряжения на рабочей обмотке.

нейном элементе

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Выделим в выражении веберамперной характе-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U ′′I ′ = U I ′′ ,

 

 

 

 

 

 

 

 

 

 

 

 

 

(1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ристики

слагаемые,

 

 

соответствующие

 

 

постоянным

в реактивном элементе

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

составляющим тока и потока

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U I ′ = −U ′′I ′′ .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

+ 3a

 

 

 

 

 

 

 

 

 

+ а Ф3

 

 

 

 

 

 

 

 

 

 

 

Запишем

ВАХ

нелинейного

 

элемента в

виде

w

 

I

 

 

a

 

 

 

 

Ф

 

 

 

 

 

 

 

(7)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I

2

=

f

2

(U ) =

(I )

2

+

 

(I ′′)

2

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

y

 

1

 

 

 

 

3

ωw

 

 

0

 

 

 

3 0

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда

для резистивного

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

элемента с учетом (1) получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и составляющим первой гармоники

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I′ = ±

 

U

f (U ); I ′′ = ±

U ′′

f (U ) ;

 

 

 

 

(3)

 

 

U

M

 

3Ф2

U

M

 

 

 

 

3

 

 

 

 

 

U

M

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

wi = − a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U

U

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

a

 

 

 

 

 

 

cos ωt .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 ωw

 

 

 

 

 

 

0

 

ωw

3

 

 

 

 

 

3

ωw

 

 

 

для реактивного элемента с учетом (2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I′ = ± U′′ f (U ); I′′ = ± U f (U ) .

 

 

 

 

(4)

ной индуктивности по действующим значениям тока

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Последнее выражение дает уравнение ВАХ нелиней-

 

 

 

 

 

 

 

 

 

 

 

U

 

 

 

 

 

 

 

 

 

U

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

Из выражения активной мощности следует для рези-

 

 

 

 

 

 

 

 

 

 

 

 

U

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

U

 

 

 

 

 

 

 

I =

 

 

 

a + 3a Ф

 

+ 1,5a

 

 

 

 

, в

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и напряжения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U I ′ > −U ′′I ′′ .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

3

 

 

 

стивного

элемента

Подставляя

сюда

 

 

 

 

 

 

 

 

 

 

 

ωw

2

1

 

 

 

 

3

 

 

 

 

 

 

 

 

 

ωw

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U

(I ′′)2 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U′′,U

 

 

из

 

(1),

 

получаем

 

 

 

U I ′ > −

котором

 

 

 

 

 

определяется из (7)

 

 

 

и зависит от тока

 

 

 

 

 

 

 

I

Ф0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U ′′

(I )2 > −U ′′I ′′ . Поэтому знаки U

и I , а также U ′′

подмагничивания. В соответствии с этим выражения

 

(5) для управляемых источников тока схемной модели

 

 

 

I′′

в выражениях (3) совпадают,

 

 

 

 

 

 

 

индуктивности принимают вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и

 

 

I ′′

и эти выражения

 

 

 

 

 

U ′′

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U

2

 

 

 

 

принимают окончательный вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U

 

 

 

 

U ′′

 

 

 

 

 

 

 

 

 

I

=

 

 

 

 

 

 

 

a

 

+

3a Ф2 + 1,5a

3

 

 

 

 

 

 

;

 

 

(8)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I′ =

 

 

 

 

I ′′ =

 

 

 

 

 

 

 

 

 

 

 

 

 

ωw

2 1

 

 

 

 

3 0

 

 

 

 

 

 

 

 

 

 

ωw

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (U ) ;

 

 

 

 

 

f (U ) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U

 

U

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I ′′ = −

I

U

; U 2 = (U )2 + (U ′′)2 .

 

 

 

 

 

 

 

 

 

 

 

Из выражения для реактивной мощности в ин-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дуктивном элементе U ′′I ′ > U I ′′ .

Подставляя

сюда

 

 

 

 

 

 

 

 

U ′′

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U ,U ′′

 

 

из

(2),

 

получаем

 

 

U ′′I ′ > −

U ′′

 

(I ′′)2 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U

(I )2 > U

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I ′′ . Поэтому знаки U ′′

и I

в выраже-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I ′′

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ниях (4) совпадают, а знаки U и

I ′′

противополож-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ны, и эти выражения принимают окончательный вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I ′ =

U ′′

 

f (U ) ;

I′′ = −

U

 

f (U ) .

 

 

 

 

(5)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U

 

 

 

 

 

 

U

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для нелинейного емкостного элемента знаки в (5) будут противоположны. Аналогичное представление

можно выполнить для нелинейных элементов с ВАХ вида U = f (I ) .

С помощью полученных схемных моделей эле- ментов строятся схемные модели цепей переменного тока, представляющие собой две цепи, объединенные управляемыми источниками. Расчет этих моделей дает вещественные и мнимые составляющие токов и напряжений в исследуемой цепи.

Рассмотрим построение характеристик последо- вательной феррорезонансной цепи с подмагничивани- ем. Примем для веберамперной характеристики ин- дуктивности полиномиальную аппроксимацию F = a1Ф+ a3Ф3 , где магнитодвижущая сила опреде-

ляется как током i рабочей обмотки

w , так и током

I y обмотки подмагничивания wy :

F = wi + wy I y .

Рис. 2. Схемная модель последовательной феррорезонансной цепи

Схемная модель последовательной феррорезо- нансной цепи изображена на рис. 2, где в соответст-

вии с (7) обозначены: g1 = a

+ 3a

 

U

2

управляе-

 

 

 

 

 

 

1

 

3

ωw

 

мая проводимость, g2 = a3Ф02 нелинейная проводи- мость, J 3 = wy I y , r сопротивление рабочей обмот-

ки, управляемые источники напряжения представля-

48

Електротехніка і Електромеханіка. 2008. 6

[4, 5, 8],

ют емкость, а сопротивление rc имитирует потери на

где ФМ (t), ϕ(t) – медленно меняющиеся амплитуда

гистерезис и принимается постоянным [8]. Управляе-

первой гармоники потока и его начальная фаза.

мые источники тока J1, J 2 определяются в соответст-

 

вии с (8). Напряжение на проводимости g1 численно

 

соответствует постоянной составляющей магнитного

 

потока.

 

Рис. 3. Вольтамперная характеристика последовательной феррорезонансной цепи

На рис. 3 изображена ВАХ цепи, полученная просчетом схемной модели при изменении величины J BX . Входное напряжение подсчитывается по выра-

жению U

BX

=

(U

)2

+ (U′′

)2 . Параметры цепи в

 

 

BX

 

BX

 

 

относительных

единицах:

r = 0,025 ;

ωw = 1 ;

xc = 0,22 ;

rc = 100 ;

a1

= 0,5 ;

a3 = 1; wy I y = 1 . Про-

вал напряжения при I = 8 возникает из-за насыщения

индуктивности при больших токах, второй провал при I = 3 объясняется искривлением ВАХ индуктивности

вследствие подмагничивания. При плавном измене- нии входного напряжения в цепи возникает двойной триггерный эффект.

Для построения АЧХ в схемной модели на рис. 2,а источник тока J BX следует заменить единич- ным источником напряжения, а на рис.2,б контакты с левой стороны схемы закоротить. На рис. 4 изображе- на АЧХ цепи, построенная просчетом измененной таким образом схемной модели при изменении часто- ты. Наклон резонансной кривой вправо (участок 2 – 4) при высоких токах объясняется повышением резо- нансной частоты с падением индуктивности, а наклон влево (участок 1 – 2) при небольших токах возникает при некотором возрастании индуктивности вследст- вие подмагничивания.

АНАЛИЗ РЕЗУЛЬТАТОВ Для оценки устойчивости режима цепи обратим-

ся к уравнению

1

uL + ri(Ф, I y ) + C i(Ф, I y )dt = U M sin ωt ,

где i(Ф, I y ) – нелинейная зависимость между мгно- венными значениями токов и магнитного потока. По- сле дифференцирования получаем

 

duL

= f (Ф, I y ) + ωU M cos ωt ,

(9)

 

 

 

dt

i

 

 

 

 

 

 

где f ( , I y ) = −r

 

d

1

i(Ф, I y ) . Перепишем

 

dt

 

 

 

 

 

 

C

 

далее выражение (6) в более общем виде:

 

(t) = 0 +

М (t) cos x ; x = ωt + ϕ(t) ,

 

Рис. 4. Амплитудочастотная характеристика последовательной феррорезонансной цепи

Дифференцируя последнее выражение по "ко- роткому времени", как это принято в методе медленно меняющихся амплитуд получаем

 

dФ

= −ωФM sin x . Дифференцирование же по "длин-

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

ному времени" дает

 

 

 

 

 

 

 

 

 

 

d

 

 

 

dϕ

 

 

 

 

 

 

 

 

cos x Ф

(t)

ω +

sin x = −ωФ sin x . (10)

 

 

 

 

 

 

dt

 

 

 

dt

 

 

 

С учетом равенства

uL (t) = w

dФ(t)

 

и соотношения

dt

 

 

 

 

 

 

 

 

 

 

 

(10) уравнение (9) принимает вид

 

 

 

 

dФ

sin x + Ф

ω +

dϕ

cos x

1

 

f (Ф, I y ) =

 

 

 

ωw

 

dt

 

 

dt

 

= − U M cos ωt . (11) w

Уравнения (10), (11) представляют уравнения для медленно меняющихся амплитуды магнитного потока (t) и начальной фазы ϕ(t) . После приведения

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

период по x ,

получаем укороченные уравнения, пра-

вые части которых не зависят явно от времени:

 

 

dФM (t)

= F1 (Ф

, ϕ);

dϕ(t)

= F2 (Ф

, ϕ) .

(12)

 

 

dt

dt

 

 

 

 

 

 

 

 

 

 

Необходимое условие устойчивости стационар-

ного режима состоит в выполнении неравенств

 

 

 

(a11 + a22 ) > 0 ; a11a22 a12a21 > 0 ,

 

где a =

F1

 

; a

= F1 ; a

21

=

F2

; a

22

= F2 .

 

 

 

11

 

12

∂ϕ

 

 

 

∂ϕ

 

 

 

 

 

 

 

 

 

 

Непосредственная проверка этих условий для иссле- дуемой цепи достаточно громоздка, однако для второ- го условия непосредственно по виду АЧХ можно оп-

ределить граничные по устойчивости точки, в кото-

рых а11а22 а12а21 = 0 .

Для этого отметим, что условия стационарного

режима

укороченных

уравнений имеют вид:

F1(

, ϕ) = 0 ; F2 (

, ϕ) = 0 . Последнее равенство

задает неявную функцию ϕ = ϕ( ) ; тогда для ста-

Електротехніка і Електромеханіка. 2008. 6

49

ционарного режима получаем уравнение относитель- но одной переменной

 

 

F1(

, ϕ(

)) = 0 .

 

(13)

Согласно

правилам

дифференцирования

 

неявной

 

 

 

 

 

F1

F1

 

 

F2

 

 

функции получаем

 

dF1

=

 

Ф

 

. Следо-

 

 

Ф

 

 

 

 

 

dФ

 

∂ϕ

F2

 

 

 

 

 

 

 

 

 

 

 

∂ϕ

 

 

вательно,

граничное

условие состоит

в

равенстве

dF1 = 0 . Поэтому для выявления граничных точек d

на АЧХ, которая является графическим представле- нием уравнения (13), нужно найти точки касания АЧХ с вертикальными прямыми, так как при задан- ном токе управления амплитуды рабочего тока и по- тока связаны однозначно.

В соответствии с выполненным анализом на рис.4 выделим граничные точки 1 – 4, отделяющие устойчивые и неустойчивые режимы. Из физических соображений можно заключить, что режимы на уча- стках 0 – 1 и 4 – 5 устойчивы. Тогда участкам 1 – 2 и 3 – 4 соответствуют неустойчивые режимы. К устой- чивым можно отнести участок 2 – 3. В соответствии с этим в цепи возможны множественные триггерные эффекты при плавном изменении частоты.

При плавном увеличении частоты с точки 1 про- исходит скачок тока вверх на участок 2 – 3; при даль- нейшем увеличении частоты с точки 3 на вершине кривой происходит скачок тока вниз на участок 4 – 5. При плавном уменьшении частоты с точки 4 происхо- дит скачок тока вверх на участок 2 – 3; при дальней- шем уменьшении частоты происходит скачок тока вниз на участок 0 – 1. Кроме того, при увеличении частоты с точки 1 скачок тока может происходить сразу на участок 4 – 5, минуя вершину графика. По- этому неустойчивые участки не могут быть достигну- ты при изменении частоты.

ВЫВОДЫ Предложены схемные модели основных линей-

ных и нелинейных электротехнических элементов в стационарном режиме переменного тока в условиях метода гармонической линеаризации. Построена схемная модель последовательной феррорезонансной цепи, состоящая из сопротивлений, линейных, нели- нейных, управляемых и независимых источников на- пряжения и тока.

Предложена методика построения вольтампер- ных и амплитудочастотных характеристик ферроре- зонансных цепей с помощью полученных схемных моделей, основанная на просчете режима при измене- нии входного тока и частоты.

Обосновано применение метода медленно ме- няющихся синусоид для выявления граничных точек на АЧХ, отделяющих участки, соответствующие ус- тойчивым и неустойчивым режимам феррорезонанс- ной цепи. Граничные точки являются точками каса- ния вертикальных прямых и графика амплитудоча- стотной характеристики.

ЛИТЕРАТУРА

[1]Соколовский Г.Г. Электроприводы переменного тока с частотным регулированием. – М.: Изд. центр "Акаде-

мия", 2006. – 272 с.

[2]Калугин Е.Н. Метод оптимальной линеаризации для расчета установившихся режимов нелинейной электри- ческой цепи // Электричество. – 1989. – 10. – С. 53-60.

[3]Орешников В.Г. Метод анализа вынужденных колеба- ний в нелинейных цепях // Известия ВУЗов. Электроме-

ханика. – 1997. – 6. – С. 9-11.

[4]Бессонов Л.А. Теоретические основы электротехники. Электрические цепи. – М.: Изд-во "Гардарики", 2006.– 701 с.

[5]Данилов Л.В Теория нелинейных электрических цепей / Л.В.Данилов, П.Н.Матханов, Е.С.Филиппов. – Л.: Энер- гоатомиздат, 1990. – 256 с.

[6]Бирюк Н.Д. Анализ колебаний в нелинейном контуре методом комплексных амплитуд / Н.Д.Бирюк, В.Н.Дальчев // Электричество. – 1988. – 8. – С. 46-51.

[7]Канов Л.Н. Схемное моделирование нелинейных элек- трических цепей переменного тока // Вестник СевГТУ. Вып. 41. Информатика, электроника, связь: Сб. науч. тр. Севастоп. нац. техн. ун-т.Севастополь: Изд-во Сев-

НТУ, 2002. – С. 151-154.

[8]Филиппов Е. Нелинейная электротехника. – М.: Энер-

гия, 1976. – 496 с.

Поступила 30.05.2008

50

Електротехніка і Електромеханіка. 2008. 6