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

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

.pdf
Скачиваний:
5
Добавлен:
23.10.2023
Размер:
5.03 Mб
Скачать

так-как

..min . win

 

175

МВт и не может бить меньше

(этим учи­

 

5.

тываются ограничения).

 

 

 

 

 

 

 

 

После того, как таблица составлена полностью, определяем рас­

ход sea sa ,

соответствующий оптимальному распределению нагрузок.

В наиеы случае

это

|s (N)

= 1030,1 Гкал.

 

 

 

 

 

Этому расходу соответствует нагрузка

= 50 МВт и

Q

(50)=

= I2X, Гкал. Расход тепла двумя первыми агрегатами при этом со­

ставляет

f2(N)=f3(.H) -

Q3 C5 0 )=

Ю 30.1

- 121 ,0

=909,1

Гкал .

Из таблицы этому значению

| (N)

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

оптимальная

нагрузка

М^ =

175

МВт. Функция

^(N) = ^{N) - Qj (I75)

=909,1

-

- 346,3 =

5 6 2 ,8

Гкал,

чему соответствует

оптима^кая

нагрузка

= 300 МВт. Таким образом, оптимальный план распределения

на­

грузки 525

МВт найден.

 

 

 

 

 

 

 

 

 

 

агрегатами

 

Однако мы решили

не только поставленную задачу.

По таблице м о ­

гут быть найдены

оптимальные

планы распределения

нагрузки во

всем диапазоне от 225 до 525

МВт, так как вместо

одной задачи

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

ются оптимальные нагрузки

 

,

Wj ,

Результаты опреде­

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

М

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

на рис. 4

-1 0 .

 

 

 

 

 

 

О п т и м а л ь н о е

к о л и ч е с т в о

з а п а с н ы х

ч а с т е й .

Задача формулируется

следующим образом: система ■

содержат

ГЦ

элементов типа I

, предназначенных для работы

в течение

t

часов ; один

элемент

i

-ого

типа

имеет интен-

90

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

X-L и

стоит CL рублей.

Необходимо определить

оптимальное количе

гво

запасных элементов

каждого типа (

) ,

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

достаточности

запас­

ных элементов d.

при

заданных ограничениях на их стоимость.Под

вероятностью достаточности запасных частей понимается то , что за

время

