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

книги / Надежность и диагностика энергетических электромашин

..pdf
Скачиваний:
2
Добавлен:
12.11.2023
Размер:
8.33 Mб
Скачать

Рис. I . Области значений моментов Л/асин­ хронных двигателей при изменении диамет­ ра J) , допустимого перегрева IT и ин­ тенсивности охлаждения.

тов или

о твета,

что в заданном габарите нужный технический пока­

затель

получить

н ел ьзя .

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

Для определенного типа электрической машины в интересующем достаточно широком, диапазоне геометрических параметров, перегревов активных частей и интенсивностей охлаждения численным эксперимен­ том на ЭВМ по разработанному алгоритму расчета, описанному в ли­ тературе94, было получено поле реализуемых электрических машин. Анализ таких расчетов показал, что зависимость момента, развивае­ мого электрической машиной, от габаритных размеров при постоянной линейной скорости бочки ротора представляет собой в логарифмичес­ кой системе координат прямые линии для различных интенсивностей охлаждения и перегревов. На рис.1 представлена зависимость момен­ та М на валу асинхронного двигателя от его наружного диаметра J)

для шести значений интенсивности охлаждения, охватывающих’ следую-

хАнпилогов Н .Г ., Бандурин В ;В ., Вакуленко К .Н ., Вишникин А.И. Алгоритм проектирования частотно-управляемых короткозамкнутых асинхронных двигателей . - Техн.электродинамика,1981, А 4, с . 49-53.

171

щвй диапазон коэффициентов теплоотдачи на основной поверхности теп­

лосъема: 1 - 250 Вт/м2 °С; 2 -

500;

3

- 750; 2 -

1000;

5

-

1500;

6

- 2000 Вт/м2 °С. Такой диапазон

соответствует

всем

возможным зна­

чениям коэффициентов теплоотдачи в системах струйного масляного

охлаждения электрических машин. Зоны,

обозначенные на рис.1 номе­

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

равно-,

го

JT = 80°С, а

оверху - перегревом

200°С. Зоны I и 2-й

систем охлаждения

ограничены на

р и с .I

сплошными линиями, зоны 3 и

4 -

штриховыми линиями, зоны 5

и 6 -

штрихпунктирными линиями. Для

электромашин транспортного исполнения

принятый диапазон по

перегре­

ву охватывает все

классы применяемой изоляции. В исходных

данных

задается температура охладителя на

входе, следовательно,

по

вели­

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

го варианта величина перегрева.

 

 

 

 

 

 

Для каждой системы охлаждения расчеты выполнялись

для

соотно-

’ шения наружного диаметра электрической машины к

ее

длине

I

.рав­

ного

л/1- =

1 при следующих габаритах: J) - 0 ,3

м;

0 ,4

м;

0,5 3 м;

0,6

м; 0,7

м и при перегревах ^7 V 8 0 °C ; 120°С;

160°С;

200°С.

 

Выбранные диапазоны указанных параметров в

приводимом числен­

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

Если задана температура на входе и максимально допустимая температура изоляция, то выбор вариантов электрической машины це­ лесообразно проводить по графику, показанному на р и с .2, где для перегрева, равного лТ= 120°С, приведены зависимости момента от

наружного диаметра. Аналогичные зависимости могут быть построены для различных перегревов по данным р и с .1 .

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

172

Рио.2 . Зависимость момента, развиваемого асинхронным дви­ гателем, от наружного диа­ метра при допустимом перегре­ ве i T = 120°С для различных условий теплосъема.
0,5
Ц4 0,5 ' Q6 Дм

статочно эффективная система ох­

МЮ\Н-м

лаждения, или индукция насыщения,

 

стали, т . е . невозможно обеспечить

 

необходимые

размеры магнитопрово-

 

да. Очевидно, что проектирование

 

таких электромашин при габаритных

 

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

сводится к нахождению

 

оптимального

соотношения между

 

электромагнитными нагрузками я теп ­

 

