Файл: Любчик, М. А. Оптимальное проектирование силовых электромагнитных механизмов.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 16.10.2024

Просмотров: 141

Скачиваний: 7

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

функции, входящие в комплекс Кв(£, х) без учета стали. Полученные значения K jbi(-yi) и принятое x = Xi дают возможность определить индукцию В\, что в свою оче­ редь дает возможность скорректировать совокупность функций

k (x u В ,)= о (-*,, В,). В■(хи В,), ср0(л-,, В,), 30(х., В ,)...}.

входящих в комплекс Ki(.v'i, Вi).

/

Р и с . 3 - 1 .

Таким образом, можно получить более точную зави­ симость Dx— fD(x, В) н тем уточнить графоаналитиче­ ский расчет, который применительно к характерным си­ стемам постоянного и переменного тока подробно изло­ жен в [Л. 50, 51].

Итерационный метод расчета определяющего размера

чЗначительные преимущества с точки зрения точности

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

Уравнения (3-45) или (3-53), определяющие значе­ ния х при принятых постоянных н заданных кратностях геометрических размеров, по сути являются функцио­

нальной зависимостью вида

Х= [а(х).

(3-56)

258

Начальное (нулевое приближение) значения корня этого уравнения можно получить, приняв К(£, х )= 1 , в этом случае

х«) = { / С 0ГОЛ.

(3-57)

Более точный результат (первое приближение) полу­ чим, подставив в (3-56) значение х<°\ т. е.

X(,>= fa [-V'(0>J-

(3-58

Повторяя этот процесс, т. е. итерируя несколько раз, можно получить (г+|1)-е значение корня

*0=*<ж )=/а[*<•■)],

(3-59)

которое будет отличаться от предыдущего на величину, не .превышающую заданную точность отклонения:

и. (3-60)

При этом процесс будет сходиться, если выполняется достаточное условие [Л. 28]

dfa (

х )

< 1

(3-61)

dx

 

на отрезке изменения х в границах xi<x<xz.

Если условие (3-61) не выполняется, то (3-56) всегда

можно представить в виде

 

х = х —су(х),

(3-62)

где у(х) = f D(x)—Dc, а константа с подобрана так, что­ бы для функции

fa (x)= x —cy(x)

(3-63)

условие (3-61) имело место. Тогда согласно (3-59) полу­ чим:

,г0=д:(,'+1).=л:^—су [х<!’1].

(3-64)

Если при расчете определяющего размера необходи­ мо учитывать насыщение участков стали, например при малых критических зазорах или ферромагнитных шун­ тах в рабочих зазорах, то процесс итерации включает расчет индукции. В этом случае исходные зависимости х и В по (3-53) при принятых постоянных и заданных кратностях | включают в себя корректирующие функции

17*

259


k = k { x, В ) и, следовательно, являются зависимостями вида

 

х =

С ° K [ k ( x ! В)]

=

В

(3-65)

 

В==

CDKB [kjx.

В)]

f,(AT.

5).

(3-66)

 

х‘Гв (1)

=

 

 

,

 

 

 

 

При

этом начальные значения

(нулевые

'приближе­

ния)

 

 

 

 

 

 

 

 

 

•»"" =

/ с т а ;

в м =

г,(бн »..ц .-

<3-67>

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

приближе­

ние

Х(')= /Й:[#),

В(О];

5(i)= fB[x(‘),

В(0)],

(3-68)

 

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

 

 

 

 

 

 

 

*<*)= /в[х(1>,

fid)]; B(2)=fB;[х(2), B(D]

 

и т. д. до (г'+1)-го приближения

 

 

 

 

 

 

 

 

 

В<*>],

 

(3-69)

когдадостигается

 

заданная

 

точность

расчета

|^i+i)—

и соблюдаются ограничения типа

 

В с т =

Я (м - 1 )=

/ в [*<<•+1),

 

 

ас.

( 3- 70)

Указанный метод итераций может быть распростра­ нен при синтезе СЭММ по (3-12) или (3-14) при после­ довательном определении и других параметров, напри­ мер обмоточных данных по функциям fs и fw. Возможен контроль сходимости итераций по комплексному крите­ рию

{С, (*(*•> _.*и-ы>]* + г;2 [s^~ - Si£+I)]2 +

+

lk!

~

'’I2+ ... + Cn

~ В « +"\}'12< и ,

(3-71)

где

£i,

tz,

£з, . . U — соответствующие коэффициенты

веса.

Если процесс итерации расходится или затруднитель­ на проверка по (3-61) в общем виде, полезным оказы­ вается метод деления отрезка пополам (метод дихото­ мии) .

260


Метод дихотомии при расчете определяющего размера

Преобразуем (3-53) или (3-54) к виду

 

' У = У ( х ) Ч в ( х ) —й с,

(3-72)

где

Ы * )= * 'В Д , х)\ Я с = СоГ(£).

Затем, выбрав реальный диапазон изменения опреде­ ляющего размера х Xi<:xsg;x2, определим значение функ­

ции у в граничных

точках — концах отрезка [xi, x j

(рис. 3-1,6) lJi = y{Xi)\

у2 = у(х2).

Если при этом одно из полученных значений по абсо­ лютной величине меньше или равно заранее заданной

степени точности поиска корня

 

 

Ы * ) |< « ,

(3-73)

то корень определен случайно.

осуществлено, проверя­

Если соотношение (3-73) не

ются знаки функций yi — y(xi)

и уг=у(хг)

на концах

отрезка [яу, хг]. Если произведение функции больше ну­

ля: (t/i,

yz)>0, то знаки

г/j и у% одинаковые (рис.

3-1,6)

и, следовательно, корня

нет пли он

не

единственный.

В этом

случае нужно изменить границы отрезка.

то ко­

Если

произведение меньше нуля:

(уи

уг) <0,

рень находится в принятых границах. В этом случае производят деление отрезка пополам и определяют:

*з= (a'i + .v2)/2.

(3-74)

Затем проверяется

абсолютное

значение функции

при х3

\y{xz) |< « .

 

У г=

 

Если заданная точность выполняется, то корень най­ ден: хо=х3, если нет, то сопоставляют попарно произ­ ведения функции на концах половинного отрезка

(У,У3) [ > ,

< ]0 ;

)

(3.75)

(У3У »){>,

< ]0 .

/

 

Если произведения указанных функций больше нуля, то корня в данной половине отрезка нет, так как функ­ ции у(х) одного знака. При значениях произведения меньше нуля (узУг) <0- корень находится в данном от­ резке, при этом вновь производят его деление пополам

х&=-(Хз+Хг)/2

(3-76)

261



и проверяют абсолютное значение функции по

| У(ха) | < ы

и т. д., пока методом деления пополам предыдущего от­ резка границы изменения определяющего размера не определяется корень хо, при котором достигается задан­ ная точность расчета

\y(xo)\ = \fD(xo)—D e \ ^ u .

(3-77)

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

Рис. 3-2.

262 -

указанных методов при синтезе отдельных групп СЭММ будут освещены ниже.

На рис. 3-2—3-4 приведены диаграммы алгоритмов, реализующие указанные методы расчета на ЭЦВМ. Как следует из анализа методов прямого синтеза, его реа­ лизация требует (по определению) предварительного за­ дания комплекса кратностей геометрических размеров ядра, следовательно, и задания совокупности основных

кратностей:

 

&={«, т'Ш е, v ...}.

(3-78)

Рис. 3-3,

263

Рис. 3-4