{

с заданной

вероятностью

d. запаса

потребуется не

бо­

л ее ,

чем

заготовлено.

 

 

 

 

 

 

Рассмотрим систему

- часть

низкого давления

(ЧНД)

турбины

 

К -300-240

; Ч1Щ - трехпоточная,

в

каждом потоке

по 5

ступеней.

Определим оптимальное количество запасных сопловых и рабочих ло­

паток (

рядов

лопаток,так как при отказе в ряду даже

одной лопат­

ки, как

правило, перелопачивается весь ряд)

для обеспечения ве ­

роятности R

= 0,999 за

{ = 6000 часов.

 

 

В табл.

4-8 приведены

частоты отказов (

l L

) ; средний

расход лопаток за время 6000 часов и относительная стоимость(С )

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

стоимости сопловых лопаток

первой

ступени.

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 4 -8

:

 

Р я д ы

-

 

!4 4

; И-ерД 1

Ci

пп

 

 

I

 

Сопловые лопатки

I

ст .

2 .2

0,0396

I

2

 

Рабочие

лопатки

I

ст.

2,8

0,0504

3

3

 

Сопловые

лопатки

2

ст .

2 ,2

0,0 3 9 6

1 .2

4

 

Рабочие

лопатки

2

ст .

5

0 ,0 9

3

5

 

Сошговые лопатки

3

ст.

2 .2

0,0396

1 ,2

6

 

Рабочие лопатки

3

ст .

6

0 ,108

4

7

 

Сопловые лопатки

4

ст.

2,8

0,0504

2

8

 

Рабочие лопатки

4

ст .

7

0 ,1 2 6

6

9

 

Сопловые лопатки

5

ст .

2 .8

0,0504

5

10

 

Рабочие

лопатки

5

ст .

10

0 ,1 8

10

В таблице значения

1^ приняты близкими к полученным в

испыта­

ниях турбин

мощностью 150 МВт [ 5 6 ] ; относительная стоимость ря -

дов лопаток

принята по

[5 7 J для турбин К -300-240 ЛМЗ.

 

Как показал анализ,

с вероятностью,большей 0 ,9 8 , закон

рас -

пределения вероятности безотказной работы лопаток является экспо­

ненциальным

[5 8 J .

 

 

Поэтому

средний расход рядов лопаток определен по формуле .

 

N-cp.i = ^t A it

= 3 .6000 l l f

(4 -7 7 )

 

 

 

91

так как в 'W jjo

три ряда лопаток каждого типа_(_три _потока

пара).

По

показателям

 

и

 

с помощью таблицы

7-2

[59.]

определяется вероятность недостаточности запасных рядов

.

а по

ней необходимая нам вероятность

 

 

. Таким образом,

варьируя Hj ,

можно определить

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

им вероятности

сЦ н затраты в

запас

лопаток определенного типа. Общая вероят -

ность

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

как

произведение вероятностей oIl

 

 

 

 

 

 

 

 

 

 

(9 -7 8 )

• Рассмотрим

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

этапы

процесса вычислений

с

по­

мощью

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

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

последовательностей '■[бО]

-Мажо­

рирующей последовательностью называется такая последовательность

{(Й0,СД

}

, для которой

переход в-состояние

с

более

высокой вероятностью

R

происходит

с минимальными затратами С.

. Такая последовательность есть не что иное, как зависимость

вероятности достаточности запаса от стоимости, полученная

 

при ■

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

 

 

На каждом этг.пе процесса вычислений объединяются попарно от­

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

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

будет по

лучена последим возможная комбинация. В нашем случае ряды лопа­

ток образуют подсистемы: I

- сопловые лопатки первой ступени,

П- рабочие лопатки первой

ступени ;

D!

- сопла второй ступени и

т .д .

 

Эти подсистемы

объединены в

комплексы, как

показано

на

рио.

4—II .

 

 

 

 

 

 

 

 

 

 

 

Irnm.

2 mm.

Зпап.

 

ta m .

5стип.

 

 

Оуаеткм, что таких.объединений можно осуществить ( I

) ,где

i - число подсистем

а общих правил объединения

не

/

 

 

92

Рабочие лопатки

существует.

 

 

 

 

 

 

 

 

 

 

Для первой пары подсистем составляется таблица,

где в

верх­

ней графе

указаны

т‘оличество запасных рядов

2^

,

вероятность

и стоимость

С-!

» получающаяся в результате добавления

одно­

го запасного ряда

в первой подсистеме С в таблице

4 -9 это подси­

стема I ) . Аналогично

заполняется левая часть

таблицы

(по

подси­

стеме

П ) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 4-9

 

 

 

 

 

I . Сопловые ,лопатки

 

 

 

 

 

* 1

:

0

 

 

 

:

2

 

 

:

3

 

d l

:

0 ,9 6 1 6

 

0 ,9 9 9 2 3

0,99999

 

 

1 ,0 0 0 0

Cl

0

 

 

I

 

2

 

 

 

3

 

0

 

0+0

 

 

1-Ю

 

2 Ю

 

 

 

3+0

 

0 ,9 5 2 4

 

0 ,9 1 5 8 0

 

0 , 9 5 1 7 0

0 ,9 5 2 3 9

 

 

0 ,9 5 2 4 0

0

 

0

Со)

 

I

C I )

2

 

 

 

3

 

I

 

0 + 1

 

 

I + I

/■•

2+1

 

 

 

. 3 + 1

 

 

 

 

 

 

 

 

 

0 ,9 9 8 8 1

0 ,9 6 0 4 6

 

/ 6 ,9 9 8 0 4

0 ,9 9 8 8 0

 

 

0 ,9 9 8 8 1

3

 

3

( 2 )

X

4

О )

5

С4)

 

 

6

 

2

 

0+2

 

 

1+2

 

2+2 ^

 

 

 

3+2

 

0 ,9 9 9 8

 

0 ,9 6 1 5 8

 

0 ,9 9 9 2 1

/ 6 , 9 9 9 9 7

 

 

0 ,9 9 9 9 8 .

6 '

 

6

 

 

7

С5) / X

8

 

 

С96 )

С 7)

3

 

0+3

 

 

1+3

 

2+3

 

 

 

3+3

|

1 ,0 0 0 0 0

0 ,9 6 1 6

 

0 ,9 9 9 2 3

0,9 9 9 9 9

 

 

1 ,0 0 0 0

9

 

В

 

.

ю

 

I I

 

 

 

