LEC-22 (Материалы к лекциям), страница 2
Описание файла
Файл "LEC-22" внутри архива находится в следующих папках: Материалы к лекциям, Lecturessemestr7. Документ из архива "Материалы к лекциям", который расположен в категории "". Всё это находится в предмете "методы решения задач механики сплошных сред" из 7 семестр, которые можно найти в файловом архиве МАИ. Не смотря на прямую связь этого архива с МАИ, его также можно найти и в других разделах. Архив можно найти в разделе "лекции и семинары", в предмете "методы решения задач механики сплошных сред" в общих файлах.
Онлайн просмотр документа "LEC-22"
Текст 2 страницы из документа "LEC-22"
3. Формирование треугольных элементов посредством соединения сгенерированных узлов пересекающимися ребрами.
4. Оптимизация сгенерированной сетки.
В некоторых программах предварительное нанесение узлов сетки выполняется вручную. Другие программы, имеющие более высокую степень автоматизации, содержат модули генерации узлов согласно заданной плотности их распределения в пределах каждой подобласти и на границе всей области.
Для объединения узловых точек в треугольные элементы используют различные модификации стратегии образования регулярных СКЭ.
Одна из возможных стратегий – алгоритмы Делоне, рассматривается в дальнейшем. Оптимизацию сгенерированных сеток также выполняют различными способами, но чаще других применяется геометрическая регуляризация. Примерами является опять-таки алгоритмы Делоне и алгоритмы Рапперта, обсуждаемые в дальнейшем.
Основным достоинством автоматизированных вариантов данного метода является непосредственное использование заданной плотности распределения узлов в процессе построения СКЭ. Однако метод требует обязательного контроля взаимного расположения сгенерированных узловых точек сетки, что намного усложняется и замедляет работу алгоритма.
9.2.2.5. Метод выравнивание - выемка
Данный метод является одним из наиболее эффективных методов формирования сетки. При использовании метода реализуются большие возможности контроля и управления построением при небольшом числе исходных данных.
Для двухмерных областей произвольного вида наиболее перспективным является треугольный элемент, поскольку треугольная сетка удовлетворительно аппроксимирует поверхности любой кривизны, не накладывает ограничений на размеры элементов и позволяет легко осуществить их контроль. Это предопределило выбор треугольного элемента в качестве базового при создании метода выравнивание-выемка.
Рассмотрим основные этапы этого метода.
Плоская многосвязная область, подлежащая триангуляции, представляется в виде непересекающегося объединения односвязных базовых подобластей, граница которых состоит из непересекающегося объединения базовых линий.
Две (или более) базовых линий могут иметь общие концевые точки (узлы). Таким образом, простейшим элементом топологической модели плоского континуума является точка, которая определяется порядковым номером и координатами. Другим элементом является базовая линия, которая также имеет свой порядковый номер с указанием номеров, образующих ее базовых узлов. Базовая подобласть, в свою очередь имеет номер с указанием номеров ограничивающих базовых линий.
Описанная структура характеризует схему соединения областей.
Поскольку линию можно задавать несколькими определенным образом согласованными точками, целесообразно разработать шаблон или логические схемы преобразования множества исходных данных.
Рассмотрим базовые элементы трех типов: прямую, дугу окружности и квадратичную параболу.
П
рямолинейный элемент может быть образован, если известно положение начального и конечного базовых узлов с расположенными между ними промежуточными точками. Для получения неравностоящих промежуточных узлов можно задать их число и функцию, характеризующую вид разбиения линии.
Базовый элемент в форме дуги окружности определяют заданием трех базовых узлов, причем средний узел не обязательно должен совпадать с каким-либо промежуточным узлом, - он предназначен для определения координат центра окружности. Дискретизация базового кругового элемента выполняется разбивкой соответствующего ему центрального узла.
Базовый параболический линейный элемент, подобно базовому элементу в форме дуги окружности, задается тремя базовыми узлами. Очевидно, что вид параболы зависит от выбора координатных осей. Для однозначного определения процедуры формирования параболического элемента первый и последний базовые узлы соединяют прямой, которая образует локальную ось OX ’.
Средняя точка 0, соединяющей линии принимается за начало локальной оси y’ . Базовый параболический элемент в координатной системе x’Oy’ строится единственным образом и имеет симметричную форму. Разделение базового параболического элемента на несколько отрезков выполняют разбивкой соединяющей линии и проведением перпендикуляров к оси Оx’, которые пересекают параболу в запроектированных точках. Неравномерную разбивку производят аналогично. Если необходимо точно определить положение промежуточных узлов, используют численное интегрирование.
Топологической модели плоского континуума соответствует непересекающееся объединение односвязных подобластей в форме многоугольников, аппроксимирующих базовые подобласти.
На основании информации о сформулированных граничных узлах (узлах дискретизации линий) производится автоматическая триангуляция базовых подобластей. Сетка КЭ строится внутрь подобласти от базовых линий с учетом локальных свойств текущей границы.
Алгоритм метода характеризуется двумя способами формирования треугольных элементов.
1) выравнивание, т.е. уменьшение текущей границы;
2
) выемкой, т.е. построение нового узла текущей границы.
Применение конкретного способа образования КЭ определяется проверкой локальных свойств текущей границы в соответствии с рядом установленных критериев.
Основным моментом процесса триангуляции для одной из модификаций алгоритма является поиск минимального угла , образованного двумя смежными отрезками прямых на текущей границе. Если 1 75о , то происходит формирование треугольного элемента способом выравнивания.
Если 75о 89о , то выполняется проверка двух других углов 2 , 3 треугольного элемента, построенного способом выравнивания.
В случае, когда эти углы отвечают требованиям регулярности сетки 2 30о и 3 30о, построение считается законченным.
При 1 > 89о или 2 < 30о или 3 < 30о происходит формирование двух треугольных элементов способом выемки.
Угол 1 делится на две равные части, а длина стороны ВД выбирается равной среднему арифметическому длин сторон АВ и ВС, умноженному на числовой коэффициент, который характеризует необходимую степень сгущения или разрежения сетки внутри базовой области.
В подобной ситуации можно использовать способ выравнивания. Тогда процесс триангуляции будет сопровождаться формированием отдельных элементов, не удовлетворяющих перечисленным ранее критериям, но такое построение не приведет к вырождению сетки.
В процессе работы алгоритма триангуляции количество построенных треугольных элементов (внутренних узлов) определяется:
1) количеством граничных узлов на базовых линиях;
2) законами распределения узлов на базовых линиях;
3) суммарной площадью базовых подобластей;
4) выбором определяющих критериев формирования треугольных элементов.
Вместе с тем, очевидно, что построение оптимальной сетки без учета функции плотности элементов невозможно.
