Лекция 4. Расширенные и глобальные вектора и матрицы. Сборка системы МКЭ
Вектора и матрицы в интегральном тождестве (3.10) определены на каждом конечном элементе и имеют размерности
и
по числу узлов в элементе. Для произвольного элементного вектора
определим расширенный элементный вектор
длиной
(по числу узлов МКЭ-сетки), у которого отличны от нуля только
компонент
. Например, для элемента 28 сетки, показанной на рис. 3.2,
(это записано в строке 28) таблицы связности. Поэтому компоненты элементного вектора
образуют ненулевые компоненты расширенного элементного вектора
Аналогично для произвольной элементной матрицы размерности
определим расширенную элементную матрицу
размерностью
, у которой отличны от нуля только
компонентов, равных
и расположенных в позициях
. Для нашего треугольника
это элементы
По существу, расширенные элементные вектора и матрицы получаются из элементных векторов и матриц рассылкой их элементов на нужные позиции в соответствии с глобальной нумерацией узлов элемента, определенной таблицей связности.
Замечательное свойство расширенных векторов и матриц определяется расположением в них нулей и выражается следующими равенствами
.
Используя эти свойства, можно переписать интегральное тождество (3.10) а терминах расширенных векторов и матриц:
(4.1)
Определим, далее, глобальные вектора и матрицы как сумму расширенных элементных векторов и матриц:
.
В этих терминах можно переписать равенство (4.1) в виде
.
Поскольку глобальный вектор в последнем равенстве может выбираться произвольно, то это значит, что выражение в фигурных скобках должно равняться нулю. Так мы приходим к глобальной системе уравнений МКЭ
(4.2)
или
(4.3)
Матрица этой системы складывается из глобальных матриц масс
, жесткости
и матрицы
, порожденной граничными условиями 3-го рода. Каждая из этих матриц, в свою очередь, представляет собой сумму вкладов от соответствующих элементных матриц. Структура правой части
уравнения (4.3) аналогична в том смысле, что также представляет собой сумму вкладов от всех элементных векторов. Такая структура системы уравнений МКЭ определяет алгоритм поэлементной сборки ее матрицы и правой части. Вначале все элементы глобальной матрицы
и глобального вектора правой части
обнуляются, затем в цикле по конечным элементам
вычисляются элементные матрицы и вектора и добавляются к соответствующим элементам
и
. Это соответствие определено в таблице 2, где для каждого элемента
указаны глобальные номера его узлов. Непосредственно из алгоритма сборки видно, что матрица
будет разреженной: в ее
-ую строку, соответствующую
-му узлу сетки, очевидно попадут вклады от элементов, инцидентных данному узлу (список таких элементов содержится в таблице 4). При этом номера столбцов
-ой строки, куда попадут эти вклады, совпадает со множеством глобальных номеров узлов инцидентных
-му узлу сетки элементов. Эти номера и перечислены в
-ой строке таблицы 3. Аналогично производится поэлементная сборка вектора
правой части системы уравнений МКЭ.
Итак, чтобы собрать систему уравнений МКЭ нужно
1. Триангулировать область , т.е. построить МКЭ-сетку.
2. Определить базисные функции на каждом конечном элементе.
3. Вычислить элементные вклады с помощью интегралов (3.4), (3.5), (3.6), (3.7), (3.9).
4. Выполнить поэлементную сборку глобальных векторов и матриц.
Проблемы 1 и 4 нами решены, а задачи 2 и 3 решаются с помощью т.н. базисного элемента, на котором легко и однообразно вводятся базисные функции и выполняются операции дифференцирования и интегрирования.