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

книги из ГПНТБ / Отраслевые автоматизированные системы управления

..pdf
Скачиваний:
12
Добавлен:
24.10.2023
Размер:
11.07 Mб
Скачать

Пусть имеется единая классификация выпуска и ввоза. Тогда, для того чтобы построить матрицу затрат-выпуска систем N предприятий, необходимо установить валовой выпуск отрасли, промежуточные поставки и конечное сальдо выпуска-затрат. Если ?/? — товарный выпуск г-го продукта на к-м предприятии, a s? — потребление, то ва­ ловой выпуск і-го продукта по всем N предприятиям бу-

N

дет следующий: yi = 2 У* (суммируются выходные величины второго квадранта). Потребление того же продук­

та исчисляется

как

сумма входных

величин

третьего

 

N

 

 

 

 

 

 

квадранта: s4 =

2 4 -

Общий баланс

і-го продукта

в си-

стеме определяется как разность zt =

yt

Если

zt^>0,

то продукт выпускается предприятием; если zt

О, то

продукт промежуточный; если zt <^ 0,

то

продукт

вво­

зится.

 

 

 

 

 

 

 

Полная величина затрат j-ro продукта на

производство

/-го определяется, очевидно, из равенства:

 

 

 

JV

N

 

 

 

 

 

 

Уи=1>Уц=2ьиу*,

 

 

 

 

 

 

(4.45)

K-=l

fc=l

 

 

 

 

 

 

где Ьи — коэффициенты полных затрат і-й продукции на производство единицы 7-й продукции. В результате имеем

балансовое соотношение по

отрасли:

yi = 2i/y + z i - .

( 4 - 4 6 )

і

 

Величина zt в данном случае равна определенной выше ве­ личине сальдо:

ИіьШ: = 2 s i : = si-

Однако

Уі = 2 Уи,

Ki

і

поэтому

2 УІЗ s i -

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

Ввоп

 

 

ач =

хи

 

 

 

 

гк

 

*+ і і

1

 

 

 

\k = N

А, =о ' Нет

 

 

 

 

 

 

Да'

 

 

 

 

Обращение

 

1 вкалывание

!>-*,*]

 

 

 

 

 

 

Л * = 0

' Н ет

 

 

 

д«

[ - 4

 

 

 

 

 

 

 

Е-Ак\

 

 

 

 

->

 

Рис. 13

 

 

 

 

Алгоритм

включает

два обращения

матриц типа

(Е — А), поэтому целесообразно исходные

матрицы вы­

разить в стоимостной форме. В этом случае при использо­

вании

итеративного алгоритма

Вп

= Е + АВп_г; В0 = Е

может быть гарантирована высокая точность обратной матрицы (Е — А)'1 (через В обозначена матрица

(Е - А)-\

5* 131

Пусть построены

матрица |] уц

 

zt и

определена

матрица

прямых затрат

где аи

= уи1У]

Первый

квадрант

матрицы прямых

затрат подразделя­

ется на 9 блоков. Диагональные блоки относятся к основ ному, вспомогательному и конкурирующему производст вам. Первый квадрант матрицы выглядит следующим об разом:

Аи

Ліз

Аіз

 

-•12,

и 2 3

 

-•Ьз

Ass

А,

Ач

Аз

Вспомогательное производство легко исключить, вос­ пользовавшись уравнением баланса:

А2іУі + А^Уч + А*-зУз = Уі-

Откуда:

у2 = [Е — А22]'1А21у1

+ [Е —

А^^А^уз.

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

{Аи

+ Ла

[Е -

Л ^ ] - 1 А21}

ух

+

13

+

А12

-

 

 

— А^]-1

А23} ys + z =

y1;

 

 

 

 

 

 

{A3l

+ Азг

[Е -

Ам]-1

А21]

У і

+

33

+ AS2

-

 

 

Аи]

~1

Аз} Уз =

Уз +

