Линейные базисные функции базисного треугольного элемента.

Линейные базисные функции базисного треугольного элемента. - student2.ru (4.11)

Вторая строчка в определении (4.11) может быть использована в качестве признака принадлежности произвольной точки Линейные базисные функции базисного треугольного элемента. - student2.ru базисному треугольнику Линейные базисные функции базисного треугольного элемента. - student2.ru . Произвольная функция Линейные базисные функции базисного треугольного элемента. - student2.ru со значениями Линейные базисные функции базисного треугольного элемента. - student2.ru может быть представлена на элементе Линейные базисные функции базисного треугольного элемента. - student2.ru разложением по базису:

Линейные базисные функции базисного треугольного элемента. - student2.ru . (4.12)

Градиенты линейных базисных функций (4.11) – это векторы

Линейные базисные функции базисного треугольного элемента. - student2.ru (4.13)

Обратим внимание, что градиенты линейных базисных функций не зависят от координат.

Коэффициенты элементных матриц выражаются через интегралы от базисных функций и их произведений. Поэтому вычислим следующие интегралы

Линейные базисные функции базисного треугольного элемента. - student2.ru ,

а также интегралы от базисных функций и их произведений по границе.

Линейные базисные функции базисного треугольного элемента. - student2.ru

Аналогично получаем

Линейные базисные функции базисного треугольного элемента. - student2.ru .

Интегралы от произведений.

Линейные базисные функции базисного треугольного элемента. - student2.ru

Аналогично вычисляются интегралы

Линейные базисные функции базисного треугольного элемента. - student2.ru .

Оставшиеся интегралы от произведения разных базисных функций равны

Линейные базисные функции базисного треугольного элемента. - student2.ru

Аналогично

Линейные базисные функции базисного треугольного элемента. - student2.ru .

Итак,

Линейные базисные функции базисного треугольного элемента. - student2.ru (4.14)

Теперь, используя разложение (4.12) произвольной функции Линейные базисные функции базисного треугольного элемента. - student2.ru по базису Линейные базисные функции базисного треугольного элемента. - student2.ru , легко подсчитать интеграл (4.9) этой функции по области:

Линейные базисные функции базисного треугольного элемента. - student2.ru (4.15)

Вторая формула (4.14) определяет матрицу масс (3.4):

Линейные базисные функции базисного треугольного элемента. - student2.ru (4.16)

Используя формулу (4.7), разложение (4.12) и градиенты базисных функций (4.13), выведем формулу дифференцирования произвольной функции Линейные базисные функции базисного треугольного элемента. - student2.ru .

Линейные базисные функции базисного треугольного элемента. - student2.ru (4.17)

Вспомним, что согласно (3.5) элементы матрицы жесткости определялись как

Линейные базисные функции базисного треугольного элемента. - student2.ru

Вычислим вначале

Линейные базисные функции базисного треугольного элемента. - student2.ru

Теперь нетрудно подсчитать элементы матрицы жесткости, попарно перемножая скалярно выписанные вектора:

Линейные базисные функции базисного треугольного элемента. - student2.ru (4.18)

Итак, все элементные матрицы и векторы, необходимые для сборки системы МКЭ, определены. Напомним, что в линейной алгебраической системе уравнений (4.3)

Линейные базисные функции базисного треугольного элемента. - student2.ru (4.19)

матрица масс Линейные базисные функции базисного треугольного элемента. - student2.ru собирается из элементных вкладов (4.16), матрица жесткости Линейные базисные функции базисного треугольного элемента. - student2.ru состоит из элементов (4.18), компоненты матрицы Линейные базисные функции базисного треугольного элемента. - student2.ru определены в (3.7) интегралами по границе, вектора Линейные базисные функции базисного треугольного элемента. - student2.ru и Линейные базисные функции базисного треугольного элемента. - student2.ru состоят из компонент (3.8) и (3.9). Вспомним, как это выглядело для треугольника, у которого узлы 1 и 3 граничные:

Линейные базисные функции базисного треугольного элемента. - student2.ru (4.20)

Линейные базисные функции базисного треугольного элемента. - student2.ru (4.21)

Лекция 5. Программирование МКЭ

Структура данных

Нам необходимы глобальная матрица Линейные базисные функции базисного треугольного элемента. - student2.ru и глобальный вектор Линейные базисные функции базисного треугольного элемента. - student2.ru системы МКЭ, а также вектор неизвестных узловых значений Линейные базисные функции базисного треугольного элемента. - student2.ru . Здесь Линейные базисные функции базисного треугольного элемента. - student2.ru ‑ число узлов МКЭ сетки. Сразу отметим, что матрица Линейные базисные функции базисного треугольного элемента. - student2.ru разрежена, поэтому на практике хранят толко ее ненулевые элементы. Вопросы компактного хранения и алгебры разреженных матриц будут рассмотрены в следующем курсе, посвященном решению сеточных уравнений. Пока формально, для ясности понимания, будем использовать заполненные матрицы.

Type

float=extended;

TGlobMat=array[1..N;1..N] of float;

TGlobVec= array[1..N] of float;

sGlobMat=^TGlobMat;

sGlobVec=^TGlobVec;

Var

B:sGlobMat;

U,F: sGlobVec;

Сама сетка, как уже отмечалось, хранится в двух таблицах.

Type

TTable1=record

x,y:float;

b:integer;

end;

TTable2=record

n1,n2,n3:integer;

end;

arTab2=array[1..M] of TTable2;

arTab1= array[1..N] of TTable1;

sTab1=^arTab1;

sTab2=^arTab2;

Var

Tab1:sTab1;

Tab2:sTab2;

Наши рекомендации