книги из ГПНТБ / Ильинский, Д. Я. Обоснование решений при проектировании и эксплуатации машин и линий легкой промышленности учебное пособие
.pdfРис. 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х 2— F, |
(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