УДК 658.512.012.011.56:665.6
Автоматизированное проектирование схем размещения объектов предприятий из условия минимизации коммуникационных затрат
И.М. Зуга, В.Г. Хомченко
ОАО «Омскнефтехимпроект»,
Омский государственный технический университет
Реферат
Предложены математическая модель оптимизационного синтеза и алгоритм автоматизированного проектирования схем размещения объектов предприятий при незакрепленных местах их возможного расположения из условия минимизации затрат на коммуникационные связи с учетом ограничений на минимально допустимые расстояния на просвет между этими объектами. Ил. 3. Библ. 7.
Ключевые слова: схемы размещения объектов, автоматизированное проектирование, уровни удельных затрат
При проектировании предприятий различных отраслей промышленности: машиностроительной, нефтехимической и других - важной задачей является рациональное размещение тех или иных объектов структурных подразделений на выделенной для них площади: станков и другого оборудования производственных участков, цехов и отдельных производств в рамках предприятий в целом.
Известно решение задачи о размещении объектов на предварительно выделенные места [1,2]. В данной работе рассматривается задача автоматизированного проектирования схем размещения объектов того или иного производственного подразделения при незакрепленных местах возможного расположения объектов.
На первом уровне абстракции будем представлять объекты в виде прямоугольников либо окружностей, охватывающих в плане контуры объектов.
В коммуникационные затраты включим все возможные затраты, связанные с созданием (проектированием, изготовлением, монтажом и т. п.) и эксплуатацией коммуникаций между объектами и необходимые для выполнения этими объектами их функционального назначения.
Так как задача о размещении объектов решается, как правило, на начальной стадии проектирования, когда известны лишь ориентировочные оценки затрат на реализацию коммуникаций между объектами, то целесообразно перейти на этом этапе от оценки затрат в тех или иных абсолютных единицах к некоторым условным уровням затрат, представляющим собой экспертную оценку удельных (приходящихся на единицу длины коммуникаций) затрат в принятой предварительно системе баллов. Принятая система баллов должна адекватно отражать ранжированную последовательность, например, по убыванию предполагаемых удельных затрат на коммуникационные связи. Тогда условные коммуникационные затраты между каждой парой объектов можно выразить следующей зависимостью:
, (i=1,…, n-1; j=i+1,…, n) (1)
где
- уровни удельных затрат на реализацию коммуникационных связей между i-м и j-м объектами; - длина коммуникаций между i-м и j-м объектами;X и Y - векторы координат
и центров i-го j-го объектов размерностью n 1;n - число объектов рассматриваемой производственной структуры.
Условные коммуникационные затраты рассматриваемой производственной системы с учетом (1) будут равны:
(i=1,…, n-1; j=i+1,…, n) (2)Поскольку, как отмечалось, решение о размещении объектов принимается на начальной стадии проектирования в условиях неопределенности [3], то целесообразно, не нарушая строгости постановки задачи, в качестве длины коммуникаций
принять кратчайшее расстояние между центрами объектов, а именно:, (i=1,…,n-1; j=i+1,…,n) (3)
где
, и , - координаты центров прямоугольников либо окружностей, представляющих соответственно i-й и j-й объекты в принятой системе координат.На последующих этапах проектирования можно использовать уточненные данные как по удельным затратам на коммуникационные связи, так и по способу расчета длины коммуникаций.
На взаиморасположение объектов производственных подразделений, как правило, накладывается целый ряд ограничений, характерных для конкретной отрасли промышленности. В данной работе без потери общности в качестве таких ограничений примем минимально допустимое расстояние между объектами на просвет, как наиболее характерное для большинства отраслей промышленности.
При решении задачи о размещении объектов удобно использовать в качестве дополнительных условий синтеза минимально допустимое расстояние между центрами i-го и j-го объектов, определяемое зависимостью
, (i=1,…,n; j=1,…,n) (4)
если объекты представлены в виде прямоугольников, и выражением
(i=1,…,n-1; j=i+1,…,n) (5)если - в виде окружностей (здесь:
- регламентированное минимально допустимое расстояние на просвет между i-м и j-м объектами; , , , и , - длина, ширина и радиус соответственно прямоугольников и окружностей, охватывающих в плане i-й и j-й объекты).С учетом выражений (3), (4) и (5) допустимая область возможного расположения центров объектов будет ограничиваться системой неравенств:
, (i=1,…, n-1; j=i+1,…, n) (6)
Для выполнения условия (6) в ходе автоматизированного проектирования схем размещения объектов воспользуемся так называемой функцией штрафа, а именно:
, (7)где
- частная штрафная функция, равная 0, если ;
если .
На основе зависимостей (2) и (7) запишем целевую функцию вида
Z=minW(X,Y), (8)
X,Y
Vгде
;V - область возможных (допустимых) значений векторов X и Y;
- параметр, позволяющий регулировать влияние функции штрафа (7) на свойства целевой функции (8).
Отметим, что векторы X и Y являются в данном случае множеством свободных параметров целевой функции Z.
Поиск минимума целевой функции (8) представляет собой типичную задачу нелинейного программирования, для решения которой можно воспользоваться теми или иными известными методами [4-7]. В данной работе для автоматизированного поиска оптимального решения задачи о размещении объектов использован градиентный метод [4], алгоритм которого представлен на рисунке 1.
В блоке 1 вводятся значения уровней удельных затрат; размеры прямоугольников либо окружностей, представляющих контуры объектов в плане; минимально допустимые расстояния между объектами на просвет и, при необходимости, другие данные. В блоке 2 свободным параметрам синтеза
и (i=1,…,n) присваиваются значения, соответствующие начальному приближению оптимизационной задачи. Блок 3 предназначен для расчета целевой функции при начальном приближении свободных параметров. В блоке 4 численно определяются частные производные по параметрам и (i=1,…,n), а в блоке 5 – новые значения свободных параметров в направлении антиградиента. После расчета целевой функции (8) при новых значениях свободных параметров в блоке 6 ее значение сравнивается в блоке 7 со значением, полученным на предыдущем шаге l итерации, и в зависимости от результата сравнения расчет либо продолжается с запоминанием нового значения целевой функции в блоке 8, либо завершается. В последнем случае в блоке 9 выводятся в той или иной форме результаты решения данной оптимизационной задачи.