Аппроксимация решения кусочно-квадратичными функциями

Для решения той же задачи ( 3 .0) – ( 3 .1) воспользуемся квадратичной аппроксимацией в пределах одного конечного элемента с центральной точкой x_k. Как и в предыдущем случае, решение раскладывается по пробным функциям

, (3.12)

имеющим вид

.

Невязка уравнения ( 3 .0), получаемая на решении ( 3 .12), взвешивается с использованием тех же функций _i, _j и _k,

(3.13)

Преобразуем первое из этих уравнений:

,

.

.

Учитывая, что , и используя разложение ( 3 .12), приходим к выражению

.

Выполняя аналогичные преобразования с оставшимися выражениями в ( 3 .13), приходим к системе уравнений

Подсчитаем значения интегралов в полученных выражениях.

;

,

,

,

,

,

,

,

,

.

Подстановка найденных значений приводит к системе уравнений

Эта же система в матричной форме принимает вид

. (3.14)

Суммируя все уравнения этой системы получаем

,

уже известное условие теплового баланса ( 3 .11).

Пример 3.2. Рассмотрим задачу из примера 3.1 с теми же исходными данными. Пусть весь стержень аппроксимируется одним конечным элементом. Будем считать, как и в предыдущем случае, что на его левом конце задана температура , а на правом – граничные условия третьего рода

.

Для рассматриваемого случая система уравнений приводится к виду

,

Для принятых L, W, , и  эта система уравнений принимает вид

и имеет решение T_i=100 (левый конец стержня), T_j = 153/2 (правый конец), T_k = 619/7 (центр стержня). С учетом вида пробных функций ( 3 .2) решение запишется в виде

.

Полученное выражение является точным решением этой задачи.