Решение задач методом конечных разностей

Метод конечных элементов (МКЭ) — численный метод решения дифференциальных уравнений с частными производными, а также интегральных уравнений, возникающих при решении задач прикладной физики. Метод широко используется для решения задач механики деформируемого твёрдого тела (сопромата), теплообмена, гидродинамики и электродинамики.

, (1)

 

При начальном условии (2)

При краевых условиях

U(S,t) = (3)

Из уравнений (3) мы видим, что данная область 0 ≤S; 0≤t≤T заменим данную область сеткой, в каждом внутреннем узле вычислим значение функции.

В каждом внутреннем узле сетки значение функции будет определять из конечно-разностного уравнения. Для этого диф-е уравнение аппроксимирует следующим образом:

Примем

(4)

Уравнение (4) называется явной схемой. Данная система является устойчивой при

Рассмотрим вариант не явной схемы:

3. Решение краевой задачи методом Фурье:

Метод разделения переменных или метод Фурье, является одним из наиболее распространенных методов решения уравнений с частными производными. Изложение этого метода мы проведем для задачи о колебаниях струны, закрепленной на концах.

Итак, будем искать решение уравнения

,

(1)

удовлетворяющее однородным граничным условиям

U (0, t) = U (l, t) = 0

(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) должна удовлетворять дополнительным условиям

X (0) = X (l) =0,

(9)

так как иначе мы имели бы T (t)≡0 и U (x, t)≡0, в то время как задача состоит в нахождении нетривиального решения.

Таким образом, в связи с нахождением функции X (x) мы приходим к простейшей задаче о собственных значениях: найти такие значения параметра λ, при которых существуют нетривиальные решения задачи:

(10)

а также найти эти решения. Такие значения параметра называются собственными значениями, а соответствующие им нетривиальные решения – собственными функциями задачи (10).

Итак, найдем знак :

1 случай , например, .
Запишем характеристическое уравнение для (10):

.

Общее решение уравнения может быть записано в виде
.

Граничные условия дают:
,
то есть и .
Но в рассмотренном случае - действительно и положительно, так что .
Поэтому , и, следовательно, , а мы ищем нетривиальное решение.

2 случай Пусть .

При также не существует нетривиальных решений. Действительно, в этом случае общее решение уравнения (7) имеет вид
.
Граничные условия дают:

то есть A =0 и B =0 и, следовательно, .

3 случай , например .

Характеристическое уравнение имеет вид

.

Общее решение уравнения:
.

Граничные условия дают:
.
Если , то . Поэтому

, где n - любое целое число. Обозначим p через ,

.

- нетривиальное решение задачи (10),

(11)

определяемое с точностью до произвольного множителя, который мы положили равным единице. Этим же значениям соответствуют решения уравнения (8).

,

(12)

где и - произвольные постоянные.

Возвращаясь к задаче (1) – (3), заключаем, что функции
являются частными решениями уравнения (1), удовлетворяющими граничным условиям (3) и представимыми в виде произведения (4) двух функций.

Обратимся к решению в общем случае. В силу линейности и однородности уравнения (1) сумма частных решений

(13)

также удовлетворяет этому уравнению и граничным условиям (2).

Начальные условия позволяют определить и . Потребуем, чтобы функция (13) удовлетворяла условиям (3):

.

(14)

Если функции и удовлетворяют условиям разложения в ряд Фурье, то

.

(15)

Подставив (15) в (13), мы удовлетворим краевым условиям и получим решение уравнения.