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

книги / Решение некоторых многоэкстремальных задач методом сужающихся окрестностей

..pdf
Скачиваний:
3
Добавлен:
12.11.2023
Размер:
12.71 Mб
Скачать

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

Вэтих условиях единственным способом получения обоснован­ ного решения является использование математических методов и ЭВМ.

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

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

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

1)проектируемое предприятие размещается на однородной (или предварительно выравниваемой) территории;

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

Если перечисленные условия выполнены, то задача разработки

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

151

В задаче разработки архитектурно-планировочного решения про­ изводственного здания объектами являются отдельные поме­ щения.

В первом случае под областью подразумевается площадка пред­ приятия, во втором — абрис здания'в плане.

Размещаемые объекты могут быть связаны друг с другом комму­ никациями, которые на первом этапе фактически не проектируются, а только оцениваются по стои­

мости.

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

тать, что стоимость всей системы коммуникаций равна сумме стои­ мостей отдельных коммуникаций.

Оценкой варианта схемы генерального плана можно считать сум­ му затрат (единовременных, эксплуатационных или приведенных) на территорию и коммуникации. Возможны и другие варианты це­ левой функции в виде взвешенной суммы различных затрат.

В общем случае задачи построения схемы генплана и компоновки одноэтажного производственного здания являются многокритери­ альными [39].