z;

 

 

 

 

 

 

[А, + А2 -

4 м ] - 1 А21х

 

+ {As +

А2

[Е -

А22]^

X

 

xA23}y3

 

=

s.

 

 

 

 

 

 

 

 

 

Обозначая

 

величины

в

аъщрешы^.

скобках

через

^л» ^ізі ^ з и

-^зз! Вг

и В3,

 

получаем

преображенную

схему первого квадранта матричной модели:

Bv.

Віз

Ви

В-л

Вх

В3

Конкурирующее производство вводится в модель при условии, что продукт уя производится в количестве, недостаточном для потребности системы, и разницу при­ ходится возмещать за счет поставок со стороны в коли­ честве, равном z. Общее потребление конкурирующего продукта в системе определяется из уравнения баланса

В31Уі = IE — В33] у3 + z.

Итак, величина конкурирующего продукта исчисля­ ется по формуле:

У3 + 2== [Е — В33] 7/з +Z s".

Конкурирующий продукт можно рассматривать как ввоз со стороны. Складывая конечный продукт системы с по­

ставками в конкурирующее производство Z =

z + В13

+

у3

и вычитая из общего ввоза поставки в конкурирую­

щее

производство s' = s — В3у3, мы можем

далее

про­

анализировать систему без конкурирующего производства. Конкурирующее производство в свою очередь можно рассматривать как самостоятельную систему.

Таким образом, схема матричного моделирования обес­ печивает подготовку и организацию в экономичной форме потока данных для комплекса экономико-математических моделей ОАСУ.

6.Моделирование процессов оптимизации производственной структуры

Первые этапы разработки и внедрения ОАСУ пред­ полагают механизацию потоков информации и оптимиза­ цию проектно-плановых решений в основном на базе традиционно сложившихся производственных структур и организации управления. Между тем ОАСУ — это

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

Исследования, направленные на выработку новых про­ изводственно-организационных форм, представляются весь­ ма актуальными и привлекают внимание различных коллективов специалистов. В исследованиях путей совер­ шенствования производственных и организационных струк­ тур, как и оптимального планирования и управления, большое место отводится методам и средствам экономикоматематического моделирования. Однако математическая постановка задач улучшения структуры производства и уп­ равления, как правило, значительно сложнее постановки задач оптимального планирования при задаиных произ­ водственных структурах. Задачи оптимизации структу­ ры — это задачи синтеза, тогда как задачи оптимального планирования при заданной структуре являются задачами анализа систем управления. Из общей теории управления известно, что задачи анализа, как правило, проще задач синтеза.

В настоящее время разрабатываются комплексы эко­ номико-математических моделей функциональных подсис­ тем ОАСУ, в то же время совершенствование производст­ венных структур с помощью экономико-математических методов ограничивается в основном постановкой отдель­ ных задач математического программирования. Несмотря на известную ограниченность этих задач, в большинстве случаев в них анализируются наиболее актуальные вопро­ сы поставленной проблемы, поиск вариантов совершенст­ вования использования оборудования, выбор оптималь­ ных способов размещения, строительства и специализа­ ции производственных комплексов и т. п.

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

квыработке программы перспективной реорганизации

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

В качестве примера рассмотрим задачу оптимизации загрузки оборудования, постановка которой определяется следующим. По ряду групп станков народное хозяйство испытывает дефицит. Между тем учет использования

станков на предприятиях свидетельствует о значительной неравномерности их загрузки 8 . Для различных видов станков загрузка колеблется от 18—25 до 80—90% общего фонда рабочего времени. Одна из главных причин такой неравномерности заключается в традиционно сложившей­ ся жесткой технологической привязке деталей к станкам и практической реализации в связи с этим неоптимальных

технологических

маршрутов.

 

 

Цель задачи

оптимизации

загрузки

оборудования —

— установление для заданного

парка оборудования и произ­

водственной

программы оптимальных

технологических

маршрутов

