Файл: Любчик, М. А. Оптимальное проектирование силовых электромагнитных механизмов.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