CC BY
84
6. Forrest, S. Self-Nonself Discrimination in a Computer / S. Forrest, A. S. Perelson, R. Cherukuri, L. Allen // Research in Security and Privacy : IEEE Symposium, Oakland, 16-18 May 1994 : proceedings. - Los Alamitos: IEEE, 1994. - P. 202-212.
7. D 'haeseleer, P. An immunological approach to change detection: algorithms, analysis, and implications / P. D’haeseleer, S. Forrest, P. Helman // Computer Security and Privacy : IEEE Symposium, Oakland, 6-8 May 1996 : proceedings. - Los Alamitos: IEEE, 1996. - P. 110-119.
8. Hofmeyr, S. Architecture for an artificial immune system /
S. Hofmeyr, S. Forrest // Evolutionary Computation. - 2000. -№ 8. - P. 443-473.
9. Balthrop, J. Revisiting LISYS: parameters and normal behavior / J. Balthrop, S Forrest, M. Glickman // Evolutionary Computation: congress CEC’02, Honolulu, 12-17 May 2002 : proceedings. - Los Alamitos: IEEE, 2002. - P. 1045-1050.
10. Balthrop, J. Coverage and generalization in an artificial immune system / J. Balthrop, F. Es-ponda, S. Forrest, M. Glickman / / Genetic and Evolutionary Computation: conference GECC0-2002, New York, 9-13 July 2002 : proceedings . -San Francisco: Morgan Kaufmann, 2002. - P. 3 10.
11. Farmer, J. The immune system, adaptation, and machine learning / J. Farmer, N. Packard, A. Perelson // Physica D : Nonlinear Phenomena. - 1986. - № 2. - P. 187-204.
12. Frank, A. SPECT heart data set [Electronic resource] / A. Frank, A. Asuncion. - Irvine : University of California, 2010. - Access mode : http://archive.ics.uci.edu/ml/datasets/ SPECT+Heart.
13. Ji, Z. Negative selection algorithms: from the thymus to V-detector : dissertation ... doctor of philosophy / Z. Ji. -Memphis : The University of Memphis, 2006. - 337 p.
14. Герасимчук, Т. С. Использование искусственных иммунных систем для прогнозирования риска развития рекку-рентных респираторных инфекций у детей раннего воз-
УДК [519.674:519.688]:62-11
раста / Т. С. Герасимчук, С. А. Зайцев, С. А. Субботин // Діагностика та лікування інфекційно опосередкованих соматичних захворювань у дітей: міжрегіональна науково-практична конференція, Донецьк, 10-11 лютого 2011 : матеріали. - Донецьк : Норд-пресс, 2011. - С. 27-29.
Стаття надійшла до редакції 04.05.2011. Зайцев С. О., Суботін С. О.
ПОБУДОВА ДІАГНОСТИЧНИХ МОДЕЛЕЙ НА ОСНОВІ ПАРАДИГМИ НЕГАТИВНОГО ВІДБОРУ З ВИКОРИСТАННЯМ ПРИНЦИПУ МАСКУВАННЯ ДЕТЕКТОРІВ Дістали подальшого розвитку модель негативного відбору з використанням маскованих детекторів та метод її навчання. Досліджувалися різноманітні критерії зупину для методу навчання моделі. Розроблено критерій зупину, що дозволяє виконувати своєчасний зупин процесу навчання, і таким чином прискорити даний процес, а також скоротити використання обчислювальних ресурсів.
Ключові слова: негативна селекція, маскований детектор, критерій зупину
Zaitsev S. A., Subbotin S. A.
THE DIAGNOSIS MODEL BUILDING ON THE BASIS OF NEGATIVE SELECTION PARADIGM USING THE PRINCIPLE OF DETECTOR MASKING
A new negative selection model using masked detectors with training method has been developed. Various stopping criteria for the training method have been investigated. The stopping criterion is proposed. It allows to terminate training procedure, to speedup training process and to reduce computational resources.
Key words: negative selection, masked detector, stopping criterion.
Чопоров С. В.
Магистр Запорожского национального университета
ПОСТРОЕНИЕ НЕРАВНОМЕРНЫХ ДИСКРЕТНЫХ СЕТОК ДЛЯ ФУНКЦИОНАЛЬНЫХ МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ НА БАЗЕ ТЕОРИИ Я-ФУНКЦИЙ
В работе рассмотрена проблема математического моделирования геометрических объектов на базе теории ^-функций. Предложен подход к построению неравномерных сеток шестигранных элементов.
Ключевые слова: ^-функция, сетка, шестигранный элемент, математическая модель.
ПОСТАНОВКА ПРОБЛЕМЫ в пространстве занимают определенный объем. Напри-
мер, сложные инженерные конструкции, сооружения, Исследование прочности и долговечности проекти- тт
„ „ машины, механизмы и прочее. По сути, речь идет о ре-
руемых инженерных конструкций является важной со- „ ^ „
шении так называемой обратной задачи аналитической ставляющей современной техники. В настоящее время „ „ -
геометрии, когда для имеющейся геометрической обла-
для этих целей широко применяется компьютерное мости строится ее аналитическое описание.
делирование, позволяющее заменить длительное и до- Множество физических объектов реального мира
рогостоящее испытание опытного образца изучением
можно определить как сплошные тела, геометрическое
соответствующих математических моделей. Следователь-
моделирование которых можно рассматривать как про-
но, актуальной является проблема построения матема-
^ - цесс формализации представления геометрии существу-
тических моделей геометрических объектов, которые
© Чопоров С. В., 2011
ющего или воображаемого объекта. Наиболее общим методом определения множества точек, принадлежащих объекту, является введение предиката, являющегося индикатором принадлежности объекту точки пространства (неявное определение). Простейшей формой такого предиката будет ограничение на знак некоторой действительной функции в виде
О = {(х, у,1) I Е(х, у, 2) > 0},
где неявная функция Е (х, у, 2) больше нуля внутри области, соответствующей телу О, равна нулю на границе и меньше нуля вне тела. Например, для построения мате-магической модели тела в форме эллипсоида с центром в точке (Хо, Уо, 2о) необходимо определить коэффициенты а, Ь и с в формуле
F (xvz )= 1 (x - x0)2 Х - у°)2 (z - z0)2
FellipsoidVx> v>z) 1 7 2 2 ■
abc
Математические модели геометрических объектов нетривиальной формы могут быть определены при помощи более сложных функций конструктивно, используя логические комбинации более простых функций. Такие действия с математическими функциями логически эквивалентны стандартным операциям над множествами. В работах Владимира Логвиновича Рвачева [1-3] разработан достаточно общий подход, получивший название теории ^-функций, в котором предложены принадлежащие классу Cm функции для описания теоретико-множественных операций. Наиболее распространенная на практике система ^-функций имеет вид (1). Например, тело в форме объединения двух шаров (рис. 1), один из которых радиуса r0 с центром в начале координат, второй - радиуса r с центром в точке (0, r0, 0) можно описать формулой (2).
C = const,
x = - x,
xi л x2 = xl + x2 - д/ x2 + x2 ,
xl V x2 = xl + x2 + )/x2 + x2 .
(1)
Среди конечных элементов, используемых в пространстве, наиболее распространенными являются элементы тетраэдрической и шестигранной форм. Преимуществом первого типа элементов является большая топологическая гибкость, второго - возможность уменьшения размерности систем линейных алгебраических уравнений за счет использования меньшего количества элементов (необходимо использовать пять тетраэдров для представления объема элемента, представленного одним шестигранником, рис. 2).
При этом функции формы линейного тетраэдрического конечного элемента имеют вид
(х,у,2) = ао + а1х + а2У + аз2, а^ е Я,1 = 0,3, (3)
а функции формы шестигранного конечного элемента -
Ын (х,у,2) = ао + а^х + а2У + а3 2 + а4>ху + а5 Х2 + а^уг + а^хуг,
аI е Я,1 = о,7. (4)
Наличие нелинейных членов в соотношении (4) приводит к тому, что градиенты шестигранного конечного элемента (в отличии тетраэдрического) не постоянны и изменяются вдоль одной из координатных плоскостей. Таким образом, применение шестигранных конечных элементов является более предпочтительным в вычислительном плане.
На практике в основном требуется построение нерегулярных сеток, где количество конечных элементов сгущается в местах, где расположены так называемые сингулярности (особенности) объекта (острые углы, отверстия, трещины и т. п.). С другой стороны, построение первоначального приближенного разбиения требует последующей его оптимизации.
Fballs (x У, z) = (r02 - x2 - y 2 - z 2 )v [r12 - x2 - (y - r0 )2 - z2 J .(2)
Рис. 1. Тело в форме объединения двух шаров: Г0 = 0,8,
Гі = 0,4
Поскольку математические модели геометрических объектов в основном используются при численном анализе различных параметров состояния сложных инженерно-технических объектов, на практике требуется построение конечно-элементных дискретных математических моделей. Для этого при использовании ^-функций, в первую очередь, требуется разработка методов построения дискретных моделей областей, описанных с помощью неявных функций.
Рис. 2. Представление шестигранника пятью тетраэдрами
Таким образом, разработка методов построения неравномерных дискретных моделей трехмерных объектов, заданных при помощи Л-функций, на базе шестигранных конечных элементов является сложной и актуальной задачей.
АНАЛИЗ ПУБЛИКАЦИЙ ПО ТЕМЕ ИССЛЕДОВАНИЯ
Методы построения нерегулярной сетки шестигранных конечных элементов можно разделить на прямые и непрямые. Прямые методы используют непосредственное построение шестигранных ячеек внутри заданной границы, в то время как непрямые базируются на преобразовании некоторого исходного дискретного представления области.
Среди прямых методов можно выделить фронтальные методы (наиболее распространенной тут идеей является исчерпание области пласт за пластом, формируя таким образом и дискретизацию границы, и дискретизацию внутренней части области [4]) и методы, использующие блочную декомпозицию [5-7] или суперпозицию [8-Ю].
Методы декомпозиции и блочной декомпозиции основаны на идее (полуавтоматической или автоматической) декомпозиции тела на кубоподобные блоки, для каждого из которых дискретизация может быть получена отображением на дискретизацию единичного куба, с последующим объединением дискретизаций частей. Для автоматизации процесса декомпозиции на блоки используются срединные поверхности. Достоинством данной группы методов является их высокая скорость и качество получаемой сетки. Основным недостатком - затруднительность их применения для дискретизации сложных нестандартных областей.
Основная идея методов суперпозиции заключается в использовании начальной сетки, которая может быть более или менее просто создана для части пространства, содержащей объект, и последующей адаптации сетки к границе области. Ключевым шагом данной группы методов является генерация сетки высокого качества около границы объекта. Для адаптации начальной сетки могут быть использованы техники изоморфизма [8, 9] и проекции. При всей прагматичности такого подхода результирующие алгоритмы данного типа методов является весьма трудоемкими для реализации [11].
Непрямые методы можно разделить на фронтальные, в основе которых идея послойной конвертации начальной сетки тетраэдров [12, 13], методы, использующие двойственное представление (8ТС-представление [14]) и методы поэлементной конвертации (например, на основе шаблонов). Преимуществом данной группы методов является отсутствие вырожденных элементов. Недостатком - необходимость построения качественной начальной дискретизации, что само по себе является трудоемкой задачей.
МЕТОД СУПЕРПОЗИЦИИ ДЛЯ ФУНКЦИОНАЛЬНОГО ПРЕДСТАВЛЕНИЯ НА БАЗЕ ТЕОРИИ й-ФУНКЦИЙ
Пусть тело О представлено неявной функцией на базе тории Л-функций. Использование методов теории Л-функций при неявном представлении тела О требует учета особенностей такого представления при разработке алгоритмов получения дискретного представления. Неявное представление дает правило для проверки принадлежности точки О, однако, не дает правило генерации системы точек и топологии элементов, образующих О. Это делает более предпочтительным использование в качестве алгоритмической базы метода суперпозиции совместно с сеточным методом и техникой изоморфизма [8-1о] на шаге адаптации пограничного слоя элементов. Тогда общий алгоритм получения дискретного представления тела заданного неявной функций можно представить схемой, приведенной на рис. 3.
На первом шаге алгоритма для формирования начальной сетки достаточно покрыть область равномерной сеткой с шагом к и определить значение Л-функции в узлах сетки. Затем последовательно отсеивать элементы, в узлах которых есть внешние точки или которые не соответствуют критерию близости к границе области (на
расстоянии А = о,5к есть внешние точки).
Граница полученной таким образом начальной сетки является многогранником, каждая грань которого является четырехугольником. Следовательно, ее можно рассматривать как неструктурированную сетку четырехугольных элементов, для которой может быть получена изоморфная (однозначно соответствующая) сетка на границе тела путем установления для каждого узла соответствующего граничного, а для каждого ребра границы начальной сетки - ребро на границе. Таким образом,
Рис. 3. Общая схема построения дискретной модели
каждому четырехугольнику будет соответствовать четырехугольная грань на границе тела, которые совместно определяют шестигранный элемент.
На втором шаге к каждому граничному узлу начальной сетки определяется нормаль как среднее арифметическое нормалей смежных в узле граней (рис. 4). Полученные нормали определяют направления для поиска соответствующих граничных точек на третьем шаге.
На четвертом шаге происходит учет характерных точек (геометрических особенностей, концентраторов напряжений и т. д.), так как утеря таковых может приводить к значительной утери точности при использовании полученного дискретного представления численными методами. При этом характерные точки могут быть учтены как при помощи техники перемещения необходимого граничного узла, так и при помощи коррекции направления нормали соответствующего узла начальной сетки и последующей репроэкции.
На последнем шаге алгоритма при помощи техники изоморфизма формируются элементы, соединяя узлы начальной сетки с соответствующими узлами на границе (рис. 5).
Следует отметить, что в результате применения описанного подхода генерируются элементы приблизительно одинакового размера, что позволяет применять для последующего улучшения качества сеток методы на основе локального сглаживания Лапласа.
УНИВЕРСАЛЬНЫЙ ШАБЛОН ДЛЯ ПОСТРОЕНИЯ НЕРАВНОМЕРНЫХ СЕТОК
Наиболее распространенным подходом к получению неравномерной начальной сетки для методов на основе суперпозиции является использование множества шаб-
Рис. 4. Опеделение нормалей к узлам начальной сетки
Рис. 5. Построение граничного слоя шестигранных элементов
лонов сгущения с таблицами соответствия шаблонов [8-1о]. Идея алгоритма заключается в рекурсивном разбиении области, в которую помещен исходный объект, на 27 элементов до тех пор, пока в узлах сетки удовлетворяет некоторая управляющая функция (пример такого сгущения на плоскости приведен на рис. 6), с последующим восстановлением топологии сетки. При этом в одном узле не должны быть смежными части, у которых разность уровней рекурсии больше единицы.
При построении неравномерной адаптивной сетки для трехмерной области с учетом вращения получается 22 различных конфигурации (256 без учета вращения, рис. 7), покрывающих все возможные ситуации взаимного расположения частей с разными уровнями рекурсии. Такое количество необходимых конфигураций делает трудным их построение, верификацию и программирование.
В работах [8, 9] предложено использовать упрощенное множество шаблонов, которое изображено на рис. 8.
В работе [1о] разработано альтернативное множество шаблонов (рис. 9), которое позволяет генерировать неравномерные сетки с меньшим использованием шаблона 8 (разбиения куба на 27 частей).
Однако и первое, и второе множество шаблонов в ряде случаев приводит к необходимости использования шаблона, соответствующего конфигурации 8 (разбиению на 27 элементов), что приводит к избыточному сгущению сетки. Для решения этой задачи, в таких случаях, можно использовать универсальный шаблон (рис. 1о), конфигурации граничных элементов которого будут соответствовать одному из шаблонов множества 2.
Для корректного определения координат внутренних узлов может быть использовано изопараметрическое преобразование. Такой подход может быть использован также для сгущения построенных дискретных моделей.
Рис. 6. Идея сгущения сетки
Рис. 7. Возможные конфигурации в трехмерном пространстве
Рис. 8. Множество шаблонов 1
Рис. 9. Множество шаблонов 2 ВЫВОДЫ Е-функций. Предложен подход к построению неравно-
Таким образом, в работе разработан подход к полу- мерных дискретных моделей на базе разработанного
чению дискретных моделей, соответствующих неявным универсального шаблона. Результаты работы реализо-
аналитическим моделям, построенным на базе теории ваны в виде подсистемы САПР, которая может быть ис-
пользована при моделировании и проектировании различных технических деталей и объектов (некоторые примеры приведены на рис. 11).
Е Н
в с в с
Рис. 10. Универсальный шаблон
Рис. 11. Примеры неравномерных сеток
СПИСОК ЛИТЕРАТУРЫ
1. Рвачев, В. Л. Теория .R-функций и некоторые ее приложения / В. Л. Рвачев. - К. : Наукова думка, 1982. - 552 с.
2. Рвачев, В. Л. Новые подходы к построению уравнений трехмерных локусов с помощью .R-функций / В. Л. Рвачев, А. В. Толок, Р. А. Уваров, Т. И. Шейко // Вісник Запорізького державного університету. - 2000. - № 2. - С. 119-130.
3. Рвачев, В. Л. Введение в теорию .R-функций / В. Л. Рвачев, Т. И. Шейко // Проблемы машиностроения. - 2001. -Т. 4, № 1-2. - С. 46-58.
4. Staten, M. L. Unconstrained Paving & Plastering: A New Idea for All Hexahedral Mesh Generation / M. L. Staten,
S. J. Owen, T. D. Blacker // Proceedings, 14th International Meshing Roundtable. - Sandia National Laboratories: Springer-Verlag, 2005. - P. 399-416.
5. Armstrong, C. G. Medials for Meshing and More / C. G. Armstrong, D. J. Robinson, R. M. McKeag, T. S. Li, S. J. Bridgett, R. J. Donaghy, C. A. McGleenan // Proceedings, 4th International Meshing Roundtable. - Sandia National Laboratories, 1995. - P. 277-288.
6. Holmes, D. J. Generalized Method of Decomposition Solid Geometry into Hexahedron Finite Elements / D. J. Holmes / / Proceedings, 4th International Meshing Roundtable. - Sandia National Laboratories, 1995. - P. 141-152.
7. Liu, S. -S. Automatic Hexahedral Mesh Generation by Recursive Convex and Swept Volume Decomposition / S.-S. Liu, R. Gadh // Proceedings, 6th International Meshing Roundtable. - Sandia National Laboratories, 1997. - P. 217-231.
8. Schneiders, R. A Grid-based Algorithm for the Generation of Hexahedral Element Meshes / R. Schneiders // Engineering with Computers. - 1996. - № 12. - P. 168-177.
9. Schneiders, R. Octree-based Generation of Hexahedral Element Meshes / R. Schneiders, R. Schindler, F. Weiler // 5th Annual International Meshing Roundtable. - 1996. - P. 205-216.
10. Ito, Y. Octree-based reasonable-quality hexahedral mesh generation using a new set of refinement templates / Yasushi Ito, Alan M. Shih, Bharat K. Soni // International Journal For Numerical Methods in Engineering. - 2009. - № 77 (13). - P. 1809-1833.
11. Thompson, J. F. Hand book of grid generation / Joe F. Thompson, Bharat K. Soni, Nigel P. Weatheril. - New York: CRC Press, 1999. - 1200 p.
12. Owen, S. J. H-Morph: An Indirect Approach to Advancing Front Hex Meshing / S. J. Owen, S. Saigal // International Journal for Numerical Methods in Engineering. - 2000. -№ 1 (33). - P. 289-312.
13. Owen, S. J. Hex-domain mesh generation using 3D constrained triangulation / S. J. Owen // Computer Aided Design. - 2001. -№3(33). - P. 211-220.
14. Murdoch, P The spatial twist continuum: A connectivity based method for representing all-hexahedral finite element meshes / P. Murdoch, S. E. Benzley, T. D. Blacker,
S. A. Mitchel // Finite Elements in Analysis and Design. -1997. - № 28. - P. 137-149.
Стаття надійшла до редакції 11.05.2011.
Choporov S. V.
GENERATION OF NONUNIFORM HEXAHEDRAL ELEMENT MESHES FOR FUNCTIONAL MODELS ON THE BASIS OF R-FUCTIONS
The problem of mathematical modeling of geometrical objects on the basis of R-functions is described in the article. Author propose an approach for generation of nonuniform hexahedral element meshes.
Key words: R-function, mesh, hexahedral element, mathematical model.
Чопоров С. В.
ПОБУДОВА НЕРІВНОМІРНИХ ДИСКРЕТНИХ СІТОК ДЛЯ ФУНКЦІОНАЛЬНИХ МАТЕМАТИЧНИХ МОДЕЛЕЙ НА БАЗІ ТЕОРІЇ R-ФУНКЦІЙ
В роботі розглянута проблема математичного моделювання геометричних об’єктів на базі теорії R-функцій. Запропоновано підхід до побудови нерівномірних сіток шестигранних елементів.
Ключові слова: R-функція, сітка, шестигранний елемент, математична модель.
