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

книги из ГПНТБ / Ильинский, Д. Я. Обоснование решений при проектировании и эксплуатации машин и линий легкой промышленности учебное пособие

.pdf
Скачиваний:
4
Добавлен:
19.10.2023
Размер:
3.38 Mб
Скачать

Рис. 7. Схема для нахождения оптимальной загрузка машин

Проведение

модернизации

целесообразно, если

максималь­

но возможная

прибыль,

получаемая

при

II варианте, боль­

ше максимально возможной прибыли, получаемой при I варианте.

Условия задачи

сведены в

табл.

10.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т а б л и д а

10

 

 

 

 

 

 

 

 

Продукция* '

 

 

Величина

 

 

 

А вида'

В' вида

С вида

 

 

 

 

 

 

11

I

. п

п

Заданные постоянные

 

 

 

 

 

 

 

 

 

Производительность,

тысяч штук

 

 

 

 

 

 

в месяц (тыс. шт./мес.), цеха:

 

 

 

 

 

 

 

штамповочного

 

 

 

25

20

35

30

 

12

механического

 

 

 

33

30

16,7

10

10

сборочного

 

 

 

 

22,5

20

15

12

8

Прибыль,

получаемая от одного

15

15

 

 

14

14

изделия,

коп.

 

 

 

 

12

12

Ограничения

 

 

 

 

 

 

 

 

 

 

 

Общая длительность

работы

Т

1

1

1

1

1

1

Искомые (управляющие) переменные

 

 

 

 

 

 

Количество

изделий

определен­

Л

 

 

 

 

 

ного вида, выпускаемое в

месяц

 

 

О

 

<3

31

Частное от деления количества продукции Х\, х2 и х3 (в штуках)

соответствующего вида на производительность Q шт' (см. табл. 10)

мес.

определит длительность работы (в долях месяца) по выпуску про­ дукции этого вида.

Рассмотрим вариант I. Ограничения по ресурсу времени:

 

 

 

5 1

 

 

(I)

 

 

 

35-Юз

^

 

 

 

 

 

 

Л*2

^

1.

 

(И)

 

 

 

16,7-li з 44

 

 

 

 

< 1;

 

 

 

 

(III)

 

 

 

- —

 

<

1-

(IV)

 

 

 

 

15-Юз

 

.

 

Условия неотрицательности x'i>0; л'о>-0.

 

от выпуска л-j штук

Целевая

функция — суммарная прибыль

продукции А

и хъ штук продукции В.

 

 

 

 

 

 

F= 15 .v'i+ 12 х2 ■* max.

 

 

Строим

многоугольник

допустимых

 

значений

переменных

(рис.

8 ) и находим точку М (20,3; 6,4),

принадлежащую этой об­

ласти

и лежащую на прямой

12 х2.

 

 

 

Прямая

F =

15 х'! +

 

определяющую опти­

F, проходящая через точку М,

 

мальный объем каждого вида продукции А и В, соответствует мак­ симально возможной прибыли.

F = 15 • 20,3 • 1 0 3 + 1 2 • 6,4 -,1 0 3 = 381,3 ■1 0 3 коп. = 3813 руб.

Рассмотрим вариант II.

В отличие от варианта I количество искомых величин в вариан­ те II увеличилось, так как необходимо найти максимальную при­

быль от выпуска а'! штук продукции А,

х2 штук продукции В и хэ

штук продукции С.

имеем:

 

 

 

 

 

 

 

 

Аналогично варианту I

 

 

 

 

 

 

 

 

— +

Хо

 

+

 

 

 

<

1 ;

