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

книги из ГПНТБ / Кирпатовский, С. И. Периодические процессы в нелинейных цепях учеб. пособие

.pdf
Скачиваний:
11
Добавлен:
19.10.2023
Размер:
3.4 Mб
Скачать

 

 

 

-4 0

-

 

 

В принципе вычисление lie X X

X )

и h s(£ s, £ч)

мокко

выполнять

и другими

способами -

методом Ньютона (обратного якобиана)

и простой

итерации.

Последняя в

сочетании с методом определяющих вели­

чин может выполняться по следующей

простой

схеме:

 

из первых двух уравнений

на третьего и четвертого

Т - (К)

_ Гл - )

іJ/-с - — f-Jzcг ( к )

) х

1 ZS

г

j(trH)

Л2с(*>

1 to

(К) X

т(ы)

ІА2 S

J-IS

В случае расходимости итерационного процесса придѳтоя перейти к уравнениям узловых напряжений (при некотором виде характеристик НЭ на­ до начинать о этого метода):

 

 

ff Uic + bU3s~ CJ) Е/с

6/E)s

ffz Егс

EzEzs

О)

 

 

 

где

 

b &3g~ff

Е/с ~~ff/Е/s+ бг Ezc ~ffz Ezs ~0,

 

 

 

 

ff^fff+ffSffi

,

6=ö,+Bz+ö3 .

 

 

 

Г(ЛГ)

 

ГТ (к)

,

Задавшись

К

-ми приближениями

определяющих величин

U3c

 

U3S

вычислим

К

-тые приближения значений

всех

параметров

согласно

схеме:

 

Ѵ зс -и з

 

ff3

Elc

U3c Ufc

U,

 

g ,

I

 

E2c U3c= D~lc '-Ui '-ffг ,

X

 

X

 

 

X X

 

 

 

 

 

X X

 

U3s—^ff3~*83

, E ,S ~U3S-U/s-^U,~^6,

;

 

EZS-U 3S=U3S—•-ffi —

öi .

 

 

 

 

 

 

Uffo

-тых

приближений

 

Шс

.

ѵ £

 

