Рассматривается задача оптимизации упаковки гомотетичных одинаково ориентированных эллипсоидов в контейнере минимального объема. Строится математическая модель в виде задачи нелинейного программирования. Ограничения непересечения эллипсоидов и их включения в контейнер построены с использованием метода phi-функций В качестве контейнера рассматривается либо прямоугольный параллелепипед переменной длины, ширины и высоты, либо эллипсоид с переменным коэффициентом гомотетии. Предлагается алгоритм поиска локально оптимальных решений. с использованием гомотетических преобразований эллипсоидов и оптимизационной процедуры, позволяющей свести задачу с большим числом неравенств к последовательности задач с меньшим числом неравенств. Для поиска локальных минимумов задачи используется подход, в основе которого лежит метод мультистарта и оптимизационная процедура, включающая поиск допустимых стартовых точек и локальную оптимизацию. В качестве локально-оптимального решения выбирается наилучший из полученных локальных экстремумов. С целью минимизации числа нелинейных неравенств, формирующих область допустимых решений, предложена процедура LOFRT, которая позволяет значительно сократить вычислительные ресурсы. Приводятся результаты численных экспериментов.
The paper studies the packing problem of homothetic the same oriented ellipsoids into a container of minimal volume. The container can be a rectangular parallelepiped or an ellipsoid. We formulate the model in the form of a nonlinear programming problem. To constract the non-overlapping and containment constraints using of phi-function technique. We propose the efficient algorithm, which employes a homothetic transformation of ellipsoids and the optimization procedure Local Optimization with Feasible Region Transformation (LOFRT), which allow us to reduce considerably the dimension of the problem and computational time. Our algorithm also involves generating a number of random starting points. We choose the best local minimum as the solution of the problem. Our model can be realized by the current state-of-the art local or global solvers. A several computational results are provided.
