![](/user_photo/_userpic.png)
книги из ГПНТБ / Баженов, Ю. М. Перспективы применения математических методов в технологии сборного железобетона
.pdfСтрока 1 в табл. V.2 является направляющем; из плана исклю чается вектор лг5, имеющий в этой строке коэффициент 1. Выделим направляющую строку в табл. V.2 жирной линией.
На пересечении строки / и направляющего столбца находится ключевой элемент аи, в нашем примере его значение равно 25.
3. Проведем пересчет симплексной таблицы таким образом, что бы на месте ключевого элемента получить 1, а все остальные эле менты направляющего столбца обратить в нули.
Преобразование элементов симплексной таблицы производится по двум формулам: одна — для пересчета элементов направляющей строки, другая — для пересчета всех остальных элементов таблицы.
Новое значение элементов направляющей строки вычисляется по формуле
ИОВ ___ |
XI і |
(Ѵ.12) |
|
хи — |
аік |
||
|
Это правило применимо и для элемента столбца .ѵ0 н направ ляющей строке.
Для преобразования всех остальных элементов матрицы коэф фициентов, столбца д'о н целевой строки применимо одно общее правило:
ХІІ
'/ ч Xlk
x . k для і Ф І . |
(Ѵ.13) |
Введем вспомогательный коэффициент ß = — , который остает-
Xlk
ся постоянным при вычислении новых значений коэффициентов каж дой строки. Тогда
(ѵ-14)
4. Повторяем операции 1—3 до тех пор, пока в целевой строк не останется отрицательных чисел (в основных векторах-столбцах), что будет означать получение оптимального решения.
Повторяемость операций при пересчете симплексной таблицы позволяет эффективно использовать ЭВМ.
В настоящем параграфе для понимания метода приведен пример ручного счета. Для расчета более сложных таблиц можно исполь зовать стандартные программы для ЭВМ.
Рассмотрим порядок расчетов па примере пересчета табл. Ѵ.2. Вычисляем значение коэффициента ß:
160 |
= 6,4; |
250 |
= 10; |
28 |
= 1,12; |
—1 |
— |
— |
— |
— = — 0,04, |
|||
25 |
|
25 |
|
25 |
|
25 |
Запишем значения коэффициента ß в крайний правый столбец |
||||||
симплексной |
таблицы. |
Затем |
рассчитаем |
новую таблицу (первая |
||
итерация), |
причем |
первую строку — коэффициенты эффективности |
||||
Cj — можно |
опустить. Первые два |
столбца |
новой таблицы заполня |
ются аналогично исходной, кроме х5, исключенного из плана. На его место записывается вектор Х\ и соответствующее этой переменной значение с, = 1 вписывается в столбец cj.
Остальные элементы таблицы пересчитываются.
170
Пересчет начинаем |
с |
направляющей строки — все элементы ее |
||
ключевоп элемент: |
|
|
||
1100 |
25 |
— I • |
10 |
9 |
----- _ 44 1 |
■— |
— |
— о ,4; — |
|
25 |
25 |
|
25 |
25 |
|
15 |
|
_ |
1_ |
----- = 0 ,6 ; |
25 |
0,04. |
|
|
25 |
|
~ |
|
|
Элементы второй строки определяются по формуле |
(Ѵ.14): |
|||
20 000— 1100.6,4 |
= 12960; |
160— 25-6,4 = |
0; |
|
150 — 10.6,4 = |
86; |
1 6 0 - 9 -6 ,4 = 102,4; |
|
|
1 9 0 - 15-6,4 = |
94; |
0 — 1-6,4 = — 6,4. |
|
Подобным образом пересчитываются элементы третьей, четвер той II целевой строк.
Полученные значения вписываются в соответствующие клетки симплексной таблицы первого приближения (табл. Ѵ.З).
Значение функции цели задачи получаем в целевой строке в столбце Хо, оно равно 44. Однако план может быть улучшен, так
как в целевой строке имеются отрицательные элементы. Проводим вторую итерацию. В план вводим вектор х3, которому соответствует наибольший отрицательный элемент в целевой строке.
Вновь определяем вектор, подлежащий исключению, для чего вы
числяем элементы столбца а: |
|
|
|
|
|||
44 |
122- |
12 960 |
10 000 |
= 47,6; |
970 |
98. |
|
----------- = |
---------------- |
126,6; |
= |
|
|||
0,36 |
’ |
102,4 |
210 |
|
9,92 |
|
|
Минимальным является а в строке 3, |
следовательно, |
вектор |
х 2 |
||||
исключается |
из |
плана. |
Таким образом, |
направляющая |
строка |
3, |
ключевой элемент 210. Дальнейшие расчеты и составление новой симплексной таблицы проводятся по рассмотренному выше алго ритму.
Результаты расчетов показаны в таблице второго приближения
(табл. Ѵ.4).
Затем проводим третью итерацию (табл. Ѵ.5): вводим вектор х2,
исключаем вектор Л'8.
Среди элементов целевой строки в табл. Ѵ.5 есть только один отрицательный, но он находится в столбце дополнительного векто ра *5, ранее исключенного из плана. Поэтому его введение в план
нецелесообразно. Следовательно, план является оптимальным:
X — (л^; х2; .Ѵд) = (13; 57; 23,1).
Оптимальным планом в заданных конкретных условиях являет
ся выпуск конструкций: |
тип а....................... 13 |
шт. |
|
I |
|||
II |
» |
...........................57 |
» |
III |
» |
...........................23 |
» |
Найденные значения переменных -Ѵі= |
13 |
шт., х2 — Ы- шт., |
х3 = |
= 23 шт. подставляем в уравнения ограничений и проверяем расход
171
172
Т а б л и ц а Ѵ.З. Первая итерация
Т а б л и ц а Ѵ.5. Третья итерация
ч
ч
ч
ав
ч
ч
н
©
оГ
с?
О
m Cl
о
о*
1 —0,0002
О
О
о
о
Ч
-
—0,233 —5,0!
сч
о
1
о
о
о
о
о
ю
ч
о
г?
ю
1
0,0072
• о
ю
ю
о
о*
1
■5-
т
о
о
<М
ч
-
0,12 |
со |
|
+ |
||
|
||
сГ |
0,0012 |
|
to |
|
|
8 |
|
|
1 |
I 1 |
|
|
оо
ок |
со |
о |
|
о |
о |
1 |
1 |
СО |
|
1’ |
o ’ |
|
оо
о
оо
‘О
О)
£'
-
173
Т а б л и ц а V.6. Исходная симплексная таблица
|
|
|
Марка бетона |
|
|
Наибольшее |
|||
Вид |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
количество |
|
материала |
1 |
|
2 |
|
/ |
|
|
n |
песка данного |
|
|
|
|
|
типа |
||||
1 |
“и ьп |
си |
°І2 Ь12 |
а 1/ |
Ь1І |
СІІ |
aln |
bln |
cin |
Х1І |
*Ѵ1Э |
|
x lf |
|
|
x ln |
|
||
|
|
|
|
|
|||||
2 |
°21 Ь21 С21 |
а 22 Ь22 С22 |
a2j |
К / |
сг/ |
°2n |
П C2П |
||
|
+ |
1 |
Л*22 |
|
А2/ |
|
|
X2„ |
|
і |
|
|
|
аІІ |
bii |
cii |
|
|
+ |
|
|
|
|
х іІ |
|
|
|
||
|
|
|
|
|
|
|
|
|
|
ш |
|
|
|
|
|
|
amn bmn cmn |
||
|
|
|
|
|
|
|
x mn |
А ,п |
|
|
|
|
|
|
|
|
|
|
|
Потреб |
|
|
|
|
|
|
|
|
|
ность |
|
|
|
|
|
|
|
|
|
в бетоне |
|
|
|
|
P i |
|
|
P n |
|
в ж3 |
|
|
|
|
|
|
|
|
|
ресурсов при изготовлении конструкции |
по |
рассчитанному плану: |
|||
0 ,25-13 + |
0,1-57 + 0 ,0 9 - |
2 3 = |
11 |
чел.-дн.; |
|
1,6-13 |
+ |
1, 5 - 57+ 1, 6 |
- 23= |
144,5 м3; |
|
2,5-13 |
+ |
1,9-57 + 3 -23 = 210 |
чел.-дн.; |
||
0,28-13 + 0,24-57 + 0,2-23 = |
21,9 т. |
Из уравнений видно, что план оптимален, так как полностью исчерпаны возможности по изготовлению арматурных каркасов и из делий и максимально возможности автотранспорта. Производство бетона имеет резерв, и бетонная смесь может отпускаться непосред ственно па строительство.
Пример V.2. Рассмотрим еще одну возможность решения зада чи по планированию производства бетона методом линейного про граммирования [67].
Имеется завод, где изготовляются железобетонные конструкции. Для них требуется бетон различных марок. Песок для бетонов по ступает из различных карьеров, щебень для всех бетонных смесей используется один и тот же.
Применяя песок из разных карьеров, для каждой марки бето на можно подобрать смеси нескольких составов — по числу карьеров песка. Каждому составу будет соответствовать свой расход материа лов и компонентов смеси. Имея заданные потребности в бетонной смеси каждой марки и ограниченный объем отпускаемого песка од ного из карьеров, необходимо определить объем применяемой бетон ной смеси каждого вида при условиях:
174
tr |
и |
|
|
“ |
Ч) |
|
|
Ч с |
|
|
|
со |
о |
|
|
|
|
|
|
|
О? |
|
|
|
м |
со |
со |
|
І4» |
|
|
Ѵ.2 |
О |
TJ* |
о |
|
О |
||
|
|
|
|
примерукДанные |
|
|
со |
|
|
|
оГ |
Ѵ.7. |
|
|
|
Т а б л и ц а |
|
|
|
\Sо оt-
йСиЮ> «и
расход песка на I ж3 бетонной смеси в г; b — расход цемента на 1 ж3 бетонной смеси
а \о
\г
а
о
к
о
ь
о
\о
О
1/5
а) суммарный расход цемента па весь объем бетона наименьший (первый вариант);
б) стоимость всех материалов, израсходованных на приготовле
ние смеси, |
минимальная (второй вариант). |
||||
Эту задачу можно представить в виде табл. V.6. |
|||||
Обозначения: |
|
бетона |
марки /, приготовляемого па |
||
|
Xij |
— количество |
|||
1 , 2 , . . , / , . . . , а |
песке |
/; |
|
|
|
— марки |
бетона; |
|
|||
A lt А.,.......... |
|
|
|
|
|
Аі, . . . . |
Агл— наибольшее количество песка вида /, отпускаемое |
||||
Ръ Р*....... |
для бетона; |
причем А т—>-оо (т. е. не ограничено); |
|||
|
|
|
|
||
Pj, . . |
Рп — потребность в бетоне марки / в м3; |
||||
|
с і ; і |
— расход |
песка вида |
і на 1 м 3 бетонной смеси ви |
|
|
|
да /; |
|
|
1 лі3 бетона марки / при при |
|
Ьц — расход цемента па |
||||
|
|
готовлении ее па песке вида /; |
|||
|
c ; j — суммарная |
стоимость материалов, составляющих |
|||
|
|
1 м3 бетона марки / |
при приготовлении ее па песке |
||
|
|
вида I. |
|
|
|
Целевая функция будет иметь следующий вид: |
|||||
1. |
При минимальном |
расходе цемента и ограничении песка вид |
|||
і (первый вариант) |
|
|
|
пт
£Е Ьц' х-J = мни.
£=1 /=1
2.При наименьшей суммарной стоимости бетонных смесей всех видов и ограничений песка вида
|
т |
п |
|
|
|
|
|
Е |
Е Cjjxti = mm. |
|
|
|
|
|
■=1 /=і |
|
|
|
|
|
|
При этом должны быть выполнены следующие условия: |
|
||||
|
/ 2 |
|
|
|
|
|
|
Е а,-1 лу/ |
Л,- |
(/ = 1 , 2 , . . . . , |
///); |
|
|
|
/=1 |
|
|
|
|
|
|
:п |
|
|
|
|
|
|
S а,-/ — Рі |
( / = 1 , 2 , . . . /!); |
|
|
||
|
/—I |
|
|
|
|
|
Xij |
в обоих случаях ^ 0 . |
|
|
|
|
|
|
Для конкретного примера примем следующие данные: |
|
||||
|
песка /1| имеется не более |
500 и/3, у, = 2,6 |
т/м3, -уі = 1,56 |
т/м3. |
||
Сі=4,8 руб./т; |
|
|
|
|
|
|
|
количество песка /12 не ограничено, у2= 2,6 |
т/м3, |
у 9= 1,153 |
т/м3, |
||
с2 = |
2,9 руб/т; |
|
|
|
|
|
|
щебень гранитный у щ= 2,6 |
т/м3, у1Ц= 1,48 т/м3, |
Сщ= ] 1 руб/т; |
|||
|
цемент марки 500, сц= 1 2 |
руб/т; |
|
|
|
176
бетон требуется |
готовить: марки 100 — в объеме |
1000 м3, мар |
|
ки 300— 1000 л 3 и марки 500— 500 лі3. |
|
||
Бетонная |
смесь |
жесткая. |
|
Проведя |
расчет |
по составу бетона, получаем все |
необходимые |
данные, которые можно записать в виде табл. Ѵ.7. Анализ данных показывает, что эта задача может быть решена симплексным мето дом с искусственным базисом.
Для составления симплексном таблицы напишем целевую функ цию и условия ограничения для конкретных данных.
Первый вариант — при минимальном расходе цемента: 0 ,195лц + 0,254*12 + 0,378*із + 0,2І4*21 +
+ 0 ,180*22 + 0,4 1 3*2з — L = мин.
Условня orраничеиня:
0,63*іі + |
0,58*12 + 0,478*із |
|
+ |
780; |
|
0,58*2і + |
0,525*22 + |
0,416*23^5000; |
|
*и |
+ *2і |
|
=2000; |
|
|
* 1 2 |
+ * 2з |
= 1 0 0 0 ; |
|
|
* і з |
|
+ * 23= |
5 0 0 . |
Второй вариант— при минимальной стоимости бетона:
19,9*іі + 20,46*1а + |
21,59*із + 19,07*2і + |
+ 19,87*22 + 21,26*23 = L = мин. |
|
Условия ограничения те же. |
Для обоих вариантов *ü > 0 . |
Из уравнений ограничений |
и коэффициентов целевой функции |
составляем симплексную таблицу и методом последовательных ите
раций отыскиваем оптимальное |
решение: |
|
|
|
Для первого варианта |
|
|
*ті=3’20 м1 |
*,„=1000 мл |
780 т |
|
+ = 2 0 0 т |
580 т |
Л'із—0 |
|
л-2,=1680 лП |
л'„=0 |
.v„3=500 лО |
1187 т |
/1.=980 т |
— |
207 г |
|
2000 |
1000 |
500 |
- |
Значение целевой функции равно 839 г.
Анализ решения показывает, что при любых других соотношени ях расход цемента увеличивается. Например, если готовить бетоны только на песке Л2, то расход цемента составит 914 т. Таким обра зом, экономия составляет 75 т, т. е. 8,2%.
При других вариантах, кажущихся оптимальными, также полу
чается перерасход цемента в пределах 3—4%. |
|
12— 1023 |
177 |
Д л я в то р о го в ар и ан та
даі=0 |
ду., = |
о |
АУ,=0 |
- |
дг., = 2000 др |
а-., = |
1000 др |
д\3 =.= 5000 лГ1 |
1893 г |
Л. = 1 1 6 0 т |
|
52о т |
208 г |
|
|
|
|||
2000 |
|
1000 |
500 |
- |
Значение целевой функции равно 68 640 руб. |
составляет |
|||
Экономия |
по сравнению с |
другими вариантами |
до 2,5%.
Ѵ.З. Графический метод решения задач линейного программирования. Если в задаче имеются всего два неизвестных, то ее решение может быть найдено графи чески.
Рассмотрим в качестве примера задачу производства при ограниченных ресурсах.
Пример Ѵ.З. Завод выпускает изделия двух видов. Для их вы пуска требуются определенное количество материалов (цемента и ар матуры) и определенные затраты труда и энергии. Ресурсы мате риалов, рабочей силы и энергии ограничены. Необходимо устано вить, какое количество изделий каждого вида надо изгото вить, чтобы стоимость продукции была наибольшей при данных про изводственных возможностях.
Исходные данные сведены в табл. Ѵ.8.
Т а б л и ц а Ѵ.8. Данные к примеру Ѵ.З
|
|
Показатель |
Тип изделия |
Имею- |
||
|
|
I |
И |
іднеся |
||
|
|
|
|
ресурсы |
||
Стоимость одного, изделия в руб. |
250 |
400 |
_ |
|||
Расход |
цемента в |
кг . . . . |
80 |
170 |
20 000 |
|
Расход |
арматуры |
в кг . . . |
67 |
55 |
8 000 |
|
Трудоемкость |
изготовления |
30 |
40 |
6 000 |
||
в |
ч е л .- ч ................................... |
|||||
Расход |
энергии в квт-ч . . . |
52 |
67 |
12 000 |
||
Введем |
следующие обозначения: |
|
|
|
||
ду — число изготовляемых изделий /'-того типа; |
|
|
||||
с/ — стоимость одного изделия /-того типа; |
і-того вида; |
|||||
ö |
— затраты на одно изделие /-того типа ресурса |
Ві — общий объем имеющихся ресурсов t-того вида.
178
В общем виде уравнение целевой функции
п
£ |
с . х . = В-»макс. |
(V. 15) |
/= 1 |
' |
|
Система ограничений, отражающая требования, чтобы возмож ности ми по одному виду ресурсов не были превышены, имеет вид
п
£ b. X . <й,(£ = 1 , 2 , 3 , 4). |
(V. 16) |
/= 1
Таких неравенств должно быть столько, сколько имеется видов ресурсов.
Неизвестные не могут также быть отрицательными величинами. Поэтому добавляется еще система ограничений
* , ><) ( / = 1,2). |
(V.17) |
Рис. Ѵ.1. Схема решения задачи линейного программирования гра фическим методом
В конкретных условиях задачи имеем следующую систему урав нений и неравенств:
Целевая функция 250^ + 400„ѵ2 = Z - >макс. |
(V. 18) |
|||
Ограничения: по расходу цемента 80хг + |
1 70а'2 < |
20 000; |
(Ѵ.19) |
|
» |
» арматуры 67^х + |
55лг2 < |
8 000; |
(Ѵ.20) |
12* |
|
|
|
179 |