1 2

( 8 )

В каждом квадрате сумма цифр означает общее количество запас­ ных рядов С первая - сопловых лопаток, вторая - рабочих); вторая

цифра

вероятность*достаточности запаса

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

произ­

ведение

< *i

С0,9158=

0,9616x0,9524) ;

последняя цифра - затраты

в запас.

 

 

 

v

 

 

 

В таблице строится мажорирующая последовательность для подси­

стем I

и П С показана линией ) .

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

по

квадратам

проверкой наибольшего

увеличения вероятности R =Пс4.-1

на единицу

затрат.

 

 

 

 

 

Аналогично

строится

таблица

(т а б л ,4-10)

для подсистемы

2

С второй ступени ) .

 

 

 

 

 

На следующем этапе

объединяются данные

о стоимости и вероятно-

93

 

 

Hi.

Сопловые лопатки

Таблица 4 -1 0

 

 

 

2 ,

0

 

I

 

2

3

 

0 ,9 6 1 6 0

0 ,9 9 9 2 3

0 ,9999 9

1 , 0 0 0 0 -

 

0

~

1 . 2

 

2 ,4

3 , 6

0

0-Ю

 

1-Ю

 

2+0

з ю

0 , 9 1 7 7

0 ,8 8 2 4 6

0 ,9 1 6 9 9

0 . 9 1 7 6 9

0 , 9 1 7 7 0

0

0

Со)

. 1 , 2

C I )

2 , 4

3 , 6

1 I

0 + 1

 

-

/

2+1

3 >1

 

I+ I

 

ста R , взятые из ранее полученных мажорирующих последовательно­ стей, Такое объединение для подсистем I и 2, дающее мажорирую ~ щую последовательность для подсистемы б приведено, в таблице

а-П.

Далее эта процедура проводится для подсистем 3 , 4 , 5 . 7 ,8 -; объединяя 8 и б, получаем мажорирующую последовательность для всей системы, С помощью окончательной таблицы "обратным ходом"

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

необходимой

вероятности

С в нашем случае R.

= 0, 999 )

достаточ­

ности запаса.

 

 

 

 

 

Результаты

расчетов

приведены в таблице

д -12 . Значение

=

0,9994

удовлетворяет задач е,'так как

и з-за целочисленно-

сти

^

нельзя

добиться

точного значения

R - 0 .9 9 9 .

Значение

вероятности

R

= 0,4876

соответствует случаю отсутствия запас­

ных рядов

лопаток,

 

 

 

'

Как и

в

задаче по оптимальному распределению нагрузки, в .

рассмотренном примере решена целая серия однотипных задач,Опре­

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

иых рядов лопаток разных типов

"во всем диапазоне требуемой ве­

роятности достаточности запаса

от 0,4076 до 6 ,9 9 9 4 .

 

 

 

 

 

 

 

 

 

 

Таблица 4 - Ц

 

 

 

 

Подсистема I

 

 

 

 

 

k L

: 0 , 91580

: 0 , 4 5 1 7 0

: 0 ,9 6 0 4 6

: 0 ,9 9 8 0 4

: 0 ,9 9 8 8 0

: 0 ,9 9 9 2 1

: 0

, 99997

: 0

, 99998 : 1 ,0 0 0 0

С (

0

1

3

4

5

7

8

 

9

12

Подсистема

4 .2

0,99912

7 ,2

0.99998

8 .4 .

0,99999

I I .4

1,00000

12,6

*

8 ,2

(6 )

0,99716 \ 0,99792(8)

11,2

(7 )

12,2

 

 

 

 

 

 

-

■ I---------------------

0,99909

0,99985

 

 

 

0,99868

 

 

 

13,4 (9)

15, 4 (10)

16,4

( И )

 

 

 

 

 

 

0.99996 |

0,99997

 

 

 

 

 

19,4

(1 2 )

Р П , 4 CI3)

 

 

 

 

 

 

 

0,99998 I

1,0000

 

 

 

 

 

 

21, 6 (14)

(1 5 )

 

 

 

 

 

 

 

 

Таблица

u - I 2

 

 

Необходимее

количество

запасных

рядов

 

 

 

& ;

с

1 I

ступ. : 2

ступ.

: 3

ступ.

:

4 ступ.

:

5

ступ.

 

 

.Сопл. :Р.л:Сопл. *.£>.Л. •.Сопл. :Р .л . :С о п л .:Р .л .: Сопл: Р-Л