(«деталь — станок»), при реализации которых

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

гических

маршрутов производства.

 

 

Итак, рассматривается

производственное подразделе­

ние

(цех, участок,

предприятие), оборудованное

т стан­

ками). На

этих

станках

в

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

времени

Т

необходимо обработать

п

деталей. Технологические

маршруты прохождения деталей через станки необходимо

^.^прстроить так,

чтобы

минимизировать

сумму

штрафов

за

непроизводительное

применение станков

(частично

или

полностью

загруженных), которые

участвуют в про­

изводственном процессе. Под непроизводительным ис­ пользованием станков понимаются их простои и перена­ ладки.

 

Введем

следующие

обозначения: і — номер

станка;

/ — номер детали;

к — номер производственной операции

по обработке /-Й детали

= 1, 2, . . ., pj); а і 1

К

— время

выполнения

/с-й операции по

обработке /-й

детали на

г"-м

станке;

т;^- — время

переналадки г-го станка

после

выполнения

на нем Z-й

операции обработки

s-й

детали

(1

p s

) для

выполнения

к-ж операции

обработки

"-Й

детали

(1 < : к

<: p j ) .

 

 

 

 

8 Н. В. М а х р о в, А. А. М о д и н, 10. А. О л е й н и к, В. В. Ю ш к о в . Вопросы организации производственных комп­ лексов. М., изд. ЦЭМИ АН СССР, 1969 (ротапрпит).

Неизвестные величины, подлежащие определению, сле­

дующие:

tjh — момент

начала

обработки j-jfc детали по

к-к операции (tjh >

0);

гцк

— булевская

переменная,

равпая 1,

если к-я операция для /-й детали осуществляется

на г-м станке, и равная

0 — в противном

случае.

Еслп

технологические

маршруты («деталь — станок»)

установлены заранее, то задача сводится к построению календарного графика загрузки станков определенными деталями. Искомые величины г,- j i t теряют при этом смысл; величина tjh обозначает по существу момент начала - обработки /-и детали на некотором i0(k) станке, так как заранее известно, что к-я операция для 7-й детали должна выполняться на г"0 станке. В более общей постановке задачи основные ограничения представляются в виде следующих условий.

Возможность выполнения к-я операции для обработки /-й детали только на одном станке определяется следую- •

щим

образом:

 

7 = ] , . . . , п;

 

 

Р,

7 >

Для каждой детали последующая операция не начи­ нается прежде, чем не окончилась предыдущая операция:

 

т

.

j

^

W >

2 «и**и*,

fc,

