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

книги из ГПНТБ / Баженов, Ю. М. Перспективы применения математических методов в технологии сборного железобетона

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

Строка 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

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