(D

30-10-3

12-103

20 -Юз

 

 

 

 

X

-1- • Л'о

f

л'3

 

< 1 ;

30- 10»

 

10-103

 

1

10-103

 

 

 

 

2 J-

 

1 ;

 

 

 

 

 

 

 

(III)

 

 

 

 

 

 

 

 

 

 

 

 

Л'г

 

 

 

 

 

 

 

(IV)

 

 

12-103 ч< 1 ;

 

 

 

 

 

 

 

 

 

А'8

 

<

1 ;

(V)

 

 

 

 

 

8 -Юз

 

 

 

 

 

 

л-2,

л'з > 0.

 

 

 

 

32

xz-m'3

Рис. 8. Схема для нахождения оптимального ассорти­

мента продукции (до реконструкции)

Неравенства (III) и (IV) из-за неотрицательности пере­ менных оказываются более слабыми ограничениями, чем соответственно (I) и (II), поэ­ тому их можно отбросить.

Целевая функция аналогич­

на целевой функции

по вари­

анту I

 

12 л'2 +

F = 15 .Vi +

+

14 л'з - max.

(VI)

В

трехмерном

простран­

стве

Ох]Х2хз

многогранник

OABCDEEG определит об­ ласть допустимых значений переменных (рис. 9).

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

Рис. 9. Схема для нахождения оп­ тимального ассортимента продукции . (после реконструкции)

3.3 -

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

Обозначим через х4, х5 и х6 разности между правыми и левыми ■частями соответственно неравенств (I), (II) и (V). Переменные х4, л'5 и Л'6 выражают неиспользованные мощности соответственно по штамповке, механической обработке и сборке продукции С. В результате получаем возможность заменить эти неравенства систе­ мой равенств:

■v’

I

'V2

1

Х *

 

I Г

1 •

< П

20-10»

1 30-10»

'

12-10»

4

 

 

 

* '

4-

■**

+

'Г;|

+

х 5

1;

( I V )

30-103

 

10-10»

 

10-Юз

 

5

 

 

 

 

 

■ П 3

+

=

1 •

