Methodology to Solve Multi-Dimentional Sphere Packing Problems
This paper discusses the problem of optimally packing spheres of various dimensions into containers of arbitrary geometrical shapes. According to the international classification, this problem belongs to Sphere Packing Problems (SPPs). The aim of this work is to create an integrated methodology fo...
Збережено в:
| Опубліковано в: : | Проблеми машинобудування |
|---|---|
| Дата: | 2019 |
| Автор: | |
| Формат: | Стаття |
| Мова: | English |
| Опубліковано: |
Інстиут проблем машинобудування ім. А.М. Підгорного НАН України
2019
|
| Теми: | |
| Онлайн доступ: | https://nasplib.isofts.kiev.ua/handle/123456789/158842 |
| Теги: |
Додати тег
Немає тегів, Будьте першим, хто поставить тег для цього запису!
|
| Назва журналу: | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| Цитувати: | Methodology to Solve Multi-Dimentional Sphere Packing Problems / G.N. Yaskov // Проблеми машинобудування. — 2019. — Т. 22, № 1. — С. 67-75. — Бібліогр.: 37 назв. — англ, рос. |
Репозитарії
Digital Library of Periodicals of National Academy of Sciences of Ukraine| id |
nasplib_isofts_kiev_ua-123456789-158842 |
|---|---|
| record_format |
dspace |
| spelling |
Yaskov, G.N. 2019-09-14T14:30:47Z 2019-09-14T14:30:47Z 2019 Methodology to Solve Multi-Dimentional Sphere Packing Problems / G.N. Yaskov // Проблеми машинобудування. — 2019. — Т. 22, № 1. — С. 67-75. — Бібліогр.: 37 назв. — англ, рос. 0131-2928 DOI: https://doi.org/10.15407/pmach2019.01.067 https://nasplib.isofts.kiev.ua/handle/123456789/158842 519.85 This paper discusses the problem of optimally packing spheres of various dimensions into containers of arbitrary geometrical shapes. According to the international classification, this problem belongs to Sphere Packing Problems (SPPs). The aim of this work is to create an integrated methodology for solving SPPs. В статті розглядається задача оптимального розміщення куль різної розмірності в контейнерах довільних геометричних форм. Згідно з міжнародною класифікацією ця задача належить до класу SPP (Sphere Packing Problems). Метою даної роботи є створення єдиної методології розв’язання задач SPP. В статье рассматривается задача оптимального размещения шаров различной размерности в контейнерах произвольных геометрических форм. Согласно международной классификации эта задача относится к классу SPP (Sphere Packing Problems). Целью данной работы является создание единой методологии решения задач SPP. en Інстиут проблем машинобудування ім. А.М. Підгорного НАН України Проблеми машинобудування Applied mathematics Methodology to Solve Multi-Dimentional Sphere Packing Problems Методологія розв’язання задач розташування багатовимірних куль Методология решения задач размещения многомерных шаров Article published earlier |
| institution |
Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| collection |
DSpace DC |
| title |
Methodology to Solve Multi-Dimentional Sphere Packing Problems |
| spellingShingle |
Methodology to Solve Multi-Dimentional Sphere Packing Problems Yaskov, G.N. Applied mathematics |
| title_short |
Methodology to Solve Multi-Dimentional Sphere Packing Problems |
| title_full |
Methodology to Solve Multi-Dimentional Sphere Packing Problems |
| title_fullStr |
Methodology to Solve Multi-Dimentional Sphere Packing Problems |
| title_full_unstemmed |
Methodology to Solve Multi-Dimentional Sphere Packing Problems |
| title_sort |
methodology to solve multi-dimentional sphere packing problems |
| author |
Yaskov, G.N. |
| author_facet |
Yaskov, G.N. |
| topic |
Applied mathematics |
| topic_facet |
Applied mathematics |
| publishDate |
2019 |
| language |
English |
| container_title |
Проблеми машинобудування |
| publisher |
Інстиут проблем машинобудування ім. А.М. Підгорного НАН України |
| format |
Article |
| title_alt |
Методологія розв’язання задач розташування багатовимірних куль Методология решения задач размещения многомерных шаров |
| description |
This paper discusses the problem of optimally packing spheres of various dimensions into containers of arbitrary geometrical shapes. According to the international classification, this problem belongs to Sphere Packing Problems (SPPs). The aim of this work is to create an integrated methodology for solving SPPs.
В статті розглядається задача оптимального розміщення куль різної розмірності в контейнерах довільних геометричних форм. Згідно з міжнародною класифікацією ця задача належить до класу SPP (Sphere Packing Problems). Метою даної роботи є створення єдиної методології розв’язання задач SPP.
В статье рассматривается задача оптимального размещения шаров различной размерности в контейнерах произвольных геометрических форм. Согласно международной классификации эта задача относится к классу SPP (Sphere Packing Problems). Целью данной работы является создание единой методологии решения задач SPP.
|
| issn |
0131-2928 |
| url |
https://nasplib.isofts.kiev.ua/handle/123456789/158842 |
| citation_txt |
Methodology to Solve Multi-Dimentional Sphere Packing Problems / G.N. Yaskov // Проблеми машинобудування. — 2019. — Т. 22, № 1. — С. 67-75. — Бібліогр.: 37 назв. — англ, рос. |
| work_keys_str_mv |
AT yaskovgn methodologytosolvemultidimentionalspherepackingproblems AT yaskovgn metodologíârozvâzannâzadačroztašuvannâbagatovimírnihkulʹ AT yaskovgn metodologiârešeniâzadačrazmeŝeniâmnogomernyhšarov |
| first_indexed |
2025-11-25T23:28:48Z |
| last_indexed |
2025-11-25T23:28:48Z |
| _version_ |
1850581609290399744 |
| fulltext |
APPLIED MATHEMATICS
ISSN 0131–2928. Journal of Mechanical Engineering, 2019, vol. 22, no. 1
УДК 519.85
МЕТОДОЛОГИЯ
РЕШЕНИЯ ЗАДАЧ
РАЗМЕЩЕНИЯ
МНОГОМЕРНЫХ
ШАРОВ
Г. Н. Яськов,
канд. техн. наук
yaskov@ipmach.kharkov.ua
ORCID: 0000-0002-1476-1818
Институт проблем
машиностроения
им. А. Н. Подгорного
НАН Украины,
61046, Украина, г. Харьков,
ул. Пожарского, 2/10
В статье рассматривается задача оптимального размещения шаров различной раз-
мерности в контейнерах произвольных геометрических форм. Согласно международ-
ной классификации эта задача относится к классу SPP (Sphere Packing Problems).
Она заключается в размещении набора шаров (кругов, гипершаров) заданных радиусов
в контейнере с заданными метрическими характеристиками. Целью данной работы
является создание единой методологии решения задач SPP. Приведены основные по-
становки задачи: в виде задачи о рюкзаке; задачи с переменным размером контейнера
и соответствующие математические модели. На выбор стратегии решения влияют
вид постановки задачи, размерность пространства, в котором размещаются шары,
метрические особенности шаров (равные или неравные), их количество, геометриче-
ская форма контейнера, наличие технологических ограничений, ограничение на время
счета. Структурными элементами методологии являются математические модели,
способы построения начальных размещений, методы локальной и глобальной оптими-
зации. При разработке метода решения используется построение начальных допус-
тимых размещений случайным, решетчатым способами, с помощью жадного алго-
ритма и путем решения вспомогательной задачи нелинейного программирования. В
качестве методов локальной оптимизации рассматриваются модификации метода
возможных направлений, метод внутренней точки, метод множителей Лагранжа и
метод оптимизации по группам переменных. Для глобальной оптимизации использу-
ются метод перебора подмножеств шаров из заданного набора, метод перебора
крайних точек области допустимых решений, реализованные с помощью алгоритма
ветвей и границ, модификации методов сужающихся окрестностей, метод плавного
перехода из одного локального минимума в другой на основе увеличения размерности
задачи путем введения дополнительных переменных метрических характеристик,
метод решения, реализованный в виде последовательности задач нелинейного про-
граммирования возрастающей размерности, метод мультистарта. Предложены
стратегии решения задачи SPP для различных ее постановок.
Ключевые слова: шар, гипершар, упаковка шаров, задача о рюкзаке, задача с пе-
ременным размером, нелинейная оптимизация.
Введение
Задачи оптимального размещения геометрических объектов являются задачами геометриче-
ского проектирования [1]. Согласно международной классификации они относятся к классу задач
раскроя и упаковки (C&P – cutting and packing) [2]. Одной из наиболее распространенных задач C&P
является задача размещения шаров (SPP – sphere packing problem) различной размерности (2D – кру-
ги, 3D – шары, nD – гипершары). Она заключается в размещении шаров из набора с заданными ра-
диусами в заданном контейнере. При этом необходимо либо получить максимальный коэффициент
заполнения контейнера, либо найти минимальный возможный размер контейнера.
Задачи SPP имеют множество научных и практических применений, например, в текстильной,
швейной, автомобильной, аэрокосмической и химической промышленности [3], в ядерной энергетике при
моделировании процессов в ядерном реакторе [4], в аддитивном производстве для оптимизации геомет-
рической формы детали [5] и в медицине для планирования автоматизированного радиохирургического
лечения [6]. Размещения гипершаров используются для моделирования геометрии кристаллических со-
стояний [7]. Также такие задачи возникают при численной оценке интегралов либо на поверхности сфе-
ры, либо внутри нее [8]. Основным приложением задач размещения гипершаров является теория кодиро-
вания, цифровая связь и хранение информации, например, CD, сотовые телефоны и Интернет [8, 9].
В зависимости от постановки задачи различают два основных класса задач SPP: задача о рюк-
заке (KP – knapsack problem) и задача с переменным размером контейнера (ODP – open dimension
problem) [2]. Задача KP состоит в размещении шаров из заданного набора в контейнере заданных
фиксированных размеров с максимальным коэффициентом заполнения. В задаче ODP все шары из
Г. Н. Яськов, 2019
ПРИКЛАДНА МАТЕМАТИКА
ISSN 0131–2928. Проблеми машинобудування, 2019, Т. 22, № 1
заданного набора необходимо разместить в контейнере, все размеры которого фиксированы, кроме
одного, значение которого нужно минимизировать.
Часть исследователей в своих работах [10–14] используют формулировку SPP в виде KP. При
необходимости перейти к задаче ODP используется дихотомический поиск минимального размера
контейнера, т.е. решается последовательность задач KP. При этом для поиска начальных точек ис-
пользуются достаточно эффективные эвристические жадные алгоритмы (ЖА) [10, 11, 13, 14].
В работе [15] задача SPP формулируется как ODP. В качестве минимального размера контейне-
ра используются радиус, длина, высота, периметр, площадь, объем, площадь поверхности. Предложена
нелинейная математическая модель с дважды дифференцируемыми функциями, которая позволяет по-
лучить локальный минимум задачи. Начальные размещения шаров выбираются случайным образом.
В [16, 17] для решения задачи SPP, сформулированной в виде KP, используется математическая
модель, основанная на увеличении размерности задачи. Предполагается, что радиусы шаров временно
становятся переменными и максимизируется сумма объемов (площадей в 2D). Оптимизационный про-
цесс продолжается до тех пор, пока радиусы шаров не достигнут своих исходных значений.
Идея введения дополнительных переменных была использована и для задачи ODP для размещения
неравных шаров [18, 19]. Дополнительные переменные улучшают распределение шаров в контейнере: там,
где есть неиспользованное место, радиус шара может быть увеличен. Осуществляется переход из одной
точки локального минимума в другую, с лучшим значением функции цели. При работе такого алгоритма
важно найти связь между исходной и вспомогательной задачей с дополнительными переменными.
При небольшом количестве шаров можно использовать методы полного перебора экстре-
мальных точек задачи и теоретически получить глобальный экстремум. Однако на практике это за-
трудняется сложностью решения систем нелинейных уравнений [20].
При увеличении количества размещаемых шаров размерность задачи растет пропорционально,
а количество ограничений нелинейно. При большой размерности для решения задачи используются
либо случайные размещения шаров [21], либо приближенные методы, основанные на декомпозиции
задачи, например, с помощью метода оптимизации по группам переменных (МОГП) или ЖА [22, 23].
Таким образом, при различных постановках задачи размещения используются различные
подходы и методы их решения. Однако до сих пор нет единой методологии решения таких задач.
Поэтому целью этой работы является создание методологии решения задач размещения мно-
гомерных шаров.
Математические модели задачи размещения многомерных шаров
Очевидно, что в первую очередь на выбор стратегии решения влияет вид постановки задачи.
Рассмотрим постановки задач KP и ODP и их математические модели более подробно.
Задача KP формулируется следующим образом. Пусть имеется контейнер dС R⊂ заданной
геометрической формы с фиксированными размерами и набор шаров },...,2,1{ ,)( NIiuS N
d
ii =∈⊂ R ,
с заданными радиусами, где iu – вектор трансляции шара Ni IiS ∈ , , 2≥d
– размерность пространст-
ва. Задача заключается в размещении шаров из набора Nii IiuS ∈ ),( (всех или части из них) без
взаимных наложений в контейнере С с максимальным коэффициентом заполнения.
Предположим, что имеется K размеров шаров из набора Nii IiuS ∈ ),( с радиусами kr ,
},...,2,1{ KIk K =∈ . Обозначим количество шаров с радиусом kr , которые могут быть размещены
внутри контейнера С как kn , KIk ∈ , и сгенерируем кортеж ),....,,( 21 Knnnt = . В кортеже t должен
содержаться, по крайней мере, один ненулевой элемент kn , KIk ∈ . Обозначим множество всех
возможных кортежей t как T . Мощность множества T
∏
=
+
K
k
kn
1
)1( .
Сформируем подмножество t
S шаров )( jj uS , },...,2,1{ nJj =∈ из набора Nii IiuS ∈ ),( в
соответствии с кортежем ),....,,( 21 Knnnt = , где
APPLIED MATHEMATICS
ISSN 0131–2928. Journal of Mechanical Engineering, 2019, vol. 22, no. 1
∑
=
=
K
k
knn
1
.
Подмножество t
S с учетом трансляции обозначается }),({)( JjvEvS jj
t ∈= , где
dn
nvvvv R∈= ),...,,( 21 – вектор параметров размещения ),( jj vS ,Jj ∈ ),...,,( 21 djjjj xxxv =
– вектор
трансляции шара )( jj vS .
Задача KP может быть сформулирована следующим образом. Найти подмножество шаров
*)(vS
t , QIq ∈ , которое может полностью разместиться в контейнере С с максимальным
коэффициентом заполнения
∑
∈
∈
==
Jj
d
j
Tt
rtvFF )(max),( *** при условии TWtv
dn ×⊂∈ R),( , (1)
} ,0)( ,,, ,0),(:),{( JivjiJjivvtvW iijiij ∈≥Φ>∈≥Φ= , (2)
где )( ii vΦ – Ф-функция для t
i SS ∈ и объекта CRC int\
2* = [24,25],
∑
=
≥+−−Φ
d
k
jikjkijiij rrxxvv
1
22
0)()(=),( .
Решение задачи (1)–(2) может быть сведено к перебору элементов множества T , для которых
может быть найдена точка Wv ∈ . Поиск точки Wv ∈ может быть выполнен с помощью решения
следующей оптимизационной задачи нелинейного программирования, в которой радиусы шаров
являются переменными:
∑
∈
=κ=κ
Jj
jrrv max),( ***
при условии Mrv ∈),( , (3)
где ),...,,( 21 nrrrr = ;
{ jiJjirrvvrvM jijiij
nd <∈,≥,,Φ:∈= +
,, 0),(),(
)1(
R (4)
0),( ≥Φ iii rv , Ji ∈ , 0≥− irr , }Ji ∈ ,
2
1
2
)(),( ji
n
k
kjkijijiij rrxxrrvv +−−=,,Φ ∑
=
.
Задача (3)–(4) является многоэкстремальной. Функция цели линейная, поэтому экстремумы на-
ходятся в крайних точках области допустимых решений M (ОДР). Функция ),( jijiij rrvv ,,Φ является
квадратичной формой. Вид функции ),( iii rvΦ зависит от геометрической формы контейнера.
Теперь рассмотрим задачу ODP. Набор шаров Nii IiuS ∈ ),( , необходимо разместить в контей-
нере С с минимальным размером (площадь, объем, или метрическая характеристика).
Математическая модель задачи ODP [16] представляется как
)(min* µ=µ f при условии λ+⊂∈µ,= dN
WuY R
~
)( , (5)
где µ – переменная метрическая характеристика (вектор переменных метрических характеристик); λ
– количество переменных метрических характеристик ( 1=λ – в случае линейной характеристики,
например, длина, высота и т.д., и 2≥λ , если размер контейнера задается несколькими метрическими
характеристиками, например, площадь, площадь поверхности, объем); )(µf – функция, определяю-
щая переменный размер контейнера С ( µ=µ)(f для 1=λ );
ПРИКЛАДНА МАТЕМАТИКА
ISSN 0131–2928. Проблеми машинобудування, 2019, Т. 22, № 1
{ }NiiNjiij
dN
IiuIjiuuYW ∈≥µΦ∈<<≥Φ∈= λ+
,0),(,0,0),(:
~
R ; (6)
∑
=
+−−Φ
d
k
jikjkijiij rrxxuu
1
22
.)()(=),(
Задача (5)–(6) является многоэкстремальной задачей нелинейного программирования.
Методология решения задачи
Методология решения оптимизационных задач размещения шаров основана на анализе
постановки задачи, исходных данных и ограничений. Она включает в себя построение
математических моделей, которые охватывают класс задач размещения шаров, исследование их
особенностей и выработку стратегий решения поставленной задачи. Изучается структура, логические
связи, предлагаемые методы и средства решения задачи.
На выбор стратегии решения влияют следующие основные факторы:
– вид постановки задачи (KP или ODP);
– размерность пространства, в котором размещаются шары ( 2=d , 3=d и 4≥d );
– метрические особенности шаров (равные или неравные);
– количество размещаемых шаров;
– геометрическая форма контейнера;
– наличие технологических ограничений;
– ограничение на время счета.
В зависимости от особенностей постановки задачи и математической модели предлагается
стратегия решения задачи, структурными элементами которой являются: способы построения
допустимых размещений (начальных точек или приближенных решений), методы локальной
оптимизации и методы глобальной оптимизации.
Предлагается использовать следующие способы построения допустимых размещений:
случайный, решетчатый, жадный алгоритм, вспомогательная задача нелинейного программирования.
Случайный способ заключается в случайном выборе координат центров шаров и проверке
допустимости размещения [16]. В решетчатом способе центры шаров совпадают с узлами решетки
[16, 23, 26]. При использовании ЖА, который является модификацией МОГП, происходит
декомпозиция задачи на подзадачи [27]. В качестве группы переменных выбираются координаты
одного шара, а функция цели может быть выбрана эвристически. Для получения допустимого
размещения может также использоваться вспомогательная задача нелинейного программирования
или последовательность таких задач [26]. В некоторых случаях целесообразно применить
комбинации этих способов [16, 26].
В качестве методов локальной оптимизации в зависимости от особенностей математической
модели и количества размещаемых шаров следует использовать модификации метода возможных
направлений (МВН) [16, 23, 26], метод внутренней точки (МВТ) [18, 19, 28, 29], метод множителей
Лагранжа (ММЛ) [30], МОГП [30]. Во всех методах целесообразно применять стратегию активного
набора ограничений [31], благодаря которой значительно сокращаются вычислительные затраты.
МВН дает возможность свести решение задачи нелинейного программирования к
последовательности задач линейного программирования. Для задач SPP учитывается специфика
ограничений: часть из них может быть линейными, матрицы первых и вторых производных являются
сильно разреженными. Данные особенности делают возможность применить специальные пакеты
программ [32, 33] и решать задачи достаточно большой размерности (5000–10000) переменных. При
этом более устойчиво такие программы работают для шаров меньшей размерности 3≤d .
МВТ разработан для решения задач нелинейного программирования и эффективно работает
для задач средней размерности (до 1000 переменных).
ММЛ применяется для поиска локального экстремума в комбинации с методом спуска и
основан на анализе множителей Лагранжа активных ограничений.
В случае размещения большого количества шаров (при количестве переменных больше 10000)
следует использовать МОГП, выбирая по 1000–10000 переменных в группе. Для получения быстрого
результата при ограничении на время счета можно выполнить локальную оптимизацию размещения
APPLIED MATHEMATICS
ISSN 0131–2928. Journal of Mechanical Engineering, 2019, vol. 22, no. 1
каждого шара отдельно, выбрав группу из d переменных, задающих координаты центра шара. Такой
способ размещения называется еще методом последовательно-одиночного размещения [1] (sequential
addition [22, 34]).
Методы глобальной оптимизации представлены: методом перебора подмножеств шаров из
заданного набора (МПШ), методом перебора крайних точек ОДР (МПКТ) на основании алгоритма
ветвей и границ [20]; модификациями методов сужающихся окрестностей (МСО) [35, 36]; методом
плавного перехода (МПП) из одного локального минимума в другой на основе увеличения
размерности задачи путем введения дополнительных переменных метрических характеристик [18, 19,
29], методом решения, реализующим последовательность задач нелинейного программирования
возрастающей размерности (МРПЗ) [16, 26]. Также может быть использована комбинация указанных
методов. Для всех перечисленных методов применим метод мультистарта (ММС), который позволяет
расширить выбор возможных вариантов размещения.
МПШ используется для решения задачи KP. Для получения решения осуществляется перебор
различных вариантов выбора подмножества шаров из заданного набора, который реализуется в виде
дерева [37]. Для сокращения количества рассматриваемых вершин дерева используются правила
отсечения, с помощью которых отбрасываются неперспективные вершины на основе анализа нижней
и верхней оценок функции цели.
В МПКТ рещаются все подсистемы из системы ограничений, т.е. исследуются все поверхности,
описывающие границу ОДР. Для этого используется схема метода ветвей и границ. В случае линейной
функции цели решение находится в одной из крайних точек ОДР. Если в задаче все ограничения
обратно-выпуклые (например, размещение гипершаров в гиперпараллелепипеде), то крайняя точка
определяется системой уравнений, количество которых равно количеству переменных задачи.
С помощью модификаций МСО осуществляется оптимизация функции цели на различных
перестановках шаров и дереве решений на основе вероятностных свойств функции цели.
Теперь рассмотрим основные постановки и стратегии решения задачи.
Задача KP. В случае, если рассматривается задача, сформулированная в виде KP, то для
каждого выбранного варианта необходимо решить задачу (3)–(4) с помощью МВН и МВТ. При
небольшом количестве шаров и размерности шаров 3≤d возможен и полный перебор крайних точек
ОДР (МПКТ). Если количество шаров больше 10, то используется либо усеченное дерево решений,
либо МРПЗ. При большом количестве шаров (больше 10000 переменных) используются различные
модификации МОГП.
Способ получения начального размещения зависит от размерности задачи, соотношения
размеров шаров и контейнера и вида контейнера. Если решается задача размещения равных шаров,
размер контейнера значительно больше размера шаров, а размерность задачи не больше четырех, то
применяется решетчатый метод построения начального размещения. В других случаях используется
либо случайный способ, либо специальный ЖА [26], например, при пошаговом процессе МРПЗ.
Если контейнер имеет сложную геометрическую форму (например, при наличии большого
количества зон запрета), то для получения начальной точки используется вспомогательная задача
нелинейного программирования, решение которой позволяет восстановить допустимость из любой
случайно выбранной точки, не принадлежащей контейнеру [26].
Наличие технологических ограничений, например ограничений на минимально и
максимально допустимые расстояния, сужает область допустимых решений, а следовательно, и
количество вариантов размещения. Обычно в этом случае не применяется решетчатый способ
получения начальных размещений.
Задача ODP. Для решения задач с переменным размером контейнера, сформулированных в
виде (5)–(6), стратегия решения также зависит от числа размещаемых шаров и размерности задачи.
Обычно в задачах ODP контейнер с зонами запрета не рассматривается.
Если количество шаров не больше 10, а размерность шаров 3≤d , то используется МПКТ,
реализованный в виде модификации алгоритма ветвей и границ [20]. Для решения систем уравнений
используется метод Ньютона. В этом случае множество крайних точек включает в себя множество
локальных экстремумов и нет необходимости применять метод локальной оптимизации.
ПРИКЛАДНА МАТЕМАТИКА
ISSN 0131–2928. Проблеми машинобудування, 2019, Т. 22, № 1
При увеличении количества шаров или увеличении размерности пространства необходимо
использовать методы локальной оптимизации (МВН, МВТ), МСО и МПП, которые хорошо работают
для задач средней размерности (10–300) шаров. Для задач большой размерности в качестве метода
локальной оптимизации применим только алгоритм с использованием МВН, состоящий в решении
последовательности задач линейного программирования (5000–10000 переменных). Для получения
приближенного решения задачи необходимо применять МОГП. Этот же метод может быть
использован для получения приближенных решений.
В качестве способа получения начальных размещений используется решетчатый (для
одинаковых шаров) и ЖА (для разных шаров).
Следует отметить, что задача KP может быть решена как задача ODP, в которой
минимизируется коэффициент гомотетии контейнера. Такой переход целесообразен при решении
задач размещения неравных шаров. Он позволяет применять МПП для решения задач KP.
Выводы
Предложена единая методология решения задач размещения многомерных шаров.
Методология является развитием теории геометрического проектирования и может использоваться
специалистами в этой области для выбора стратегии решения задачи. С помощью разработанной
методологии имеется возможность решать задачи размещения шаров, сформулированные как KP, так
и ODP. Методология ориентирована на современные разработки в области геометрического
проектирования и на использование мощных пакетов программ для решения задач линейного и
нелинейного программирования.
Литература
1. Стоян Ю. Г., Яковлев С. В. Математические модели и оптимизационные методы геометрического проек-
тирования. Киев: Наук. думка, 1986. 268 с.
2. Waescher G., Haussner H. An improved typology of cutting and packing problems. Europ. J. Operational Re-
search. 2007. Vol. 183. Iss. 5. P. 1109–1130. https://doi.org/10.1016/j.ejor.2005.12.047
3. Hifi M., M'Hallah R. A. literature review on circle and sphere packing problems: models and methodologies. Ad-
vances in Operations Research. 2009. Vol. 2009. P. 1–19. doi: https://doi.org/10.1155/2009/150624
4. Snoj L., Ravnik M. Effect of packing fraction variations on the multiplication factor in pebble-bed nuclear reactors.
Kerntechnik. 2006. Vol. 71. No. 4. P. 208–213. https://doi.org/10.3139/124.100295
5. Leary M., Merli L., Torti F., Mazur M., Brandt M. Optimal topology for additive manufacture: A method for ena-
bling additive manufacture of support-free optimal structures. Materials and Design. 2014. Vol. 63. P. 678–690.
https://doi.org/10.1016/j.matdes.2014.06.015
6. Blyuss O., Koriashkina L., Kiseleva E., Molchanov R. Optimal Placement of Irradiation Sources in the Planning of
Radiotherapy: Mathematical Models and Methods of Solving. Computational and Math. Methods in Medicine.
2015. Vol. 2015. P. 1–8. https://doi.org/10.1155/2015/142987
7. Agapie S. C., Whitlock P. A. Random packing of hyperspheres and Marsaglia's parking lot test. Monte Carlo
Methods and Appl. 2010. Vol. 16. Iss. 3–4. P. 197–209. doi: https://doi.org/10.1515/mcma.2010.019
8. Conway J. H., Sloane N. J. A. Sphere Packings, Lattices, and Groups. New York: Springer-Verlag, 1999. 706 p.
https://doi.org/10.1007/978-1-4757-6568-7
9. Torquato S., Stillinger, F. H. Exactly solvable disordered sphere-packing model in arbitrary-dimensional Euclidean
spaces. Physical Review E. 2006. Vol. 73. Iss. 3. P. 031–106. https://doi.org/10.1103/PhysRevE.73.031106
10. Hifi M., Yousef L. A dichotomous search-based heuristic for the three-dimensional packing problem. Cogent Eng..
2015. Vol. 1. P. 1–15. https://doi.org/10.1080/23311916.2014.994257
11. Zeng Z., Yu X., He K., Huang W., Fu Z. Iterated Tabu Search and Variable Neighborhood Descent for packing un-
equal circles into a circular container. Europ. J. Operational Research. 2016. Vol. 250. Iss. 2. P. 615–627.
https://doi.org/10.1016/j.ejor.2015.09.001
12. Sutou A., Dai Y. Global Optimization Approach to Unequal Global Optimization Approach to Unequal Sphere
Packing Problems in 3D. J. Optimization Theory and Appl. 2002. Vol. 114. Iss. 3. P. 671–694. doi:
https://doi.org/10.1023/A:1016083231326
13. Zeng Z. Z, Huang W. Q., Xu R. C., Fu Z. H. An algorithm to packing unequal spheres in a larger sphere. Advanced Mate-
rials Research. 2012. Vol. 546–547. P. 1464–1469. https://doi.org/10.4028/www.scientific.net/AMR.546-547.1464
14. Kubach T., Bortfeldt A., Tilli T., Gehring H. Greedy Algorithms for Packing Unequal Spheres into a Cuboidal Strip
or a Cuboid. Asia-Pacific J. Operational Research. 2011. Vol. 28. No. 06. P. 739–753.
https://doi.org/10.1142/S0217595911003326
APPLIED MATHEMATICS
ISSN 0131–2928. Journal of Mechanical Engineering, 2019, vol. 22, no. 1
15. Birgin E. G., Sobral F. N. C. Minimizing the object dimensions in circle and sphere packing problems. Computers
& Operations Research. 2008. Vol. 35. Iss. 7. P. 2357–2375. https://doi.org/10.1016/j.cor.2006.11.002
16. Stoyan Yu., Yaskov G. Packing congruent hyperspheres into a hypersphere. J. Global Optimization. 2012. Vol. 52.
Iss. 4. P. 855–868. https://doi.org/10.1007/s10898-011-9716-z
17. Яковлев С. В., Яськов Г. Н., Коробчинский К. П. О методах переменного радиуса в задаче упаковки шаров в кон-
тейнеры. Питання прикл. математики i мат. моделювання: зб. наук. пр. Дніпро: ЛІРА. 2017. № 17. С. 265–272.
18. Стоян Ю. Г., Шайтхауер Г., Яськов Г. Н. Упаковка неравных шаров в различные контейнеры. Кибернетика
и систем. анализ. 2016. № 3. С. 97–105.
19. Yaskov G. N. Packing non-equal hyperspheres into a hypersphere of minimal radius. Проблемы машинострое-
ния. 2014. Т. 17. № 1. С. 48–53.
20. Stoyan Y., Yaskov G. Mathematical model and solution method of optimization problem of placement of rectan-
gles and circles taking into account special constraints. Intern. Transactions in Operational Research. 1998. Vol. 5.
No. 1. P. 45–57. https://doi.org/10.1111/j.1475-3995.1998.tb00101.x
21. Visscher W. M., Bolsterli M. Random Packing of Equal and Unequal Spheres in Two and Three Dimensions. Na-
ture. 1972. Vol. 239. P. 504–507. https://doi.org/10.1038/239504a0
22. Mueller G. E. Numerically packing spheres in cylinders. Powder Technology. 2005. Vol. 159. Iss. 2. P. 105–110.
https://doi.org/10.1016/j.powtec.2005.06.002
23. Яськов Г. Н. Упаковка большого числа конгруэнтных кругов в цилиндре. Доп. НАН України. 2009. № 12. С. 43–49.
24. Chernov N., Stoyan Yu., Romanova T. Mathematical model and efficient algorithms for object packing problem.
Computational Geometry. 2015. Vol. 43. Iss. 5. P. 535–553. https://doi.org/10.1016/j.comgeo.2009.12.003
25. Stoyan Y., Pankratov A., Romanova T. Quasi-phi-functions and optimal packing of ellipses. J. Global Optimiza-
tion. 2016. Vol. 65. Iss. 2. P. 283–307. https://doi.org/10.1007/s10898-015-0331-2
26. Stoyan Y., Yaskov G. Packing equal circles into a circle with circular prohibited areas. Intern. J. Computer Mathe-
matics. 2012. Vol. 89. Iss. 10. P. 1355–1369. https://doi.org/10.1080/00207160.2012.685468
27. Яськов Г. Н. Метод решения задачи размещения разных кругов с перспективным выбором начальных то-
чек. Зб. наук. праць Харк. ун-ту повітряних сил. 2010. Вип. 3. №. 25. С. 119–122.
28. Wächter A., Biegler L. T. On the implementation of a primal-dual interior point filter line search algorithm for
large-scale nonlinear programming. Math. Programming. 2006. Vol. 106. Iss. 1. P. 25–57.
https://doi.org/10.1007/s10107-004-0559-y
29. Stoyan Yu., Yaskov G. Packing unequal circles into a strip of minimal length with a jump algorithm. Optimization
Letters. 2014. Vol. 8. Iss. 3. P. 949–970. https://doi.org/10.1007/s11590-013-0646-1
30. Yaskov G. N. Random packing of identical spheres into a cylindrical container. Системи обробки інформації.
2008. Вип. 1. №. 68. С. 128–130.
31. Stoyan Y., Chugay A. Packing cylinders and rectangular parallelepipeds with distances between them into a given
region. Europ. J. Operational Research. 2009. Vol. 197. Iss. 2. P. 446–455. https://doi.org/10.1016/j.ejor.
2008.07.003
32. Gondzio J. HOPDM (version 2.12) – a fast LP solver based on a primal-dual interior point method. Europ. J. Opera-
tional Research. 1995. Vol. 85. Iss. 1. P. 221–225. https://doi.org/10.1016/0377-2217(95)00163-K
33. Meszaros C. On numerical issues of interior point methods. SIAM J. on Matrix Analysis and Appl. 2008. Vol. 30.
Iss. 1. P. 223–235. https://doi.org/10.1137/050633354
34. Torquato S., Uche O. U., Stillinger, F. H. Random sequential addition of hard spheres in high Euclidean dimen-
sions. Physical Review E. 2006. Vol. 74. Iss. 6. P. 061–308. https://doi.org/10.1103/PhysRevE.74.061308
35. Stoyan Y., Scheithauer G., Yaskov G. Packing of various radii solid spheres into a parallelepiped. Central Europ. J.
Operations Research. 2003. Vol. 11. No. 4. 2003. P. 389–407.
36. Stoyan Y., Yaskov G. Packing Identical Spheres into a Rectangular Parallelepiped. In: Bortfeldt A., Homberger J.,
Kopfer H., Pankratz G., Strangmeier R. (eds) Intelligent Decision Support. Gabler, 2008. pp. 47–67.
https://doi.org/10.1007/978-3-8349-9777-7_4
37. Яськов Г. Н. Дерево решений для задачи размещения неравных шаров в шаре заданного радиуса. Інформацій-
ні системи і технології ІСТ-2017: пр 6-ї міжнар. наук.-техн. конф. (Коблеве, 11–16 вер. 2017). Коблеве, 2017.
С. 143–144.
Поступила в редакцию 14.01.2019
|