Формальная постановка задачи заключается в следующем. За­ данный набор объектов Tt (i — 1, 2, ..., п) необходимо разместить в области £2 так, чтобы некоторая функция цели принимала экстре­ мальное значение. Размещаемые объекты Tt (i = 1, 2, ...,я) и об­ ласть Q представляют собой прямоугольники или фигуры, состав­ ленные из прямоугольников со взаимно параллельными сторонами. Кроме того, прямоугольники, образующие объекты Г, (i = 1, 2, ...

..., л), и прямоугольники, из которых состоит область Q, также имеют взаимно параллельные стороны (рис. 36). Таким образом, на область Q и объекты Т{ (» = 1, 2, ..., п), кроме требования их аппроксима­ ции прямоугольниками, иных ограничений не накладывается, т. е. они могут быть невыпуклыми, многосвязными, несвязными. -

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

Отметим, что рассматриваемая постановка задачи проектирова­

152

ния в некотором смысле совпадает с постановкой задачи, приведенной в работе [36].

Область й свяжем с неподвижной системой координат хОу, а раз­ мещаемые объекты Т{ — G собственными подвижными системами

координат 6* 01 т/ (< = 1, 2, .... п). Пусть объекты 7\ имеют пара­ метры размещения хь yt, 04 [54]. Положим, что параметры (i = 1,

2 ,..., п) могут принимать значения, кратные -у-, т. е. в процессе раз­

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

На языке математического программирования рассматриваемая задача в общем виде формулируется следующим образом.

Найти min я (z) — к (г*), где G определяется системой ограниче-

г£0

ний

 

 

g ,(z )> 0

(t = 1, 2, . . . ,

m),

г = z (х, у, 0) = z (xlt лг2, . . .

, хп, ylt yz, . . .

, уп, 0Х, 0 2 , . . . . 0„).

Как указывалось ранее, задача проектирования генеральных планов в общем случае является многокритериальной. Предполага­ ется, что все критерии я( (г) (i — 1 ,2 ,..., п) можно объединить в еди­ ную функцию цели и следующим образом:

 

 

 

* (z) =

У etXi (г),

 

 

(4.10)

 

 

 

 

i=i

 

 

 

 

т. е. функция цели представляется в вцде линейной

комбинации час­

тных критериев. При этом считается,

что коэффициенты с( (i = К

2, ..., п) можно определить аналитически, методом

экспертных оце­

нок и т. д.

 

 

которые

формируют функцию

цели

- Введем частные критерии,

х (г)

рассматриваемой задачи.

 

 

 

 

 

1. Площадь области й, которую можно представить в виде

 

 

(*> У. 9) = шах (хс+ а, (0)) max (у{ +

р{ (0f)),

(4.11)

где

 

 

»€[1.я]

 

*€[1,«]

 

 

 

а( =

max [($ -f а {) cos 0, +

(t]f -f b\) sin 0,];

 

 

 

 

Pf =

.

,

,

.

 

(4.12)

 

max [(— g/ +

a!) sin 0i + (rjf +

b[) cos 0<];

 

a'i,

bf — размеры /-го прямоугольника i-го размещаемого объек­

та; gf, iif — координаты полюса

(левого нижнего угла) /-го размеща­

емого объекта;

kt — количество прямоугольников, образующих t-й

размещаемый

объект.

 

 

 

 

 

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

153

манхэттенской метрике [84]

 

х2(*+, у+, 6) = £ '

£

* i , ( \ x ? - x t \ + \ y ? - y t \ ) '

(4.13)

 

i = i / = / + 1

 

 

 

 

 

 

 

где а>ц— стоимость единицы длины коммуникации;

 

 

 

x t =

х( +

§£+ cos 0£ +

г|/*" sin 0£;

 

 

 

y t =

уi — & sin 0£+

v f cos 0t;

 

 

If", Tii1"— координаты точки входа

коммуникации t-ro

размещаемого

объекта (см. рис. 36).

 

 

 

 

 

 

 

 

 

 

Система ограничений, образующая множество возможных зна­

чений параметров размещения, состоит из:

объектов [54]

условий взаимного

непересечения

размещаемых

 

f u

( x i <

У и

е„ X / ,

У

Ь

0/) > 0

 

(4.14)

 

(«',/=

1,2, . . .

,

п;

/ > / ) ;

 

 

условий непересечения размещаемых

объектов с границей

обла­

сти й [54]

 

 

 

 

 

 

 

 

 

 

 

ft (*t. Vi, 9f) >

0

(t =

1,2 .......... n)\

 

(4.15)

архитектурного ограничения [31], которое для поставленной за­

дачи представляется в виде равенства

 

 

 

 

 

 

шах (х{ + а £)/шах

 

 

+

&) = L,

 

(4.16)

 

>Ф,п]

 

 

 

 

 

(4.12); L — заданное

 

где «£,

вычисляются

по формулам

число;

иными словами, требуется, чтобы выдерживалось заданное соотно­ шение длин сторон области й;

условий выполнения минимально и максимально допустимых расстояний между парами размещаемых объектов, между объектами

и границей области Й, т. е.

 

 

dH < Ре/ < Гц ( * , / = 1 , 2 ..........т;

i > /),

(4.17)

(t = 1, 2, . . . ,

п),

(4.18)

где dij, гц, dit rt — заданные числа.

Очевидно, эти неравенства определяют область возможных поло­ жений полюса /-го относительно t-ro размещаемого объекта. Геомет­ рически эта область образует некоторое топологическое кольцо. На рис. 37 такое кольцо заштриховано. Заметим, что неравенства dy < < Ру и dj < р; были учтены выше. Для этого достаточно в неравен­ ствах (4.14) — (4.15) считать, что г,7 — заданные минимальные рас­ стояния между объектами, а знак этих неравенств изменить на про­ тивоположный. Таким образом, правые части неравенств (4.17), (4.18) порождают следующую систему неравенств:

ф// (*t. Vt* 0t, */, У/, 0/) >

0

(t, / = 1, 2, . . .

, m; t > /),

(4.19)

Фt (*t, У1, 9f) >

0

( / = 1 , 2 , . . . ,

n);

(4.20)

154

условий принадлежности границ заданной совокупности разме­ щаемых объектов Tt (1 = 1,2,...,/?) границе области Q. Эти условия можно выразить требованием одновременного равенства нулю мини­ мально и максимально допустимых расстояний между объектами Т{ (i = 1, 2, ...,/?) и границей области Q, т. е. равенством р, = 0 (i = = 1, 2, ..., р). Учитывая, что эти условия являются частным случа­

ем условий (4.20)

при r{ = d( =

0, получаем

 

Ь

(xlt Уь 0,) = 0

(г = 1,2 .......... р).

(4.21)

Таким образом, решение поставленной задачи сводится

к опреде­

лению минимума функции цели (4.10) на области G, определяемой

системы неравенств (4.14) — (4.16), (4.19) — (4.20).

 

Отметим, что

система ограничений (4.19) — (4.20) может быть

настолько сильной, что нахождение хотя бы одного размещения объек­ тов, удовлетворяющего такой системе, будет сложной задачей. Этот случай нами не рассматривается. Будем полагать, что система ограничений (4.19) — (4.20) не настолько жесткая и ее выполнение не представляет особых вычислительных трудностей.

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

На первом этапе получается возможный вариант размещения объек­ тов строительства методом последовательно-одиночного размещения на базе годографа функции плотного размещения (гл. 1, § 3). Началь­ ная точка размещения задается в виде последовательности из п чисел (перестановки), согласно которой происходит размещение объектов.

На втором этапе методом сужающихся окрестностей формируются новые последовательности из п чисел (перестановки).

Идею реализации годографа функции плотного размещения в методе последовательно-одиночного размещения для получения воз­ можных размещений рассмотрим на следующем простом примере. На рис. 38 изображены два размещенных прямоугольника Тг и Т.г. Пусть размещается третий прямоугольник Т3. Обозначим у31 и у32 — годограф функции плотного размещения третьего размещаемого

155

Н ом ер

 

 

 

 

 

 

 

 

 

 

 

о б ъ ­

l

2

3

4

5

6

7

8

&

10

и

екта

 

 

 

 

 

 

 

 

 

 

 

1

2

о:б

0.5

0.5

0.5

0.5

0.5

0.5

0.5

2

 

2

2.5

1.0

0.5

0.5

0.5

0.5

0.5

3

 

10.0

7

0.5

0.5

0.5

0.5

0.5

0.5

4

 

 

 

1

0.5

0.5

0.5

0.5

0.5

5

 

 

 

 

1

0.5

0.5

0.5

0.5

6

 

 

 

 

 

1

0.5

0.5

0.5

7

 

 

 

 

 

 

1

0.5<

0.5

8

 

 

 

 

 

 

 

1

0.5

9

 

 

 

 

 

 

 

 

1

10

 

 

 

 

 

 

 

 

 

11

 

 

 

 

 

 

 

 

 

12

 

 

 

 

 

 

 

 

 

13

 

 

 

 

 

 

 

 

 

14

 

 

 

 

 

 

 

 

 

15

 

 

2.0

 

 

 

 

 

 

16

 

 

 

 

 

 

 

 

 

17

 

 

 

 

 

 

 

 

 

18

 

 

 

 

 

 

 

 

 

19

 

 

 

 

 

 

 

 

 

20

 

 

 

 

 

1

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

21

 

 

 

 

 

 

 

 

 

22

 

 

 

 

 

 

 

 

 

23

 

 

 

 

 

 

 

 

 

24

 

 

 

 

 

 

 

 

0.50.5

0.50.5

0.50.5

0.50.5

0.50.5 |

0.5

0.5 |

о ся

о Ьп

0.5 [

0.5

0.5

0.5

10.5

1

-и.

4

> %

156

 

 

 

 

 

 

 

 

 

 

Т аблица

4ЛЗ

12

13

14

15

16

17

18

19

20

21

22

23

24*

0.5

0.5

0.5

0.5

0.5

0.5

0.5

0.5

0.5

0.5

 

 

0.5

0.5

0.5

0.5

0.5

0.5

0.5

1.0

1.0

1.0

1.0

 

 

0.5

0.5

0.5

0.5

0.5

0.5

0.5

1.0

1.0

1.0

1.0

 

 

0,5

0.5

0.5

0.5

0.5

0.5

0.5

0.5

0.5

0.5

0.5

 

 

0.5

0.5

0.5

0.5

0.5

0.5

0.5

0.5

0.5

0.5

0.5

 

 

0.5

0.5

0.5

0.5

0.5

0.5

1.0

0.5

0.5

0.5

0.5

 

 

2.0

0.5

0.5

0.5

0.5

0.5

0.5

0.5

0.5

0.5

0.5

 

 

0.5

0.5

0.5

0.5

0.5

0.5

0.5

0.5

0.5

0.5

0.5

 

 

0.5

0.5

0.5

0;5

0.5

0.5

0.5

0.5

0.5

0.5

0.5

 

 

0.5

0.5

0.5

0.5

0.5

0.5

0.5

0.5

0.5

0.5

0.5

 

 

0.5

0.5

0.5

0.5

0.5

0.5

0.5

0.5

0.5

0.5

0.5

 

 

0.5

1

0.5

0.5

0.5

0.5

1.0

0.5

0.5

0.5

0.5

 

 

1.0

 

1

0.5

0.5

0.5

0.5

0.5

0.5

0.5

0.5

 

 

0.5

 

 

К

0.5

0.5

0.5

0.5

0.5

0.5

0.5

 

 

------ 0~

 

 

 

 

0.5

 

 

 

1

0.5

0.5

0.5

0.5

0.5

0.5

 

 

0.5

 

 

 

 

1

0.5

0.5

0.5

0.5

0.5

 

 

0.5

 

 

 

 

 

1

1.0

1.0

1.0

1.0

 

 

0.5

 

 

 

 

 

 

3

0.5

0.5

0.5

 

 

0.5

 

 

 

 

 

 

 

2

0.5

0.5

 

 

0.5

 

 

 

 

 

 

 

 

2

0.5

 

 

0.5

 

 

 

 

 

 

 

 

 

1

 

 

0.5

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

1

167

объекта относительно первого и второго, размещенных с учетом мини­ мально допустимых расстояний, а Г31 и Г32 — годограф функции плотного размещения с учетом максимально допустимых расстояний. Область abed (показана двойной штриховкой) является областью воз­ можной постановки полюса третьего размещаемого объекта. При этом гарантируется непересечение третьего объекта с первым и вто­ рым, а также выполнение ограничений на минимально и максималь­ но допустимые расстояния. Точка а на рис. 38 является точкой плот­ нейшего расположения трех рассматриваемых объектов. Если в задаче

речь идет о плотной упаковке (а это, как правило, именно так), то естественно из всего множества точек области abed выбрать точку а.

Для 6-го шага последовательно-одиночного размещения рас­ сматриваются все возможные пары пересекающихся областей, огра­ ниченных годографами функций плотного размещения у& и Г** (I =

=1, 2, ..., 6 - 1 ) .

Пусть Ak — множество точек плотной упаковки 6-го размещае­

мого объекта. Тогда из всех точек Ak выбирается подмножество A'k точек, которые не являются внутренними для областей, ограниченных годографами функций плотного размещения у&- (I = 1,2,..., 6 — 1).

Из множества А* выбирается точка, в которой наилучшим образом выполняется архитектурное ограничение (условие 7). Если таких точек более двух, то из них выбирается та, в которой функция цели (4.10) принимает меньшее значение.

Для размещаемых объектов, которые аппроксимируются наборами прямоугольников, годографы функций плотного размещения у// и 1\/, как и в случае размещения прямоугольников, ограничивают области, изображенные на рис. 38. Однако при этом увеличивается информационность задачи и несколько усложняется логика решения. Процедура поиска точек возможного положения полюса очередного размещаемого объекта ничем не отличается от рассмотренной выше.

П р и м е р . Заданы объекты строительства, показанные на рис. 39. В табл. 4.13 приведены минимально (наддиагональные эле­

158

менты) и максимально допустимые (поддиагональные элементы) рас­ стояния между размещаемыми объектами и границей области раз­ мещения. В табл. 4.14 указаны стоимости единицы длины коммуни­ каций, связывающих объекты. Стоимость единицы площади равна

Примечание. В каждом прямоугольнике цифры над диагональю означают пары номеров це­ хов, а под диагональю — удельную стоимость, тыс. р./100 м.

144 тыс. р. Требуется разместить объекты строительства в прямо­ угольной области о соотношением сторон 2 : 1 так, чтобы стоимость генерального плана была наименьшей.

Результат решения задачи на ЭВМ «БЭСМ-6» методом сужа­ ющихся окрестностей показан на рис. 39. Время счета составило около 30 мин.

15»

§7. Плотная упаковка набора параллелепипедов

впараллелепипеде с областями

запрета

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

Рассмотрим следующую задачу оптимизации размещения гео­ метрических тел [9, 58].

Пусть в параллелепипеде Q о основанием А х В и о заданными в пространстве Q областями запрета Kt ( t= 1 ,2 .......т), имеющими форму параллелепипедов, требуется разместить п параллелепипедов Pt (t = 1, 2, ..., п) заданных размеров так, чтобы высота х занятой части Q была наименьшей. Параметры (размеры и координаты полю­

сов)

областей запрета Kt : at X bt

х ht, х", у и

г* (t = 1 , 2, ...

(i =

размеры размещаемых параллелепипедов

Pt : at X be X ht

1, 2, ..., n). Соответствующие

грани параллелепипедов Рс (i =

=1, 2, .... п), Kt (t — 1,2, ..., m) и Q взаимно параллельны. Параллелепипеды Р{ (i — 1, 2.......п), Kt (t = 1, 2, .... т) вло­

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

Будем считать, что полюсы параллелепипедов Рс, Kt и Q совпа­ дают с началом их собственных систем координат. Эго приводит к некоторым упрощениям при программном построении годографа век­ тор-функции плотного размещения [56] двух параллелепипедов.

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

зом.

 

 

целевой

функции х (Z)

 

 

Найти минимум

 

 

 

 

 

 

 

П

 

 

 

 

 

min х (Z) = min V<*fo + hi),

 

 

 

 

 

zee

 

/=i

G — множество

где

 

— символ

операции

^-дизъюнкции [521;

точек,

удовлетворяющих системе неравенств, определяющих:

 

а)

П

Pt с

£2 (i =

1, 2, ..., п), т. е.

 

 

условие U

 

 

 

 

 

 

r%(Qio) — x2i — y f> 0 ,

(4.22)

где

0«) = arcsin

 

f

;

 

 

 

 

 

V

x f +

Уi

 

 

160

Соседние файлы в папке книги