ловыми проводимостями, приводящи­

 

ми к максимальному моменту на в а ­

 

лу. Предлагаемый ход проектирова­

 

ния с выбором по графику оптималь­

 

ного соотношения между габарита­

 

ми, энергетическими и тепловыми

 

параметрами

позволяет оператив­

 

но решать задачу выбора вариан­ та электрической машины для транс­ 0,2

портного средства, где компоно­

вочные решения могут изменяться

0J

в зависимости от ее выходных па­ раметров .

На основании полученных зави ­ симостей по задаваемым величине момента и габаритным размерам вы­ бираются параметры системы охлаж­

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

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

173

Ш: 621.3.012.6.001.24 : 681.3

Г .Л.Баранов

МОДЕЛИРОВАНИЕ ОБЪЕКТОВ С ПЕРЕДАТОЧНОЙ ФУНКЦИЕЙ ' ВТОРОГО ПОРЯДКА

Динамика объектов различной природы описывается дифференци­ альным уравнением или передаточной функцией. Известно, что переда­

точные функции высокого порядка представимы в воде произведения

передаточных функций первого

и второго

порядка,

моделирование

ко­

торых на ЭРМ выполняется с помощью рекуррентных соотношений.

 

 

При моделировании непрерывных объектов,

у

которых имеется

пе­

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

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

 

 

 

 

В первом случае корни характеристического

уравнения комплекс­

но-сопряженные

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

 

 

 

 

 

 

 

W,<P)

Ггр + ГР + (Г+б)г Г2)

 

 

 

 

 

(1)

 

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

тогда

модель

 

объекта ш еет

вид

 

 

 

к(Г,р+ f)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2)

 

 

 

 

 

 

 

Гр + I

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В третьем случае обычно корни разные действительные

числа

и

передаточная функция объекта эквивалентна передаточной функции

 

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

соединенных звеньев первого порядка

 

 

 

 

 

 

 

w ,

\

 

«

ггР + 0

 

 

'

 

 

 

 

О )

 

 

 

 

 

 

(Гр*1)(Тл р+1)

 

 

 

 

 

В качестве дискретного аналога непрерывной модели объекта в

этих случаях целесообразно

применять рекуррентную

формулу

 

 

где

А, В,

С, J> - постоянные

eYi-2i’CXi "-DXM ’

 

 

 

от

 

(4)

коэффициенты, зависящие

физичес­

ких параметров, передаточной функции объекта

и от

шага

А ,

опре­

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

(7 - /;

-

/ -

 

и

/ -м -

( / +

I)-м

состояниями

объекта

в

переходном

процессе;

 

/ / ,

-

входные

8начения

воздействия

 

соответственно

в

 

 

и ( / - 1)-й моменты времени;

Yh

 

Y;_,,

'

Y}_2

 

выходные

значения

реркции объекта на входное воздействие

соответственно

в

/',

U-D,

( i - 2)-й моменты времени.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Дискретная модель вида (4) полностью определяется численны­

ми значениями коэффициентов А,

3,

С, D ,

которые

зави сят

от

пере-

174

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

даточной функции объекта моделирования и от споооба описания вход­

ного

воздействия.

 

 

В зависимости от целей моделирования численные значения А .

S,

С,

V

меняются..

 

 

Определим численные значения постоянных коэффициентов А , S,

с

, J) ,

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

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

предварительного

интегрирования, примененному в работе / I /

для

объектов первого

порядка.

 

 

 

 

 

В качестве примера дадим аналитический вывод для объекта с

передаточной

функцией

 

 

 

 

 

 

 

. ___________ *

 

Y(P)

 

 

 

 

Шр)~ Тгр+2аГр +г

Х(р) '

(5)

где

р =

- оператор дифференцирования.

 

 

В этом случае дифференциальное уравнение переходного процес­

са

на интервале

/ < t* t.

имеет

вид

 

 

 

 

 

 

УСв

= к х С И ,

( 6) .

 

 

 

 

 

 

 