(V ')

 

 

 

 

8-10»

 

ь

 

 

Очевидное решение xi= x2==x3=X4 —xs — xe = Fc=0 соответствует вершине многогранника OABCDEFG, лежащей в начале коорди­ нат.

Выберем переменную, которая входит в целевую функцию F (VI) с наибольшим положительным коэффициентом xb

Сравним величины отношений свободного члена уравнений (Г), (1 Г) и (У') к коэффициенту при переменной ху в каждом из ука­ занных уравнений и выберем наименьшее положительное из инк:

1 / ( 2 0 -1 0 3) - 1 = 2 0 -1 0 »;

1/(30-103) - 1= 3 0 - 103;

1 /о = °°.

Наименьшее отношение 2 0 1 0 3 соответствует уравнению (I'). Составим новую систему уравнений:

*1 + -

jc, + -

*з + 20 -10-3 х 4 = 20-10";

 

(Г)

3

3

 

 

 

 

 

•X,

90-109

f x 4+ x 6 =

- i;

( i n

90-109

и

О

 

 

 

■/гг +

-^б — ^ •

 

(V')

 

 

8-109

 

 

 

Первое уравнение системы получено умножением уравнения (Г) на 2 0 - 1 0 3; второе- - сложением уравнения (I"), предварительно умноженного на (—30- 103)-1, с уравнением (1Г)'; у|равнение (V ) фигурировало выше.

Целевая функция

2 х 2 — 11 л-з-ЗО-Ю* х 4 = F — 30-Ю4

(VIх)

получена сложением уравнения (I"), предва1рительно умноженного на (—15), с выражением (VI).

Новое,

улучшенное

решение,

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

точке А (см.

рис. 9)

 

 

 

 

 

* 2 = х 3 =

= 0 ; лу = 2 0 1 0 3; х 5 — — .

 

 

 

3

При этом Fa = 3000 руб.

 

 

Сравним величины

отношений

свободного члена уравнений

(I"), (II")

и (V') к коэффициенту при переменной,

входящей в це­

левую функцию с наибольшим лоложительным коэффициентом. В целевую функцию. (VI') с наибольшим положительным коэффи­ циентом входит переменная х2.

..Сравним величины отношений свободного члена уравнений (I"),

(II")

и (V') к коэффициенту при переменном х2 в этом же уравне­

нии

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 30-103,

 

 

 

 

 

 

 

 

 

 

 

 

 

2/3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

30-10»

 

 

 

 

 

 

 

 

 

 

 

 

90 - Ю»

7

 

 

 

 

 

 

 

 

 

 

 

 

1 / 0 =

©о.

 

 

 

 

 

 

 

 

, ,

 

 

 

 

 

 

 

 

 

30-108 .

 

 

 

Наименьшее положительное отношение— ^- соответствует

уравнению (II").

 

 

 

 

 

 

 

 

 

 

 

 

 

Составим новую систему уравнений:

 

 

 

 

 

 

 

 

 

 

 

 

60-108

90-1Q8

 

30-,108

»

 

(1П

 

х 3 + тт х з — "" 7- х 4 -jr — Г— х-а

 

 

 

 

, 9

 

 

 

 

 

 

 

 

/

 

 

 

 

.

170-Юз

60-103

A'i =

120- Юз

 

( П

 

л'1 +

— х3 н------

----- X j -|------7

 

----

----

 

 

 

 

 

 

8

+ X6 — 1 -

 

 

 

 

 

 

(V0

 

 

 

 

 

-Юз

 

 

 

 

 

 

 

 

 

Первое уравнение системы получено

 

умножением

уравнения

(II")

на 9°-' 10- .

второе — сложением уравнения

(I"'),

предвари-

 

7

I

 

на

(—2/3),

с уравнением

(II');

уравнение

тельно умноженного

(V0

фигурировало выше.

 

 

,

 

 

 

 

 

 

 

 

Целевая функция

 

 

 

 

 

 

 

 

 

 

 

 

 

85

 

22-1 0-'

180-ИГ

х 5

— 308 • 1 03

(VI")

 

7

 

 

7

 

 

7

 

 

 

 

 

 

 

 

получена сложением уравнения (II"), предварительно умноженно­ го на (—2 ), с выражением (УГ).

Новое, улучшенное решение соответствует точке В (см, рис. 9):

х

—х4 — х5

= 0;

 

120-103

17 000;

 

7

 

 

*2

30-1Q3 =

4300.

 

7

 

При этом Fb =3080 руб.

В точке В многогранника допустимых значений OABCDEFG целевая функция достигает максимума, так как в выражении (VI") все коэффициенты отрицательны и переход от вершины В к лю­ бой соседней вершине G или С (в вершине А величина целевой функции уже определялась) только уменьшит значение целевой функции F.

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

Сопоставляя варианты I и II видим, что максимальная прибыль при варианте I больше максимальной прибыли при варианте II (так как 3813>3080). Следовательно, задача расширения ассорти­ мента выпуска продукции должна решаться не частичной модерни­ зацией производства, выпускающего продукцию А и В, а, например применением специализированного оборудования, предназначен­ ного для выпуска продукции С.

Рассмотренные задачи содержали всего две или три управля­ ющие переменные, поэтому легко решались вручную. Когда чис­ ло переменных в задачах значительно (:в частности, задачи поряд­ ка 50x100, т. е. задачи, содержащие 50 ограничительных условий, в которых фигурирует 10 0 переменных, обычно считаются задачами среднего размера), применяются универсальные и специализиро­ ванные ЭВМ.

Например, в результате решения одной задачи такой величины на универсальной ЭВМ «Минск-22» (машинное время решения составило около пяти минут) был создан оптимальный план рас­ кроя рулонной трансформаторной стали для сердечника трансфор­ матора. При этом раскройный коэффициент (отношение суммарной площади выкраиваемых деталей к суммарной площади исходного материала), характеризующий оптимальный план раскроя, вырос с 91 до 97%.

Использование специализированной ЭВМ — электронной маши­ ны ЭМРТ, предназначенной для определения-, оптимального пла­ на раскроя настила ткани (по критериям либо величины концевого остатка, либо количества длин настилов, на которые раскраивается кусок ткани и др.), по данным швейной фабрики имени СмирноваЛасточкина (Киев) обеспечивает экономию 11 тыс. руб. в год.

36

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

НЕЛИНЕЙНОЕ ПРОГРАММИРОВАНИЕ

Задачи оптимизации, которые не укладываются в модель линей­

ного программирования,

относят к

нелинейным (в том числе и те

задачи, которые имеют

линейную

целевую функцию и линейные

ограничения, но на управляющие переменные их наложены условия целочисленности),

Эта область [.математического программирования

разработана

менее других.

значения ско­

З а д а ч а I. Требуется определить оптимальные

ростей элементов машины.

 

Искомые переменные: скорости х\ и х2 рабочих органов соответ­ ственно I и II, функционирующих последовательно.

Ограничения, накладываемые на искомые переменные:

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

зов прямо пропорциональна скорости)

 

2хг + Зл:, < 6 ;

(I)

из условий выполнения заданной циклограммы

 

-*Ч > 0,5;

(II)

* 2 >0,5;

(III)

из конструктивных рекомендаций

 

 

(IV)

х 2 .< 2.

(V)

В качестве целевой функции в данном случае

принимается

функция

 

F = 2хг х I + х а— min,

(VI)

которой прямо пропорциональны приведенные затраты на изготов­

ление единицы продукции.

Построим в координатах х\Ох2 многоугольник ACDF, являю­ щийся областью допустимых значений искомых переменных Xi и х2

(рис. 1 0 ).

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

бол

 

х 2= х] 2х 2F,

(VII)

ось симметрии которых отстоит от оси абсцисс на расстоянии «I»,

37

равным

где а п b — коэффициенты трехчлена у = ах2 + Ьх + с.

Оптимальное решение определяется координатами точки Мх, в которой линия CD касается той из семейства парабол, которая соот­ ветствует искомому минимальному значению целевой функции. Найдем координаты точки М\.

Рис. 10. Схема для нахождения оптимальных режимов работы рабо­ чих органов машин

Из уравнения линии CD (см. выражение I), записанной в виде

■**=—-1*1 — 6 ,

(VIII)

величина тангенса угла наклона линии

CD к оси Ох\ определит­

ся как— 2/з.

(VII)

С другой стороны, продифференцировав уравнение

^ - = 2 ^ — 2,

(IX)

дхх

 

получаем угол наклона линии CD, являющейся касательной к па­ раболе

х2 = х\ 2хг F .

Всоответствии с выражениями (VIII) и (IX) имеем:

откуда Я]'=0,66, а из уравнения (I) :лс2—1,55.

38

В точке М, (0,66; 1,55) целевая функция составляет:

F = 2 - 0 ,6 6 - 0,662 + 1 >55 = 2>44 = ^ ш(п.

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

2*2 + 3*2 ч< 6 ,

(Г).

то область допустимых значений искомых переменных х\ и х2 изо­ бразится в виде многоугольника ABEF (см. рис. 10).

Оптимальное решение, соответствует точке М2. В этой точке касательная к кривой BE, имеющей уравнение

2 * 2 + 3*2 = 6 ,

совпадает с касательной к тому семейству парабол, которое соот­ ветствует искомому минимальному значению целевой функции. Кроме того, точка М2 является общей для кривой

2 А. 2 —|—3*2 = 6

и для указанной параболы. Найдем координаты точки М2. Приравняем производные от уравнения (К), дифференцируе­

мого как неявная функция, и от целевой функции (IX).

- £2 = 2‘ х 1- 2 . (X)

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

темы уравнений

 

— ~ = 2*i — 2

(XI)

2 * 2 + 3*2 = 6 ,

 

откуда *1=0,80; *2=1,26. В точке М2 (0,80; 1,26) целевая функция:

составляет

Е = 2 • 0,80—0,802 + 1,26=2,22=Fmln.

З а д а ч а 2 (метод поочередного изменения величины управ­ ляющих переменных или метод Гуасса-Зайделя). Определить оп­ тимальные (по критерию приведенных затрат) величины внутрен­ него диаметра D и толщины теплоизоляции трубы водопровода для подачи горячей воды.

Целевая функция имеет вид

F = — + 6ln[(D + h)jD]-1+ cD + dli

Z55

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

Параметры задачи: а=10,0; 6=1,0; с=0,5; d=l,0.

39

Таблица II

40

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