Метод конечных элементов (МКЭ) — численный метод решения дифференциальных уравнений с частными производными, а также интегральных уравнений, возникающих при решении задач прикладной физики. Метод широко используется для решения задач механики деформируемого твёрдого тела (сопромата), теплообмена, гидродинамики и электродинамики.
, (1)
При начальном условии (2)
При краевых условиях
U(S,t) = (3)
Из уравнений (3) мы видим, что данная область 0 ≤S; 0≤t≤T
заменим данную область сеткой, в каждом внутреннем узле вычислим значение функции.
В каждом внутреннем узле сетки значение функции будет определять из конечно-разностного уравнения. Для этого диф-е уравнение аппроксимирует следующим образом:
Примем
(4)
Уравнение (4) называется явной схемой. Данная система является устойчивой при
Рассмотрим вариант не явной схемы:
3. Решение краевой задачи методом Фурье:
Метод разделения переменных или метод Фурье, является одним из наиболее распространенных методов решения уравнений с частными производными. Изложение этого метода мы проведем для задачи о колебаниях струны, закрепленной на концах.
Итак, будем искать решение уравнения
![]() ![]() | (1) |
удовлетворяющее однородным граничным условиям
![]() | (2) |
и начальным условиям
![]() ![]() | (3) |
Уравнение (1) линейно и однородно, поэтому сумма частных решений также является решением этого уравнения. Будем искать решение уравнения в виде
![]() ![]() | (4) |
где X (x)- функция только переменного ,
T (t)- функция только переменного
.
Подставим (4) в уравнение (1), получим:
![]() ![]() | (5) |
Чтобы функция (4) была решением уравнения (1), равенство (5) должно удовлетворяться тождественно, то есть для всех значений независимых переменных ,
. Правая часть равенства (5) является функцией только переменного x, а левая- только
.
Фиксируя, например, некоторое значение х и меняя t (или наоборот), получим, что правая и левая части (5) при изменении своих аргументов сохраняют постоянное значение, то есть
![]() ![]() | (6) |
Из соотношения (6) получаем обыкновенные дифференциальные уравнения для определения функций X (x) и T (t).
![]() ![]() | (7) (8) |
Граничные условия (2) дают:
.
Отсюда следует, что функция X (x) должна удовлетворять дополнительным условиям
![]() | (9) |
так как иначе мы имели бы T (t)≡0 и U (x, t)≡0, в то время как задача состоит в нахождении нетривиального решения.
Таким образом, в связи с нахождением функции X (x) мы приходим к простейшей задаче о собственных значениях: найти такие значения параметра λ, при которых существуют нетривиальные решения задачи:
![]() ![]() | (10) |
а также найти эти решения. Такие значения параметра называются собственными значениями, а соответствующие им нетривиальные решения – собственными функциями задачи (10).
Итак, найдем знак :
1 случай , например,
.
Запишем характеристическое уравнение для (10):
.
Общее решение уравнения может быть записано в виде
.
Граничные условия дают:
,
то есть и
.
Но в рассмотренном случае - действительно и положительно, так что
.
Поэтому ,
и, следовательно,
, а мы ищем нетривиальное решение.
2 случай Пусть .
При также не существует нетривиальных решений. Действительно, в этом случае общее решение уравнения (7) имеет вид
.
Граничные условия дают:
то есть A =0 и B =0 и, следовательно, .
3 случай , например
.
Характеристическое уравнение имеет вид
.
Общее решение уравнения:
.
Граничные условия дают:
.
Если , то
. Поэтому
, где n - любое целое число. Обозначим p через
,
.
![]() ![]() | (11) |
определяемое с точностью до произвольного множителя, который мы положили равным единице. Этим же значениям соответствуют решения уравнения (8).
![]() ![]() | (12) |
где и
- произвольные постоянные.
Возвращаясь к задаче (1) – (3), заключаем, что функции
являются частными решениями уравнения (1), удовлетворяющими граничным условиям (3) и представимыми в виде произведения (4) двух функций.
Обратимся к решению в общем случае. В силу линейности и однородности уравнения (1) сумма частных решений
![]() ![]() | (13) |
также удовлетворяет этому уравнению и граничным условиям (2).
Начальные условия позволяют определить и
. Потребуем, чтобы функция (13) удовлетворяла условиям (3):
![]() ![]() | (14) |
Если функции и
удовлетворяют условиям разложения в ряд Фурье, то
![]() ![]() | (15) |
Подставив (15) в (13), мы удовлетворим краевым условиям и получим решение уравнения.