- •А.В. Румянцев
- •Содержание
- •Глава 6. Программная реализация метода
- •Глава 1. Краевые задачи теории поля
- •1.1 Уравнение переноса в обобщенной криволинейной системе координат
- •1.2 Краевые условия задачи
- •1.3* Краткая характеристика методов решения краевой задачи
- •Глава 2. Метод конечных элементов в краевых
- •2.1 Методы взвешенных невязок
- •2.2 Основная концепция метода конечных элементов
- •Глава 3. Геометрические аспекты мкэ
- •3.1 Типы конечных элементов. Базовый каталог элементов
- •3.2 Дискретизация области на элементы
- •А) с разными (вид сверху); б) с одинаковыми.
- •Цифры – это номер элемента по каталогу
- •3.3 Нумерация элементов и узлов
- •3.4 Индексация узлов и формирование таблицы входных данных
- •Осесимметричной детали
- •Геометрическая часть таблицы входных данных
- •Глава 4. Математическое описание элемента
- •4.1 Метод Крамера
- •4.2 Метод Лагранжа
- •4.3 Обобщенный метод Крамера-Лагранжа
- •4.4 Эрмитовы элементы
- •4.5 Свойства базисных функций элемента
- •Глава 5. Вычислительные аспекты мкэ
- •5.2 Матричное представление элементного вклада
- •Производные базисных функций
- •5.3 Формирование глобальных матриц для исследуемой области
- •А) Сокращенная
- •5.4 Стандартизация матриц элементов
- •5.5 Естественная система координат
- •5.6 Средние температуры элемента
- •Глава 6. Программная реализация мкэ
- •6.1 Задание краевых условий задачи
- •6.2 Решение системы динамических уравнений
- •Временная циклограмма q(τ)
- •6.З Учет температурной зависимости теплофизических параметров
- •6.4 Радиационный компонент теплообмена
- •6.5 Сходимость, полнота и согласованность, точность
- •Базовый каталог объемных элементов
- •Осесимметричные объемные элементы
- •Базовый каталог одно- и двумерных элементов
3.4 Индексация узлов и формирование таблицы входных данных
При математическом описании элемента необходимо узлы элемента индивидуализировать, т.е. сделать их отличимыми друг от друга.
С этой целью каждому узлу присваивается индекс в зависимости от количества узлов элемента. Выбор-го узла элементапроизволен, индексация остальных узлов выполняется последовательно в направлении против часовой стрелки. Проиндексировав узлы, составляют таблицу, которая ставит в соответствие индексы элемента глобальным номерам его узлов. С помощью этой таблицы впоследствии осуществляется включение матрицы элемента в соответствующие глобальные.
Фиксация положения элемента в пространстве осуществляется заданием координат его проиндексованных и занумерованных узлов.
Процесс дискретизации области определения задачи завершается формированием геометрической части таблицы входных данных.
На рис. 3.5 изображена двумерная (для наглядности) область,
Рис. 3.5 Пример дискретизации сечения
Осесимметричной детали
разбитая на два элемента. Узлы, выбранные в качестве i-х, обозначены звездочками. В скобках указаны номера элементов и координаты узлов . Геометрическая часть таблицы входных данных с указанием номера элемента по каталогу имеет следующий вид.
Таблица входных данных должна содержать, естественно, сведения о свойствах материала элемента – его механических, теплофизических, электрических и т.п. свойствах, соответствующих физической природе задачи. Эти сведения целесообразно объединить в физический каталог, в котором каждому материалу присвоен номер . Например, девятый номер присвоен ниобию, а третий - нержавеющей стали. Введение в программу постоянно пополняемого физического каталога унифицирует ее и делает компактной таблицу входных данных, поскольку громоздкие, занимающие несколько разрядов данные о свойствах материала элемента (притом для каждого, даже если они одинаковые) заменяются не более чем тремя разрядами номера материала по физическому каталогу.
Задание 3
3.1 Дискретизируйте объект – треугольную прямую призму – на два неодинаковые по размерам элемента. Задайте координаты узлов; составьте таблицу входных данных.
3.2 Дискретизируйте делением по углу и по длине на 4 элемента объект,
Таблица 2
Геометрическая часть таблицы входных данных
Номер элемента |
Номер элемента по физическому каталогу |
Номер элемента по базовому каталогу |
Индекс и номер узла |
Координаты узлов | |
(e) |
eфк |
ekei |
i j k l |
Ri Rj Rk Rl |
Zi Zj Zk Zl |
(1) |
9 |
42 |
1 3 4 2 |
0,0 2,0 2,0 0,0 |
1,0 1,0 3,0 3,0 |
(2) |
3 |
52 |
3 5 6 4 |
2,0 3,0 3,0 2,0 |
1,0 1,0 3,0 3,0 |
представляющий собой четверть (по углу) полого цилиндра. Задайте координаты узлов; составьте таблицу входных данных.
3.3 Какие свойства материала элемента следует внести в физический каталог в задачах теплопроводности?