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

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

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

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

Добавлен: 16.10.2024

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

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

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

Если сечение намагничивающей катушки прямоуголь­ ное с размерами Ак— па и Нк=та=$па, то при принятой кратности коэффициентов теплоотдачи ут— спра­ ведливо

=

2рV (1 +

J - ) ; *°к = р/г2.

(2-139)

Приложение

полученных

результатов

(по модели

«0О») к проектированию реальных исполнений электро­ магнитных механизмов с втяжным и внешним прямохо­ довым и поворотным якорем подробно изложено в [Л. 51].

М о д и ф и к а ц и я м о д е л и «0П» (рис. 2-15,г). Теп­ ловая модель «0П» в отличие от модели «0о» учитывает неравномерное распределение превышения температуры

&*(z,y) по сечению окна S0K обмотки, а следовательно,

дает возможность определить наиболее (максимально)

нагретые участки (точки с координатами уи, zM) и пре­

вышение температуры

в ней 0*M(zM, ум),

что существен­

но при оценке степени

нагрева катушки,

а также дает

возможность

найти

среднее

превышение

температуры

 

е*СР

 

[J 9*

(У, z)dy dz,

(2-140)

 

 

 

\s)

 

 

 

которое определяет и. с. катушки.

для

модели «0П»

Согласно

принятым

допущениям

(рис. 2-16,а)

нагрев

ее определяется

уравнением тепло­

проводности