Подставив числовые значения£/*

 

 

 

 

 

 

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

(

К

+

и

£ІК>

.

Далее с

помощью одного

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

 

1)

приближения ведающих величин.

В случае простой итерации

вовможна иная

схема

 

- без вычисления

невязок:

подставить числовые значения числения невязок: подставить числовые зна­ чения К -г о приближения параметров и решить систему относительно за ­ дающих величин, что ң даст их ( К + 1) приближения.

-4 1 -

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

или

//

(£,, 6г ,

■ ,

Сі , ..

 

С ,) ,

i=-f

Cf =

min,

иное

( £і,

Ег.

Сі,

■■

 

г

 

min

 

 

 

 

 

 

 

.

 

7

1£ r l=

 

или

 

 

Сг, . .

 

 

 

 

у 1

 

J J

(Сі,

 

Ci,

. -

 

‘ - >

±

Сі -

min

Здесь

&і , . . . .

€ „)

-

называют функцией минимизации (ми­

нимизирующей функцией). Наилучшую сходимость процесса минимизации (не­ избежную для линейных систем) обеспечивает функция первого вида, а наихудшув - третьего вида. Несмотря на ѳ то , вторую и третью также применя­ ют в отдельных случаях, когда их простота оказывается выгодной.

Дадим представление о процедуре минимизации на примере решение, уравнений узловых напряжений:

д Ufr + 6 U3S - д , Ею ~ 6 іЕ/3 ~Qz Ezc ~3Z E2s - О ,

öU ic-gU is* 6,Е,с - д! Efs +ЗгЕ2сQzEzs = <?.

Последовательность

 

операций:

IrJ f

 

,

вычисляем соответствующие

1 .

Задаемся значениями

 

 

,

 

 

 

им значения

параметров

 

g f',

6,'°’, д'с),

 

6/°’, д

 

öjc>,

все подставляем

 

 

 

 

£гм

 

 

 

 

 

'°’ ,

 

 

в уравнения иттвычисляем«Ю) 7тСот

 

 

,

 

и (£/';]

+

 

/с>)г

 

,с>

 

2 . Оставляя неизменным значение

U E

,

берем

последовательно зн а-

 

77 » 0

 

Use

 

uir

I

 

 

 

 

л

 

 

 

 

 

 

 

 

 

 

чѳния

Use

,

,

и т . д . ,

вычисляем

 

каждый раз

значения параметров

IJ,

бІ1<ат)ъ

значение

 

и !'т°,

пока не

будет

получено

ее

минимальное

зна­

 

 

 

 

 

 

 

 

 

 

)

min

.

 

Значение

U

 

при котором

зтот

чение (первый частный минимум

 

 

 

ä c ,

минимум достигнут,

считаем

первым приближением

Djf’

 

 

3 .

Оставляя неизменным значение

 

 

 

,

берем

последовательность

 

значений

 

 

 

Uj°z)

,

b

. ,

Uf°m>

, . . .

 

,

вычисляем

каждый раз значе­

 

 

gli,ow),

 

 

 

ния параметров

 

 

,

 

; " om),

и значение ^ у / ^ ^ п о я а

оно не будет

ми­

нимизировано

(второй частный

 

 

)jz

 

=

min

.

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

минимум

 

 

 

 

 

 

- 42 -

значение U}s считаем нервны приближением

Uls

.

Процедура повторяется

до достижения мннжыума-миниыорума,

о чем

будет

свидетельствовать

невоз­

можность дальнейшего уменьшения

j j

=

£f+

 

ни

ивменѳниѳм

Uic

, ни

изменением

U3s .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Конечно,выполнение описанной процедуры минимизации, как и ранее описанных процедур итерации, лучше возложить на ѲЦШ . Для не самых простых вадач это неизбежно.

10. О РАСЧЁТЕ БНЭЦ, НЕПРИВОДИМОЙ К ЭЛЕМЕНТАРНОЙ

задач;

Эта задача трудней, и подчас значительно, чем аналогичная задача расчета ИНЭЦ. В зависимости от условий задачи (степени ее нелинейности, числа НЗ) и от требуемой точности может оказаться целесообразным при­ менение методов различной сложности и громоздкости. Приведем некоторые рекомендации:

а) Метод' ѳквивалѳнтных синусоид целесообразно применять, когда нелинейность слабая , а высшіе гармоники не имеют принципиального зна­ чения для действия рассчитываемого устройства. Такой подход применяет­ с я , в частности , в цепям, содержащим слабо насыщенные стальные сердеч­ ники (дроссели , трансформаторы, иногда влѳвтрическиѳ машины). Пример применения приведен в п . 11 .

б) Методы малого параметра применимы к цепям со слабой нелиней­ ностью. Для периодических процессов ив втой группы применяется метод вовиущѳний [ 3 , 4 , 1 0 ].

в) Метод гармонического анализа более точен и предложен для цепей с умеренной нѳдинѳйноотью, т .ѳ . для цепей не содержащих НЗ с резкими изломами характеристики: ключевых схем (диодов, триггеров, р ел е ), маг­ нитных влѳмѳнтов (фѳрросердечннков) с прямоугольной характеристикой, сначала выполняется расчет по первой гармонике, что повводяѳт затем провести расчет по высшим вармоняхам. Для выполнения расчетов необхо­

димо

иметь

г а р м о н и ч е с к и е

х а р а к т е

р и с т

и к и

(ГХ)

всех

НВ цепи. Представление о получении ГХ дано

в п .1 4 .

Распола­

гая гармоническими "характеристиками всех НЭ для первых гармоник, можно провести расчет по первой гармонике полностью аналогично расчету ИНЭЦ, изложенному в п . 8 ,9 . Второй этап заключается в расчете цепи по высшим гармоникам, выполняется в соответствии с [ 7 ] и изложен в п . 14. В от­

 

 

 

- 43 -

 

 

 

 

 

по первой

дельных случаях ограничиваются лишь первый этапом - расчетом

гармонике, если достаточен грубый результат.

 

 

 

того же

г)

Методы дискретных значений функций

применимый к цепям

класса,

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

спо­

собов вычислительной математики и для ввода в ВЦВМ. Иввестны попытки

применения к цепям, содержащим диоды, давшие сравнительно

грубые

резуль­

таты. Эффективен один из методов ѳтой группы, называемый

т

о

ч

е ч ­

н ы м

[1 7 , 53• Основные сведения об одной

из форм точечного

метода

приведены в п . 15.

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

д) Метод

припасовывания

аппроксимации.

Пригоден для

цепей также

я- о

сильной нелинейностью. Эф­

 

 

фективен для расчетов переходных процессов, но для периодических про­ цессов применение часто встречает затруднения. Вопросы применения раз­ работаны в [12З . Элементарные примеры приведены в п . 13.

ѳ) Графо-аналитические методы иногда позволяют существенно упрос­ тить исследование и расчет НЭЦ. Среди них следует назвать фундаменталь­ ный метод фазовой плоскости [1 5 , Ю З и частный метод Волынкина [ 3 , 4 3 .

ж) Методы частного значения применяются каждый для какого-нибудь увкого класса задач. Примером такого метода является изложенный в п .12 "Метод условной линеаризации".

Время, отводимое для изучения данного раздела, позволяет изложить вдѳсь лишь некоторые ив названных методов, отобрав наиболее употреби­ тельные и доступные для учащихся.

11. РАСЧЕТ НЭЦ, СОДЕРШ Щ ШБОНАСЩЕНШЕ РЕАШИВНЫЕ КАТУШКИ И ТРАМОВІАТОРЫ

Поставленная задача наиболее просто (хотя довольно грубо) решает­ ся с помощью метода эквивалентных синусоид, упомянутого в пункте а) предыдущего параграфа. Дадим предварительную характеристику втому ме­ тоду.

Расчет НЭЦ практически всегда выполняется с погрешностями, обуслов­ ленными принятыми допущениями, облегчающими расчет. Это так называе­ мые погрешности метода, кроме которых всегда возникают погрешности собственно вычислительных операций. Последние принципиально могут быть сделаны как угодно малыми, но методические - могут быть весьма значи­ тельными, о чем иногда приходится мириться, ибо каждое последующее

-44 -

1точнеиие достигается иного более дорогой ценой, чем предыдущее. Поэтолироко црмнѳняются в грубые расчеты, кажущиеся мало обоснованными.

‘днни яз таких

примеров

явдяѳтоя м е т о д

э к в и в а л е н т н ы х

с и н у с о и д .

Этот метод допустимо применять лишь в тех случаях,

«>гпа высшие гармоники

не определяют принцип действия анализируемого

устройства. Он основывается на расчетной замене несинусоидального пе­

риодического процесса зквивалѳнтным, no действующим значениям

и мощнос­

т и ,

синусоидальным процѳосом согласно следующим формулам:

cosy*

 

 

 

 

f o s ^ = У 1, К /■■І

3

 

 

 

 

иэ

 

ей

 

 

Хза)и,~иСл

 

где

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

 

 

и

Сл

- линейное приблихеннѳ нелинейных параметров.

 

 

 

 

 

Практически вто означает, что генерирование высших гармоник нелинейными элементами учитывается косвенно и весьма приближенно, а формально рас­ чет ведется теми же методами, которые применяются для расчета ИНЭЦ.

Вподце целесообразно применение метода эквивалентных синусоид к БНЭЦ со слабой нелинейностью, в частности к НЭЦ, содержащим слабо насыщенные реактивные катушки и трансформаторы*! Чем слабей насыщение, тем более точные результаты дает метод эквивалентных синусоид.

В цепь, содержащую реактивную катушку или трансформатор, нелиней­ ность вносит фѳрросѳрцѳчник. Поэтому основная цель заключается в пост­ роении характеристик эквивалентной нелинейной электрической схемы (сх е ­ мы замещения), приближенно описывающей энергетический процесс в фѳрросѳрдечникѳ. Располагая такой характеристикой для действующих значений,

можно далее

вести расчет

ИЗЦ т а к , как показано

 

в п .

8

и 9 .

 

Электрическая схема

замещения имеет вид,

представленный ьа

 

рис. 24, а , где нелинейная активная проводимость

дс

(

) описывает

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

и вихревые

токи,

т .ѳ . условия

преобразования в теплоту, а нелинейная

индуктивная -

 

проводимость

6а { Ua) -

описывает условия

создания

в

сердечнике

маг­

нитного потока. Задача заключается в том,

чтобы построить характѳристн-

Именно такой подход повводил построить классические теории трансформаторов и ѳлѳктричѳских машин переменного тока.

- 45 -

ки lg (b~o),lt (Ub а ватем I 0Wo) , гд a I0 = j Ід ± і\ . Для ѳѳ решения необходимо располагать данными об энергетических прѳобрааованиях в ма­

териале фѳрросердѳчника, которые находятся

иа

эксперимента и приводят­

ся

в справочниках. Это

зависимости

Р'(Вт)

 

и

О'(бт) ,

 

заданные

таб­

 

 

 

Р'(бш )

 

лично или графически,

как на рис. 2 4 ,6 .

Черев

 

 

и

Q'Cß/n)

обо­

значены активная и реактивная мощности, приходящиеся на 1 кГ веса

У

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

 

д а н н о й

 

ч а с т о ­

т е

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

Вт

.

 

Учтем известные свяви:

*■

ü0ЧКф$- wSßm= üßm'

где У - вес сердечника рассчитываемой катушки (трансформатора);

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

(X

- постоянная, значение которой определяется коэффициентом формы

 

кривой

Рср=

1,11 для синусоиды,

частотой

j

, числом витков

 

первичной

обмотки

иг

Ри поперечным оечениѳы магнитопровода ^ .

 

 

Очевидно, что зависимости

'(Вт )

и

Q'(Bm)^

можно рассматри­

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

P '(U .)

 

и

O '(U0)

і

поскольку

В т - Ш

.

 

 

 

 

 

 

 

Поэтому

- 46 -

lg = P'([/„)■ - Ц - = lg (Uo) ;

Iâ-Q'(U,>%a- I e(V.)

могут быть построены по точкам, для чего надо аадаться последователь­

ностью

U0

и вычислить соответствующие вначѳния

Ід

м Tg ,

Соответствую­

щие графики покованы на р и с. 25, а .

Для последующего расчета цепи,

со ­

держащей феррооѳрдечник, удобней подьвоваться характеристиками

І 0

( К )

и

(-f0 ( U0),

которые элементарно строятся по уже известным:

 

 

 

 

 

 

 

 

 

Ь (Uo)=yll(Uo)+ll(U0) ;

% (0о)= artig

 

 

 

 

и показаны на рис. 2 5 ,6 .

 

 

 

 

 

 

Ыа рис. 26 покавана цепь, состоящая ив эквивалентной схемы транс­ форматора с нагрувкой. Нетрудно видеть, что эта цепь удовлетворяет сформулированным в п . 7 условиям и может быть рассчитана по методике, изложенной в п . 8 . В случае НЦ большей сложности, т .ѳ . содержащей бо­

лее одного эквивалентного источника или более одной определяющей вели­ чины ( Э>1,5>1 ) , расчет может быть выполнен ооглаоно п . 9 .

-47 -

12 . МЕТОД УСЛОВНОЙ ЛИНЕАРИЗАЦИИ

Применение этого метода покажем ва характерном примере цепи, представленной на рис. 27 и описываемой нелинейным дифференциальным уравнением: ^

 

 

 

 

 

 

 

+ г! = U(t),

где и (і) = ~]/zUSin cot.

 

 

Зависимость ЦЧі) нелинейна, вадана

 

 

 

 

 

 

таблицей-графиком или аппроксимирующим

 

 

 

 

 

 

выражением,

например как

 

 

 

 

 

 

 

 

 

 

 

 

L (</>)=а,<р + а}у 5+абц>5.

 

 

 

 

 

 

 

 

 

Если в уравнении один член в неко­

 

 

 

 

 

 

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

 

 

 

 

 

 

І™ 1

макс

^

\ju\ макс^ ,

то

проще

 

 

 

 

 

 

 

 

 

 

приближенное

 

 

 

 

 

 

всего

решать

уравнение

 

 

 

 

 

 

IHL± Li(t)

»

</’ ä

u(t)

и затем

найти

 

 

 

 

 

 

^

и //)

 

 

 

 

 

 

 

 

 

 

 

 

по известной аависимости.

В неко­

такой

прием допустим,

но

торых

случаях

при

очень сильном неравенстве

значительно

целесообразней

не

пренебрегать членом

гс

, а учесть его

приближенно

как линейный,

т .е . л и н е а р и з о в а т ь

 

м а л ы й

член.

Для этого при переходе

в уравнении к одной

переменной

величине

надо ограничиться линейным членом этой связи . Выражения

і

черев

(р,

 

 

Jd'

 

^ Очевидно,

что неравенство

7t '-7 ~[7?

не может выполняться,

т .к .

■ jj:

 

 

 

 

также

может проходить через

нуль.

 

 

 

 

 

 

- 48 -

ограничимся вависимостью t —Ot(p=--j— ф , т .к . коэффициент О/ имеет смысл и размерность, обратные статической индуктивности. Получаем урав­

нение

л и н е й н о е

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

п о т о к о с ц е п -

л ѳ н и я

 

 

 

 

 

Ult',

решением для тока в

решение

которого

 

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

линейной цепи

г

,

L

:

 

 

 

 

 

і ли

y>L= a fctg

V

 

 

 

 

 

1/rz+{uLtf

 

 

 

 

 

 

 

где у* - угол сдвига фею мѳаду W и U . Потокосцѳплениѳ будет ивнѳнятьоя во времени синусоидально:

5P ( t ) = V T Jfr-s tn (o jt-y L) i

а ток, найденный по зависимости

i(t)= O t< p (t) + Q3 (p3(i)+ Os (ps( t ) ,

будет содержать первую, третью и пятую гармоники, в чем нетрудно убе­ диться, выполнив подстановку и преобразования подобные тем, что пока­

заны в п .

6, в

пункте

б ) .

 

 

 

 

 

 

 

 

 

 

Если линеаризуемый член достаточно мал, то нелинейные соотношения

учитываются при решении достаточно полно.

Погрешность

расчета

возрас­

тает при неудачном выборе коэффициента

0/=

^— ,

что

зависит

от выбо­

ра предполагаемого рабочего участка, т .е .

участка

применения

аппрокси­

мации. В том случае ,

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

(далеким

от режима,полученного

в результате

р асчета), Оможет

оказаться

необходимым повторить

р асчет,

приняв лучшее

значение

/ .

 

 

когда

выполняет­

Рассмотренный метод может быть применен и тогд а,

ся неравенство

1^\иокс^ гі

ткс

. В этом

случае для малого члена лине­

 

|(і) |

,

 

аризуется

зависимость

 

например

 

<pti)= ö,L = LfijL.

 

 

 

Ф(ё)= 6,і -

,

 

как

 

 

 

 

 

 

 

 

 

 

 

 

 

 

так как коэффициент 6Г= имеет смысл и размерность динамической индуктивности. При зтом уравнение цепи, линейно относительно тока

U sin ult,

который легко определяется как

- 49 -

 

где

 

i(t )= V z

I

Sin (U)t~<p),

1 = ,

^

-- ,

■-

<J>= arctg СОил

 

Vr*+(uLt1f

 

У

Затем по линейной

еависиыоста(

можно найти

 

 

 

р (і) =

0 ,і - 5 з ь 3.

 

Заметим, что метод условной линѳаривации можно рассматривать как

линейное приближение метода вовмущѳний (например, [ lÖ ] } .

не

Если ни один ив членов нелинейного дифференциального уравнения

является малым,

это

метод

условной линѳаривации непригоден. При

втом можно воспользоваться методом возмущений, в котором расчет ведет­ ся с учетом не только линейного члена, но и членов более высоких сте ■ пеней.

13. О ПРИМЕНЕНИИ КУСОЧНО-ЛИНЕЙНОЙ АППРОКСИМАЦИИ (МЕТОД ПРИПАСОШВАНИЯ)

Использование кусочно-линейной аппроксимации позволяет заменить нелинейное уравнение цепи (а в случав системы нелинейных уравнений - каждое ив уравнений) сиотѳмой линейных уравнений, каждое из которых описывает процесс лишь на одном линейном отрезке (куске) характеристи­ ки. Написать эти уравнения и вид их решения легко , но обычно оывает трудно выполнить сопряжение этих решений на краях соседних участков, удовлетворяющие условиям задачи, в частности законам коммутации. Эта

часть решения, как и метод

в целом, называется п р и п а е о в ы в а -

н и ѳ м (подгонкой одного

к другому).

Несмотря на трудности

применения, этот метод имеет два достоинст­

ва: дает решение в аналитической форме и , главное, позволяет исследо­

вать влияние самых сильных нелинейностей,

которые и выражаются излома­

ми характеристики на большие углы.

 

 

 

 

 

Рассмотрим простые примеры применения. Легче всего зтот метод

применить в случае,

когда участки

характеристики НЭ

совпадают с осями

координат или параллельны им.

,

L(i)

по

рис.

28, а аппроксимиро­

Пусть индуктивность в цепи 7’

 

вана характеристикой

Ц](і)

согласно рис.

2 8 ,6 .

В цепи действует сину-

 

Соседние файлы в папке книги из ГПНТБ