0,9994

8 2 ,8

2

2

2

2

2 .

2

2

2

2

3

0;9929

63,8

2

2

2

2

2 '

I

2

2

I

2

0,9840

4 9 ,4

I

I

I

2

I

I

I

I

[

2

0,6553

3 2 ,2

I

I

0

0

I .

I

I

I

I

I

0,7210

20,4

I

0

I

0

I

I

I

I

I

0

0,(3044

8 ,2

I

0

0

0

I

I

I

0

0

0

0,5266

2 .2

I

0

0

0

I

0

0

0

0

0

0,4876

0 ,0

0

0

о .

0

0

0

0

0

0

0

Общие вопросы динамического программирования подробно рассмот­ рены в [4 ,6 ,1 4 ,2 6 ,2 9 ] , применительно к энергетике в [3 5 ,6 1 ] ; воп­ росы дискретного динамического программирования изложены в jj52j .

§ 4 -5 . Линейное планирование

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

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

При втором способе решения задач на первом этапе получают ус­ ловно - оптимальный план, который обеспечивает оптимум целевой функции, но не является допустимым по ограничениям. Он становится допустимым лишь после определенных преобразований. Сюда относятся: методы решения транспортной задачи, метод разрешающих слагаемых А.Л. Лурье, метод дифференциальных рент А,Л. Брудно, венгерский метод [ 4 2 ] .

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

96

еле л ойчтельное у.т шение опорного плана по

определенной схеме.

Эта схема вычислении следующая:

 

1. В опорном плане выделяются базисные

переменные, т .е . пере­

менные, входящие в этот план. Базисные переменные опорного плана

выражаются через небазисные, не входящие в

план ;

2. Выражается значение целевой функции

и через небазисные

переменные.

 

3 . Выбирается та из небазисных переменных, введение которой в план способно, улучшить значение U .

'4 . Определяется, какая из базисных переменных должна быть при

этом исключена из плана и сделана небазисной.

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

6. Все остальные

базисные переменные и значение функции Ц

выражаются через новые небазисные переменные.

7 . Повторяются операции, указанные в

пунктах 3 , 4 , 5 , б.

Если на каком-то

этапе оказывается,

что введение в план любой

из небазисных переменных не улучшает U. ? то последний план оказы­ вается оптимальным. Доказано, что если следовать приведенной схе­ ме, то через конечное ( и не слишком большое ) число шагов можно

прийти к оптимальному плану [_4о] .

Рассмотрим использование симплекс-метода к решению задачи

определения

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

с о с т а в а

с м е с и .

Имеется три сорта топлива с различной теплотворной способно­

стью, зольностью и стоимостью

С табл. 4 -13 ) .

 

 

 

 

 

 

 

 

Таблица 4-13

Сорт топлива

:

Зольность

:

Тепло творит

:

Дева

 

.

К Г / К Г '

:

Qp , ккал/кг

:

руб/т

 

 

 

I

 

0,0 В

 

8000

 

12

• П

 

0 ,1 6

 

600Q

 

8

Ш

 

0 ,3 2

 

4000

 

4

Требуется определить, сколько и какого сорта топлива необхо­

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

5000 ккал/кг,

общее количество

золы не превышало 0 ,1 4

§[Н гоы 1ва.

а затраты на топливо были минимальны.

 

Обозначим

х, » ЭСХ, ,

-

неизвестное количество

к г топлива

97

каждого сорта. Тогда задача сводится к минимизации стоимости топ­ лива, глиной

 

U.=.№I,+6Xl +l»Xs

(X L*0)

(4 -7 9 )

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

 

 

 

8000

X , + 6000 х 2 + 4000 х ъ - 5000

(4 -8 0 )

. 0 ,0 8 X, + 0 ,1 6 Х г + 0 ,3 2 Х ь « 0 ,1 4

(4 -8 1 )

Прежде всего

задачу необходимо свести

к каноническому

виду

- неравенства заменить равенствами. В нашем случае для преобразо­

вания (Л-8 1 )

в

равенство можно ввести

дополнительную переменную

Хц

-

количество фиктивного

топлива с

теплотворностью

и

стоимо­

стью, равными нулю,но с зольностью 1%,

Тогда условие (4 -0 1 )

примет

вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

_8Х , + 1 6 Xj. +32 Хъ

+ Х 6 =

14

 

 

 

(4-82)'

 

При разделении переменных на базисные и небазисные (свобод­

ные)

придерживаются общепринятой

записи:

в системе

ограничений

-базисные переменные стоят всегда в левой части, небаэисные - в

правой. Примем за базисные переменные х,

и Хц '.поэтому

выраже­

ния (4 -8 0 ) и

(4 -8 2 ) представим следующим образом

 

 

 

 

 

 

 

 

 

х { =

1

-

й х ,

 

-

 

it зсъ>

 

 

 

О -ез)

 

 

 

 

 

 

8

 

8

 

8

 

 

 

 

 

 

 

 

 

 

 

Хц = 1 4 - е - 16Хг - 32 X j

 

 

 

(4 -8 4 )

 

г эдставив

в

(4 -8 4 )

выражение

(4 - 8 3 ),

получим

 

 

 

 

 

 

 

 

 

= 9-IO

- 28

 

 

 

 

 

 

(4 -8 5 )

 

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

ные

переменные

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и М 2 ( | - | х г | х ь) + В Х г + ^ Х 5 = 7 ,5 - Х 2 - 2 х 5

 

(4 -8 6 )

Примем в качестве опорного плана:

X, з Х 5 = 0

; из

 

(4 -8 3 )

следует,

что X, = g

;

из

(4 -8 5 )

-

 

=9 и из

(4 г 8б)

-

за­

траты

U

= 7 ,5

руб/т.

 

 

 

 

 

 

 

 

 

 

 

 

Из

(4 -8 6 )

видно, что введение

в

план

ведет к

снижению

затрат;

примем

поэтому

 

за

базисную переменную.

Из

(4 -8 3 )

следует,

что

 

должна быть исключена лэ плана,

то есть

стать

небазисной переменной. В соответствии с пунктом 5 расчетной схе­ мы из (4 -8 3 ) получаем

98

 

 

 

 

 

 

 

 

 

 

 

 

(4 -8 7 )

Подставляя

(4 -8 7 )

в

(4 -8 5 ) и (4 -8 6 ) , получил

 

 

 

 

 

V

9 - W ( |

- | х , - | х ъ) - . ? 8 Х з = | ^ х г

Й Х 5 ,

 

(it_88)

И

» 7 ,5 -(| -

| Х , -

|XS) - 2 X * * 6 ,6 6 6

+ | Х Г

| Х 5

 

(4 - 8 9 )

Полученный план выглядит

следующим образом:

X , = X s = О

;

Х г = 5/6

;

Х[,

=

4/6 и

U = 40/6= 6,666 руб/т. Таким образом,

движение осуществляется

к оптимальному плану,

т .к .

U

уменьша­

ется.

 

 

 

 

 

 

 

 

 

 

 

 

Введем

вместо

 

 

в

план переменную

X j

.И з

(4 -8 8 )

полу­

чаем

 

 

 

 

 

 

 

 

 

 

 

 

Х 5 = 1/32

+ 5/8

- 3/64X1,

 

 

 

 

(4 -9 0 )

Тогда

по

(4 -8 7 )

 

 

 

 

 

 

 

 

 

а по (4 -8 9 )

 

 

 

 

 

 

 

 

 

 

 

Из (4 -9 1 )

тут же

следует, что введение в

план

X,

и

не

уменьшает

значение

U

,

а , следовательно, полученный план

Х , = 2 ц * 0

 

;

Х-,

=

39/46= 0 ,8 1 2 5 -;

XJ=

1/32

= 0,0313 ;

U = 6,625 руб/т является оптимальным.

 

 

 

 

 

 

Процесс нахождения оптимального плана значительно упрощает­

ся , если его

свести

к

заполнению стандартных,

так

называемых

симплекстаблиц [4 ] .

Рассмотрим еще одну задачу линейного программирования.Пред­

положим, что имеется

IX

пунктов отправления

Kj ,

, . . .

, в которых

сосредоточены запасы какого-то однородного гру­

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

й, ,'С Ц , , . . . , Qm .

Имеется,

кроме того,

ft

пунктов

назначения

С, , Ct , . . . ,

, которым

необходимо

подать &

,

 

единиц груза. Требуется

составить такой

плав перевозок, при котором все

заявки

на грузы

были бы выполнены, а стоимость всех

перевозок была бы минималь­

ной.

 

 

 

 

 

 

 

99

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