(9=8* ,

(9=0*

qcP

 

 

 

 

(2-141)

 

 

 

 

 

 

инулевыми граничными условиями первого рода. Практический интерес представляют также случаи, учи­

тывающие отсутствие теплоотдачи с одной (рис. 2-16,6) или двух (рис. 2-16,е) поверхностей. Как будет показано ниже, анализ указанных дополнительных вариантов (мо­ дели «0п-^‘/2» и «0n-HV4») дает возможность определе­ ния максимального и среднего превышения температуры практически при всех реальных несимметричных случаях нагрева и охлаждения катушек.

Решение поставленной задачи в части определения интеграла (2-141) будем искать в виде двойного беско­ нечного ряда

в* (У, z) = 0 — 0П= £

£ Bij sin (Pi у) cos ( S j z), (2-142)

f = i

/ = i

2 0 0


где Bij, Pi Sj — подлежащие определению постоянные, значения которых находятся с учетом граничных усло­

вий рассматриваемых моделей.

как показано на

Расположив систему

координат,

рис. 2-16,г,

 

 

гдей = -^ -# к и с = - ^ - А к,

укажем на

очевидный факт

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

ной четверти окна намотки (например, S0h=<Si234 = c6), определяющийся симметрией модели «0П». В этом слу­

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

grad Q*y |у=(;= 0 и grad 0*z|2=o= O равносильны условию отсутствия теплоотдачи на этих поверхностях и сохра­ нению принятой картины изотерм в сечении. Как видно, в этом случае максимальное превышение температуры 0*м соответствует значению координат уы — с и zM= 0.

Среднее превышение температуры 0*ср по окну на­ мотки при заданной постоянной производительности внутренних источников нагрева qCJ>, учитывая симмет­ рию, сохраняет свое значение для случая рассмотрения температурного поля в половине окна (SOK= Si265i= 2£c), а также всего (полного) окна (S0K==*Si578i=4£'c) обмотки катушки. При этом исходя из принятых условий посто­ янства и равенства температуры 0 П на наружных по-

201


 

верхпостях,

т. е. при 0 < //< с ,

z = + b

и у = О, O ^ 'z ^ b ,

 

следует 0==0П, и, если учитывать условия равенства ну­

 

лю

градиента температуры

на

внутренней

поверхности,

 

т. е. при у= + с,

O ^ z ^ + b

и Os^i/sCc,

z=0

следует,

 

например,

равенство

grad

©*у 11/=с = 0,

можно,

как это

 

было показано в '[Л. 50, 51],

найти решение задачи в об­

 

щем виде:

 

 

 

 

 

Щ!

 

 

 

 

 

 

l6gcPf2

sin (2i — 1)

c h (2 t - I ) - g -

 

0*(//, z) =

2с

\

П3\*э

Б

(2t — 1)3

 

c h (2 /-l)^

 

 

 

 

i=i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2-143)

 

В этом случае максимальное значение превышения

 

температуры определяется (рис. 2-16,г)

при уы = с и zM=

 

= 0 в следующем виде:

 

 

 

 

 

 

0*.

Q с р с ~

1,03

£

 

 

Tzb

 

(2-144)

 

2Х*П

 

 

 

 

 

 

 

 

 

1=1 (2i-

\ y c h ( 2 i - \ ) w

 

 

 

Аналогично среднее значение превышения температу­

 

ры по (2-140) для рассматриваемого случая определяет­

 

ся в виде

 

 

 

 

 

 

тzb '

 

 

 

 

 

 

 

 

 

th (2/

 

 

 

0*ср ■

<?срГа

 

1 —

0,625

Е

*) 2Г

(2-145)

 

 

 

 

 

 

зк*

 

Ь/с

(2f — !)■

 

 

 

 

 

 

 

 

/=|

 

 

 

 

 

При этом для приведенных моделей температурного

 

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

 

 

 

 

 

 

Модель «0пч-74» (рис. 2-16,а):

 

 

 

 

 

 

с= Ак\

b= HK\ Ь/с — Нк/А1{ = $.

 

(2-146а)

 

Модель «0П“i- V2» (рис. 2-16,6):

 

 

 

 

 

с= Ак; 6= Як/2; Ь/с = Нк/2Ак=$/2.

 

(2-1466)

 

В общем случае для модели «0П» (рис. 2-16,а)

 

 

с— Ак12; Ь = НК/2; 6/с = 2Як/2А1(= [3.

 

(2-146в)

 

Используя затем принятый принцип базовой и кор­

 

ректирующей функций, выражения (2-144) и

(2-145)

 

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

 

 

 

 

 

0*.

4сИк

 

_ 4с[Ик

 

 

(2-147)

 

 

8\ *

M l

ср ■ I2X*. Н'ср»

 

202


где первые приближения (базовые функции) состояния приняты равными

<?сИк .

Z q. ср — (9*ср)б

9сР^К

(2-148)

Zб.м = ( 0 * м ) б

Т2Х*7’

а корректирующие функции, как следует из (2-148), (2-147) и (2-145), (2-144), определяются с учетом (2-146)

следующим образом: для модели «0ц-ИД»

для модели «0ц -е- 1/2»

(2-149)

 

^ • = 4р - - т 1 ‘|1т Р ) ;

<2-150>

для симметричной модели «0П»

 

= 1 — 1,03 sch (3;

(2-151)

^ o p = l - ^ - t h - f p .

(2-152)

В связи с быстрой сходимостью ряда в

(2-144),

(2-145) при определении корректирующих функций при­ нято i = l.

Укажем дополнительно, что аналогичная задача, од­ нако в более общем виде и в том числе при производи­ тельности источников, зависящей от температуры для различных моделей, рассмотрена нами в (Л. 51—53], там же приведена методика исследования и эксперименталь­ ные данные по температурному полю различных испол­ нений катушек, подтверждающие полученные расчетные результаты. Рассмотренная задача может быть решена и другими методами (Л. 1], в том числе и приближенны­ ми, например, так, как это сделано в [Л. 32].

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

203

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

В связи с этим рассмотрим постановку вариационной задачи несколько подробней уже в данном разделе. Из­ вестно, что если функция 0* (у, z) доставляет экстремум функционалу типа (§ 2-3)

I [0* (у, z)\ — J f R (у, z, 0*, р, s) dy, dz,

(2-153)

Ф)

 

где / — функция от функции @*(y,z) \ D — область зна­ чений функции 0*, на границе которой указанная функ­ ция постоянна или равна нулю; р и s — принятое обо­ значение частных производных:

дв*

дв*

(2-154)

р

фГ;

s ~~~dT’

 

то отыскание указанной

функции 0*(г/, z) — экстремали

равносильно нахождению решения 0*

уравнения Эйле­

р а — Остроградского [Л. 97]:

 

 

dR

 

 

 

дв*

 

 

<2' 155>

где

 

 

 

 

 

£{£}

( 2 ' 1 5 6 )

являются полными частными производными функции R соответственно по координатам у и z.

Возможна и обратная постановка задачи, т. е. отыс­ кание решения Q*(y,z) уравнения (2-155), что равно­ сильно нахождению экстремали 0*, доставляющей функ­ ционалу (2-153) оптимальное значение. В рассматривае­ мом случае уравнение Пуассона вида

д*в*

, а*8*

,

?сР _ п

(2-157)

ду2

> ' дг*

\*э

 

при граничных условиях первого рода является уравне­ нием Эйлера — Остроградского (2-155) и противопостав­ ляется функционалу (2-153), который имеет вид:

' =Ш (тг)'+(^)’ - 20* г]

<2-158>

204


где 5 — область интегрирования, в нашем случае — сече­ ние окна намотки.

При необходимости учета анизотропии по теплопро­ водности вдоль координатных осей {'К гф 'ку) и возмож­ ной линейной зависимости удельных источников нагрева от температуры вида <7= £7о.с (1 + ат0*) задача решается при определении экстремали функционала вида

/ = j f [ яг ( ^ - ) 2+ Яу ( ^ " ) 2— 2<7о.с0*—2aTq0C(Q*)2^dy dz.

(S)

(2-159а)

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

+4^ ) ’ -2<7о,с0 * — 2<7о.са т (9 * )2 г dr dz.

(5)

(2-1596)

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

©*{у, г ) .

Из прямых методов наиболее широко используются метод Ритца и метод Канторовича {Л. 34]. При примене­ нии первого метода к функционалам /[0*(xi, х2, ..., хп)], зависящим от функции нескольких переменных, подби­ рается система функций

w i ( x u х2,... , х п), W 2 (X i, Х2,...,Х п ),...

. ...

W i {хи Х 2, . . . , х п)

 

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

 

 

i

 

0*г = 2

ак®к {хи х 2......х п),

(2-160а)

&=1

 

где коэффициенты ак — постоянны.

При использовании второго метода коэффициенты aK(*i) не постоянны и являются неизвестными функци-

205