+ і е {і,..., РЛ. ( 4 - 4 8 )

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

zUkzisi

Iізк— Vsi + aisi

"I" WJj)l [tsi

— (fjk + «Uk +

 

 

 

 

 

 

 

 

 

(4.49)

і = 1 , . . .

,m;

j,s

= 1 , . . . ,

n;

 

 

 

= 1 , . . . ,Рз-;

I =

1 , . . .

, p s

;

кфі.

 

 

Ограничения

(4.49) достаточно

сложны,

поэтому

сде­

лаем некоторые

пояснения.

 

 

 

 

 

Если

хотя бы одна

из операций — к-я

по обработке

7-й детали

или 1-я по

обработке

s-й детали (/ п s могут

совпадать)

не выполняется

на

г-м станке,

то zijh

0

пли zisi — 0. Выражение в левой части (4.49) при этом обращается в нуль, следовательно, ограничения (4.49)

удовлетворяются.

 

 

 

 

 

 

Рассмотрим случай, при котором обе операции

выпол­

няются

на і-м станке.

 

 

 

 

 

1-й

вариант.

1-я операция по обработке

s-я

детали

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

к-я

операции для /-Й детали. Если

деталеоперация (j — /с)9

следует

непосредственно

за деталеоперацней

(s — I),

то

разница

между

началом

деталеонерации

(/ — к)

tjh

и концом деталеоперации (s— l)tsi

+ aisl составляет время

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

ду окончанием деталеоперации

(s — I) и

началом

детале­

операции (/ — к)

выполняется еще хотя

бы одна

детале­

операция

{г — g),

то

 

 

 

tjk >tsi

+ аы

+

x[fr + airg +

т $ .

 

(4.50)

Ограничимся

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

технологических

процес­

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

J_z~ к и s — Z) предполагается

следующее условие:

xfrj "~Ь ""Чет ~\~ ^ІГ£=^ T4js-

(4.51)

Таким образом, из неравенств (4.49) и (4.50) следует, что

* i * > t . i + в ы + 4 ' .

С 4 - 5 2 )

Иллюстрация данного случая на оси времени (рис. 14) показывает, что при непересечении деталеоперации (/ — к) и (s — I) ограничения (4.49) удовлетворяются, так как выражения в квадратных скобках всегда имеют противо­ положные знаки.

hi

hi + aisl +

 

lJk

lJk + aijk

+ *ijs

1

1

-

1

1

a*-

Рис. І4

11 Термин «деталеоперация» вводится для краткости.

2-й вариант. Деталеоперацпя (/ — к) начинается рань - lire, чем завершается деталеоперацпя (s — Г). На рис. 15 видно, что выражения в квадратных скобках неравенства (4.49) имеют одинаковые знаки, следовательно, в этом случае данное ограничение не выполняется.

Рис. 15

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

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

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

Ограничения по срокам выполнения операций характе­ ризуются неравенствами:

 

т

 

 

 

hk +

2 aUk zijk

^

Т,

(4.53)

 

i=l

 

 

 

І = 1,

. . ., щ

к

= 1, . . .,

pj.

Ограничения по загрузке станков и длительности технологического процесса обработки деталей таковы:

 

Pj

т

 

п

Pj

z'1^ ^Г '

(4'54)

max (max ^

^ а2 Шіm a x

2

2

]

к=1

i = i

1

j=i

K-=i

 

 

Pj

т

 

 

 

 

 

 

г Д е max 2

2 aUkzijk

— наибольшая длительность

Texa­

s' к =1

І = І

 

 

 

 

 

 

1 0 «Календарное планирование». М., изд-во «Прогресс», 1966, гл. 12.

138

 

n

Pj

логии полной

обработки детали; max 2

S аШгИ1

максимальная

г j = i

k = i

загрузка станка.

 

Легко заметить, что ограничения (4.54) следуют из ограничений (4.53), поэтому они не являются независимы­ ми. Так, ограничения (4.53) показывают, что обработка каждой детали (по любой из необходимых для этой детали операций) должна закончиться не позднее Т. Следователь­

но, к моменту Т все станки должны быть

свободны.

А раз так, то и максимальная производственная

загрузка

станков, и наибольшая длительность технологии

обработ­

ки деталей не

должны превышать Т.

Это и

отражено

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

(4.54). Тем не менее мы

включаем ограни­

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

В задачу оптимальной загрузки оборудования можно ввести более специальные ограничения, такие, как необ­ ходимость завершения обработки некоторых деталей к более раннему сроку Т' (Т' <| Т); частично упорядочен­ ные технологические маршруты обработки деталей (вместо

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

к = 1, 2,

. . ., р, могут

быть вариан­

ты с перестановкой

этих номеров) и др.

 

 

 

Сформулируем целевую функцию задачи оптимальной

•—загрузки оборудования. Введем показатели

С; — штраф

за

единицу времени непроизводительного использования

і-го

станка.

 

 

 

 

 

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

 

 

 

/—>гаіп,

 

 

(4.55)

 

где

 

 

 

 

 

 

 

 

(4.56)

где х — аргумент функции Д

в выражении

(4.56), et

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

і-го

станка, при которой штраф за непроизводительное исполь­ зование этого станка не включается в целевую функцию.

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