книги из ГПНТБ / Отраслевые автоматизированные системы управления
..pdfПусть имеется единая классификация выпуска и ввоза. Тогда, для того чтобы построить матрицу затрат-выпуска систем 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 |
— |
||
показатель минимальной производственной загрузки |
і-го |
||||
станка, при которой штраф за непроизводительное исполь зование этого станка не включается в целевую функцию.
