Введение к работе
Актуальность темы. Значительный рост производства электроэнергии, возросшие требования к экономичности, создание новых технологий приводят к росту единичных мощностей и номинального напряжения оборудования, что влечет повышение требований к надежности и другим технико-экономическим показателям при создании новых конструкций и аппаратов высоковольтной техники. Методы моделирования с использованием современных компьютеров часто оказываются единственным способом определения оптимальных конструкций, так как возможности их экспериментального исследования часто ограничены методическими и техническими трудностями.
Сложность расчета электростатических полей большинства электрофизических устройств заключается: І) в большой разномасштабное элементов этих конструкций; 2) в наличии большого числа изоляционных материалов с различной диэлектрической проницаемостью; 3) в наличии разнородных граничных условий заданных на сложных границах, 4) в наличии больших градиентов решения, что может привести к большой погрешности полученного численного решения.
Для решения задачи расчета полей с учетом перечисленных осо-беннностей широко применяется метод конечных элементов. Трудоемкость дискретизации области в методе конечных элементов для реальных задач часто превосходит затраты на получение численного решения. Этап дискретизации области является недостаточно изученным и требует разработки новых алгоритмов гостро-
ения сеток, отвечающих современным требованиям автоматизации и надежности.
Поэтому вопросы построения и обоснования эффективных алгоритмов дискретизации расчетных областей, хорошо отражающих основные физические свойства процессов, имеют важное значение для численного исследования рассматриваемых задач.
Целью диссертационной работы является исследование математических основ построения дискретизаций двумерных областей с учетом неравномерного распределения точек в области на основе разбиений Б.Н.Делоне и задаваемых функций'плотности, построение алгоритмов дискретизации трехмерных областей со сложной геометрией, разработка методов численного контроля корректности произвольных трехмерных сеток, практическая реализация разработанных алгоритмов и проведение численного моделирования двумерных и трехмерных электростатических полей на основе построенных дискретизаций.
Научная новизна изложенных в работе результатов заключается в еле душем:
Разработаны алгоритмы построения триангуляции Делоне с учетом функций плотности точек. ,
Введено понятие n-мерных динамических триангуляции Делоне и доказаны теоремы о гостроении таких триангуляции.
Предложен алгоритм построения трехмерных сеток на основе
локально-модифицированного подхода и поверхностных триангуляции
Делоне.
Разработаны, подходы и алгоритмы численного тестирования і
- 5.-
трхмерных сеток на основе геометрических, топологических и комбинаторных свойств.
Предложенные подходы и алгоритмы реализованы в разработанных программах: интерактивной системе 5tella-2D и ШЇЇІ РАМЗЕС-3, предназначенных, соответственно, для моделирования двумерных и трехмерных электростатических полей. На основе практического использования выработаны направления дальнейшего развития предложенных алгоритмов.
Практическая ценность. Интерактивная система Stella-2D позволяет эффективно проводить расчеты электростатических полей сложных электрофизических конструкций. Автоматический выбор точек сгущения позволяет получать сетки для метода конечных .элементов с учетом особенностей решения и геометрии расчетной области. Поэлементные итерационные методы позволяют рассчитывать сложные реальные задачи в ограниченной оперативной памяти. ПИП РАМЗЕС-3 позволяет строить и тестировать трехмерные сетки . сетки в сложных областях и проводить расчет электростатических полей в трехмерной постановке.
Система Stella-2D и ППП РАМЗЕС-3 Енедрены в следующих организациях СНГ: НПО "Исток" г.Москва, ФТЦ ИЯФ г.Липецк. ВНИИЭФ г.С.-Петербург, СибНИИЭ г.Новосибирск, ИГиЛ г.Новосибирск.
Публикации. Основные результаты, изложенные в диссертации, опубликованы в работах CI-I21.
Апробация работы. Основные результаты работы докладывались на семинарах ВЦ СО АН СССР г. Новосибирск, Сибирской Школе
по комплексам программ з г. Томске (1984г.) Шушенском (1986г.), Всесоюзных семинарах по методам построения сеток в г. Свердловске (1987г., 1930) и Челябинске-70 (1992г.), Всесоюзной конференции "Актуальные проблемы вычислительной и прикладной математики" в г.Новосибирске (1987г.), Всесоюзных семинарах "Методы расчета электронно- оптических систем" в г.Ташкенте (1988г.) и Львове (1990г.), Сибирской школе по вычислительной математике (1989г.), Дальневосточной школе по вычислительной математике (1989г.), Всесоюзном семинаре-совещании "Автоматизация проектирования и моделирования электронно- оптических систем" в г.Винница (1991г.), университетах за рубежом: г.Павиа (1990г., Италия, проф. Э.Маженес), г.ТюОинген (1990г., ФРГ, проф. Э.Каспер), Мэриленд (1992г., США, проф. И.Бабушка).
Структура к объем работы. Диссертация состоит из введения, трех глав, заключения, списка литературы и приложений. Объем содержательной части диссертации - 148 страниц, списка литературы - 10 страниц, приложений - 22 страниц. Работа Еключает 27 рисунков, список литературы из 94 наименований.
II. СОДЕРЖАНИЕ РАБОТЫ Во введении описывается общая формулировка задач построения дискретизаций для метода конечных элементов на примере решения эллиптических задач, описываются общие требования к типам областей, дискретизациям и элементам. Для алгоритмических походов к построению дискретизаций изложена классификация алгоритмов автоматической генерации неструктурных сеток.
Для подхода к построению сеток на основе разбиений Б.Н.Делоне I
приводится обзор фундаментальных и современных результатов.
Сформулирована задачи исследования и . описана структура диссертационной работы.
Первая глава посвящена изучению математических . основ разбиений Делоне, описанию критериев сравнения триангуляции, построению алгоритмов генерации нерегулярных триангуляции Делоне для двумерных областей и построению n-мерных динамических разбиений Делоне.
В _!.!. дано методическое описание постановки задач 'для построения разбиений Вороного и Делоне, вводятся критерии сравнения триангуляции, на основе критерия оптимальности Лаусона показывается геометрическая оптимальность двумерной триангуляции Делоне, приводятся необходимые определения для n-мерного случая.
В п. I.I.I описывается разбиение Вороного на плоскости, для конечного набора р из N дискретных точек как набора ячеек вида
v(i)={xR2| d(2,p1)«J(x,p;J), 3=1,' .... N}. Триангуляция Делоне вводится как разбиение, двойственное к разбиению Вороного и приводится результат о эквивалентности данного определения и критерия описанного круга.
В п. I.I.2 на осноеє индекса триангуляции еводится критерий сравнения триангуляции, а з п. I.1.3 с использованием процедуры Лаусона локальной оптимизации углов показано, что только триангуляция Делоне удовлетворяет этому критерию..
В п. І.І.4 приводятся обще определения n-мерных разбиений Делоне и приводятся классические результаты.
В 1.2 для двумерных областей строятся алгоритмы генерации множества дискретных точек Р. в области с учетом контроля
I
минимального расстояния между ними. Предлагаются 3 подхода к генерации точек равномерно в области: прямоугольная равноугольная и хаотическая структура и для рассматриваемых подходов приводятся оценки максимального расстояния между точками.
В 1.3 вводится понятие функции плотности точек и предлагаются подхода к генерации точек в области с учетом условия неравномерного.распределения точек на основе функций плотности со значениями в заданных точках. Для получения непрерывных функций плотности по значениям в заданных точках предлагаются три подхода: кусочно-линейное восполнение на основе вспомогательной триангуляции Делоне, радиальное восполнение и восполнение на основе^ квадро-деревьев.
В 1.4 на основе модификации алгоритма Ватсона предложен алгоритм построения двумерной триангуляции Делоне. При этом, в целях повышения эффективности, используются вспомогательные структуры, использованные в п.1.2 для контроля минимального расстояния между точками. Для предложенного алгоритма приводятся апостериорные оценки на основе оценки быстродействия порядка 'od/1). В этом случае параметр а определяется
экспериментально для конкретных задач как
log (Т/Т,,)
а = log (11/) ' и в результате экспериментов получено значение а<1.2.
В 1.5 вводится понятие динамических n-мерных триангуляции Делоне как последовательности триангуляции Делоне
DT± = M(P(h±(x)) с различными функциями плотности h±(i) для заданной области. Для таких триангуляции рассматривается задача построения VfE± по DT^.,. Для эффективного решения этой задачи доказывается
теорема об удалении точек из триангуляции Делоне с сохранением условия Делоне:
Теорема. Для произвольной точки q множества Р разбиение DT(P-q) = DT(P) - N(q) + DT(N(q)), построенное по множеству P-q является разбиением Делоне, где N(q) - набор многогранников с вершиной q, и ra(N(q)) - разбиение Делоне многогранника, образованного элементами из M(q).
Вторая глава посвящена построении алгоритмов дискретизации трехмерных областей, описанию алгоритма построения квазирегулярных сеток и описанию построения методов тестирования трехмерных сеток.
В 2.1 описываются основные требования и подходы к построению алгоритмов дискретизации трехмерных областей и описаны: недостатки разбиений Делоне в трехмерном случае.
В 2.2 предлагается алгоритм построения трехмерных сеток на основе локально-модифицированного подхода и рассматриваются соответствующие геометрические задачи. В и. 2.2.Г вводится определение логический функции области и приведены общие подходы к заданию таких функций. В п. 2.2.2 описан класс рассматриваемых поверхностей, в п. 2.2.3 решается задача сдвига околограничной точки на поверхность второго порядка, а в п. 2.2.4 описано решение задачи локализации точки в области на основе введенных описаний. В п. 2.2.5 введено общее определение сетки в трехмерном пространстве и доказана его эквивалентность классическим определениям. В п. 2.2.6- п. 2.2.8 описывается алгоритм построения лока^но-модифищпхдаанных сеток и приводится классификация пересечений ячеек с границей. В п. 2.2.9 вводится понятие триангуляции Делоне на поверхности.
которое используется для разбиения граничных элементов при построении сетки.
В 2.3 строятся подходы и алгоритмы численного тестирования трехмерных сеток на основе топологических характеристик границы и условий согласованности. Для предложенных подходов вводятся необходимые определения и доказываются утверждения о необходимости и достаточности предложенных подходов. Как результат - разработано II универсальных методов, называемых тестами, для автоматического определения корректности произвольных трехмерных сеток.
В п. 2.3.1 вводится общая постановка задачи и приведены необходимые определения. В п. 2.3.2 сформулированы два условия проверки и доказано, что их выполнение является достаточным для согласованности сетки. В п. 2.3.3 описано условие проверки топологии внутренних узлов сетки. В п. 2.3.4- п. 2.3.5 на основе теории выпуклых многогранников строятся необходимые и достаточные условия проверки топологии граничных компонент сетки. В п. 2.3.6 приведены условия тестирования глобальных топологических характеристик границы сетки - ориентируемости, замкнутости и топологического класса.
Третья глава посвящена задачам автоматизации процесса моделирования двумерных и трехмерных электростатических полей методом конечных элементов в системе Stella-2D и ППП РАМЗЕС-3.
В 3.1 описаны постановка задачи и особенности реализации алгоритмов в указанных программных системах. Особое внимание уделяется описанию поэлементной реализации итерационных методов при решении систем линейных уравнений, возникающих в МНЭ. В качестве предобуславливаицих матриц предлагается
- II -
использовать различные варианты матриц методов расщеплений (переменных направлений).
В 3.2 описывается интерактивная система Stella-2D для расчета электростатических полей в плоской и осесимметричной постановках на персональных компьютерах. Система stella-2D построена на основе - графической системы меню и включает графический редактор, базу данных материалов, генератор сеток, расчетный блок и блок графического анализа решения. В качестве примеров использования системы stelia-2D приводятся результаты практических расчетов изолирующей подвески ЛЭП и полимерного опорного изолятора, имеющих сложную геометрию области с соотношением размеров порядка I02- I03.
В 3.3 описывается пакет прикладных программ РАМЗЕС-3 для решения трехмерных электростатических задач. Приводятся общее описание пакета, схема вычислительного процесса, структура входного языка АВДА для описания геометрии, описание систем возникающих уравнений. Приводятся примеры описания трехмерной ' области на входном языке, примеры сеток для различных областей и пример использования пакета для решения стационарной задачи теплопроводности.
В заключении формулируются основные результаты работы.
В приложениях приводятся таблица разработанных алгоритмов для двумерных и трехмерных задач построения сеток, пример квадродерева, варианты исходных сеток, псевдо-код и иллюстрация к классификации пересечений для трехмерного алгоритма, примеры двумерных и трехмерных сеток, результаты решения практических задач.