где x(f) - модель возмущающего воздействия на входе объекта, Y(t> -

реакция объекта на это воздействие.

Если сделать замену переменных г - 1 - t-j f • на расчетном интервале принять линейную модель входного воздействия

 

 

*(*)- *н +

 

г>

 

(7)

перейти

к изображениям по Лапласу /2 7

с оператором, то

получим

операторное уравнение

 

 

 

 

 

 

где х

Y а.

у \ =——— ■—

начальные условия или исходные

'■>

О-0)

/-V (+0)

п

 

 

 

 

 

значения соответствующих

переменных,

когда

г

стремится

к нулю

справа.

 

 

 

объекта y(S)

 

 

Искомое

изображение

реакции

из уравнения (8)

найдем

в ввде

у (3 )- ^

г*

 

,

+

 

 

+2afS+ ,

- Y(tg)

 

 

 

r*S+ 2 <xT

{ M S * * j+ H t/s -//4 s s) x ,lf

• .

+

T2S2 *2аУ5*/

r zs e -*■ PetГs +/

 

 

175

Используя соотношение

j

f ( s) +-*5

f ( t ) d r ,

 

осуществим

об­

ратное

преобразование Лапласа

/ 2 ,3 / , тогда

получим требуемые

ча­

стные

оригиналы

 

 

 

 

 

 

 

 

 

 

 

 

• S+ZaTS *1

 

 

 

 

 

 

 

 

 

 

 

r*s+ 2а Г

 

 

 

 

 

 

 

 

 

 

 

 

= в

Т

(cOS H it*

S in

< v t ) ,

 

 

 

 

 

 

r*s**taTS+1

 

 

 

 

 

 

 

 

 

 

 

s ( r * s 2 + M S + 1) = /- s

r

. (a> soir*~ r

s i n *

j r ) ,

 

 

 

$2(T*s*i-2ars*f) - = (t-2 a Г)* e r ^ZaTcostvr* (

-oiT2^

sin d/rj

 

Здесь

принято,

что &=

/- а 2 .

 

 

 

 

 

 

 

Из операторного уравнении (9) на основе полученных оригиналов

после соответствующих преобразований получим рекуррентное

соотно­

шение вида (4 ),

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

(5)

при

расчетном шаге,

равном

h «постоянные

коэффициенты запишутся в

ввде

 

 

 

 

 

 

 

 

ah

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Л" [ Ж ' в , п +

 

 

sin(ah* arcty а>г\е -~ г

 

 

 

 

 

ah

 

 

 

 

 

 

 

 

 

 

 

 

sin ah,

 

 

 

 

 

 

~T~ |

 

(Ю)

C ^ ~ ^ ( h - 2 a l)^ \(^ ~ a>Tz '}swti)h+ 2ctrcostnh j e

 

 

*

-[(Mr+w

 

 

 

 

]'

 

 

Полученные коэффициенты позволяют

при любом

//,

когда вход­

ное воздействие

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

изменяется

согласно

(7 ),

рассчитывать на основе (4) точные значения реакции объекта в диск­

ретные моменты времени. Погрешность расчёта переходного процесса

в динамическом объекте появляется только при отклонении закона

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

h

от его

линейной апроксимации. Таким образом, точность р асч ета

по данно­

му методу зависит не от параметров передаточной функции

объекта,

а 0Т точности линейной апроксимации входного воздействия

на рас­

четном шаге.*

 

 

176

. Расчетные модели интегродифференциальных объектов второго порядка

'и/(р )

Постоянные

коэффициенты. А,В, С,

вычислителя

для

данных

W(p)

передаточных

 

 

 

 

 

 

 

 

 

функций

 

 

 

 

 

 

 

 

 

д= IwTsina^

 

 

 

 

s)r>-

 

arctff ®Г) \ е~АГ*

В =

с Л/Тsirr

 

 

Щ р *п

с- ж

\

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Грг+2Тр+ (7+ДгГ2)

 

 

 

 

 

 

 

 

 

 

-р/р .

. .

2Т-(Гг-Ь)0**?Г*)

 

 

кТг \ уЛИ-ЛГ+ГГг -/*-йГ,*гГ*)г* 4 ‘Г*(АГ-*»Т7-ЪГ)я'

 

 

M~hT}

[

 

 

 

r-fa(t+o*r*)

 

 

е

*П(«*+/)+

(t+a*r*)*

 

 

 

 

/+ й>Г2

 

] '

*

~ arct9 "Щ-АГЧ- тт, - г * - ц л *г*

-*а»*9(-*т)

 

 

НтгР+1)

4= [ Г/ЛSin”,/Г + COS*/Г ] ;

 

С**11 f a *А)-Ггсо54 - /*,* £ ] ;

 

 

 

Т2р2+1

B = -T s,nT ’

 

 

 

 

 

Г,'A) MS

 

 

Гу

in Л /г]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ИГтРЧ )

. t.\h(T*h) „-Щ

П Ы ) П-Ш \

г. ягз

Гш - Ь ) - л /%

 

Г(Г,-Г) - л /г ................,7 ■

 

А~Ыт3- п е

 

 

щ - п *

\’

 

7 ^ Г

е

~сг-г3) е +<ь+П-г-Щ

 

(Гр+/)($р*Г)

п

[ Ъ

_ _ £ / _

А /п .

w t o - m

w

ш _ Ог&0>+Г3) -Щ

+<**+»!

1

 

 

в

[h(r5- r ) e

 

Щ Т )*

J ’ 3

щ

1 сггЬ).

 

 

c r -t)

 

/ ; j

к(.Гьр+1)(Т2р*Г)

* -

/ / > [ '♦ > - г * - * "

] ;

 

М

 

 

(

 

 

 

?

 

'

р(Гр*П

' ~ f ( » * " ) >

 

' ^ [ ( т

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A ( fr + n

.

А

 

 

3--Г,

 

С -к

(А г/зт - - f f

)■

 

J)*Hft2/ e r +

)

 

Гр* ..

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ш рТГ)---------

A - e h/ r

;

В*0-,

 

c^B -\-i-A (rr T)*h\-t

2 =2 -[(т-А)(ГГ Г)+РА]

 

Гр

f f

~ ~ 7 Г п

 

т у -

 

с , к {? & * -)■ ,

 

 

j > - k ( w / ; - y

* J 4

L -

 

 

 

 

Выполнив аналогичный расчет и преобразования для всех корней характеристического уравнения объекта моделирования, получим соот­ ветствующие значения постоянных коэффициентов А, В, С, В , расчетные

выражения для которых сведены в таблицу. Для сравнения в таблице кроме других пяти типов моделей второго порядка приведены две мо­ дели первого порядка, для которых так же справедливо рекуррентное

соотношение -(4 ).

 

А. В,

С, J),

Полученные расчетные выражения для

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

позволяют, зная аналоговые параметры непрерывных объектов,

такие

как к - коэффициент усиления и К

постоянные

времени, пе­

реходить к точным дискретным моделям вида (4 ) . При изменении физи­ ческих параметров объекта следует, согласно таблице, пересчитать соответствующие параметры дискретной модели. Расчет коэффициентов /ft в, С7 в осуществляется на этапе настройки модели на задачу исследования.

Расчет по рекуррентному соотношению (4) справедлив для многих конкретных объектов первого и второго порядка, что позволяет реали­

зовать эффективный интегродифференциальный вычислитель с адаптив­ ной параметрической настройкой.

Разработанная модель и метод расчета коэффициентов рекуррент­

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

• I. Баранов Г.Л. Раочет на ЦВМ нелинейных динамических систем непосредственно по структурной охеме модели. - Изв.АН СССР. Энер­

гетика и транспорт, » 6, 1969, с . 71-79.

 

преобраэова-

 

2. Деч Г. Руководство п р акти ческо м у применению

 

3 . Мартыненко В.С. Операционное исчисление. - Киев : Высш.

школа, 1973. - 267 с .

 

 

УДК 621.313.321

 

 

Гв.П.Коваленко ,| Н.А.Куличихина

 

 

 

СПЛАЙН-МЕТОДЫ В ЗАДАЧАХ ВДЕИТИФИКАЦИИ

 

 

 

ЭЛЕКТРИЧЕСКИХ МАШИН•

 

 

Рассматривается задача определения параметров дифференциаль­

ного

оператора Н(р) взаимосвязи между векторами

входного воздействия

Х(р)

и реакции объекта Y(p) , Предполагается,

что для

искомой функ-

178

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

Для решения задачи идентификации предложено много различных методов / 2 , 3 , б7, причем их большое разнообразие объясняется некор­ ректностью задачи, как следствие строгого несоответствия аналити­ ческих свойств функции ff(p)€ С " и экспериментально определяемых .

динамических характеристик,

для которых,

в лучшем случае, справед­

ливо

утверждение

об их

принадлежности к

классу непрерывных функг

ций,

т . е . /КО,

Н ( ъ ) € С ° .

 

 

 

Задача

определения

т

неизвестных коэффициентов функции ^

сводится к

решению матричного уравнения

т -го порядка

л* -в.

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

Учитывая плохую обусловленность определителя Гильберта / 7 / , попытка определить из (I ) более, пяти-шести неизвестных коэффици­ ентов операторной функции Н(р) нереалистична. Удовлетворительные

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

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

Подходящим для

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

■и аппроксимирующие

кубические

сплайн-функции / 1 , 4 / .

Опыт работы

со сплайн-функциями

показал,

что условие

/ Б /

 

 

P j'-M /fy/f: Mf/H

<5,

(2).

где М-, - величина второй производной интерполирующей сплайн-фувк- ции в /'-м узле сетки определения динамической характеристики „ яв­ ляется хорошим критерием для оценки однородности исходной инфар>- мации. Причем нарушение этого условия указывает либо на грубую ошибку измерения динамической характеристики, либо на ошибку перформа'ции исходных данных.

179

При измерении динамическая характеристика обычно задается

в виде простейшей, но наиболее распространенной регулярной (пери­ одической) дискретной последовательности. Однако известно доста­

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

ностью подавления шумов измерения. Среди них, по-ввдиыому, наилуч­ шими являются’ при анализе частотных последовательностей тригономет­ рические ряда Фурье, а для временных последовательностей - ряды Пуассона /5 7 , при вычислении которых реализуется близкая к опти­ мальной методика фильтрации исходных данных.

В работе / 3 / показано, что между коэффициентами спектров

Фурье или Пуассона и искомыми коэффициентами операторной функции

НСр) существует простая й однозначная

взаим освязь. Так,

в

частно­

сти, из теоремы о смещении изображения

непосредственно

следует,что

коэффициенты спектра Пуассона

отличаются только знаком

от

коэффи­

циентов Тейлора -

Маклорена разложения

операторной функции

Я(р)

в степенной ряд в

точке

/>**

 

 

 

 

 

f K p ) * \ e (Р

)

=

~§~(

\

e * * t* q Ct)c(t =

 

 

- 'L c -P) t Pk,

где Рк - коэффициенты спектра Пуассона.

Для вычисления коэффициентов спектров Фурье или Пуассона раз­ работано достаточно много'методов /3 ,7 7 , среди которых для полу­ чения их неомещеныых оценок целесообразно использовать метод (по идее аналогичной методу Филона), использующий кубическую сплайнфункцию для гладкого восполнения дискретной последовательности динамической характеристики объекта.

Вычисление моментов Пуассона сводится при этом к простой реккурентной процедуре:

■V *

4 * * и *

+

+ * /'* * > .

.

\

Ы н

nr

т-1 ,

Ят s J

е

dw,

где Of, bj. Cj - коэффициенты интерполирующей сплайн-функции в /-м узле сетки определения Ш )

180