- •Выполнил:
- •Алгоритм решения
- •Реализация алгоритма
- •Z ai |ф*ФА = /ф* cosxdx’ к = 0,4,
- •Разрешающие соотношения
- •Алгоритм решения
- •Результаты расчётов
- •Разрешающие соотношения
- •Алгоритм решения
- •Реализация алгоритма
- •2.1. Аппроксимация решения кусочно-линейными функциями
- •Алгоритм решения
- •Результаты расчетов
- •Разрешающие соотношения
- •Алгоритм решения
- •Реализация алгоритма
- •Разрешающие соотношения
- •Алгоритм решения
- •Реализация алгоритма
- •3.1. Явная разностная схема
- •Алгоритм решения
- •Реализация алгоритма
- •Разрешающие соотношения
- •Программа 3.2
- •Реализация алгоритма
- •Алгоритм решения
- •Реализация алгоритма
Исключённые из системы уравнения в дальнейшем, после вы числения узловых коэффициентов 7}, могут быть использованы для определения тепловых потоков 0 ° и
Алгоритм решения
Ниже приведён фрагмент программы на языке Си, реализующей процедуру построения приближённого решения дифференциального уравнения с использованием кусочно-линейных пробных функций.
Программа 2.2
// |
Стационарная задача теплопроводности |
|
|
|||||
// |
Ne |
число сегментов |
|
|
|
|||
// Nk - число узлов |
|
|
|
|||||
// |
locale |
массив |
коэффициентов |
системы уравнении |
для |
сегмента |
||
// |
localF |
массив |
правых частей |
системы уравнений |
для |
сегмента |
||
// |
globalC |
массив |
коэффициентов |
системы уравнений |
для |
объекта |
||
// |
globalF |
массив |
коэффициентов |
системы уравнений |
для |
объекта |
||
// |
х |
массив координат узлов |
|
|
|
|||
// |
хО |
начальная |
точка |
|
|
|
||
/ / |
XI |
конечная |
точка |
|
|
|
//h длина сегмента
//W0 амплитуда мощности источников
//L коэффициент теплопроводности
// |
ТО |
температура |
на левом торце (хО) |
// |
Те |
температура |
среды |
/ / А |
коэффициент |
теплоотдачи |
#define Ne 64
void main(void)
{double globalC[2*Ne+l] [ 2*Ne+l], globalF[2*Ne+l], x[Ne+l]; double locale[2] [2], localF[2], x0=0., xl=M__PI, h=(xl-xO)/Ne; double L = 7 0 .0, W=1000.0, T0=100.0, Te=20.0/ A=30.0;
int i, j, k, Nk=Ne+l;
//определение координат и номеров узлов каждого сегмента for(k=0; k<Nk; k++) x[k]=xO+h*k;
//инициализация массивов коэффициентов и правых частей
//системы уравнений
for(i=0; i<2*Ne+l; i++)
for(globalF[i]=0.0/ j=0; j<2*Ne+l; j++) globalC[i] [j]=0.0; for(k=0; k<2*Ne+l; k++)
// формирование матриц коэффициентов и правых частей для сегмента
{locale [0] [0]=7.0 *L/(3.0*h);
locale [0] [1] =localC[1] [0]= - 8 .0*L/(3.0*h) ; locale [0] [2] =localC [2] [0] =L/ (3.0*h) ; locale [1] [1]=16.0 * L / (3.0*h);
locale [1] [2]=localC [2] [1]= - 8 .0*L/(3.0*h) ;
locale[2] [2]=7.0*L/(3.0*h);
localF [0] =W0* ((3.0*sin (x [kj )+sin (x [k+1] ))/h- (4.0/ (h*h) -
1.0)*cos(x[k])+4.Q*cos(x [k+l])/ (h*h) );____________
l o c a l F [1]= W 0 * (8.0*(cos(x[k])-cos(x[k+1]))/ (h*h)- 4 . 0 * (sin(x[k])+sin(x[k+1]))/ h ) ;
localF [2]= W 0 * ((sin(x[k])+3.0*sin(x [k+1]))/ h + (4.0 / (h*h)- 1.0)*cos(x [k+1])-4.0*cos(x[k])/ (h*h));
// формирование матриц коэффициентов и правых частей для объекта globalC[2*k] [2*k]+=1оса1С [0] [0];
globalC [2*k] [2*k+l]+=localC [0] [1]; globalC[2*k] [2*k+2]+=localC [0] [2]; globalC[2*k+l][2*k]+=localC[1][0]; globalC[2*k+l][2*k+l]+=localC[1][1]; globalC[2*k+l][2*k+2]+=localC[1][2]; globalC[2*k+2] [2*k]+=localC [2] [0]; globalC [2*k+2] [2*k+l]+=localC [2] [1]; globalC[2*k+2] [2*k+2]+=localC [2] [2]; globalF[2*k]+=localF[0]; globalF[2*k+l]+=localF [1]; globalF[2+k+2]+=localF [2];
}
// определение граничного условия 1-го рода на левом торце globalC [0] [0] =1.0;
for(k=l; k<2*Ne+l; k++) g l o balC[0] [k]=0.0; globalF [0]=T0;
// определение граничного условия 3-го рода на правом торце globalC [2*Ne] [2*Ne]+=А;
globalF[2*Ne]+=A*Te;
//решение системы линейных алгебраических уравнений GAUSS(globalC, globalF, 2*Ne+l);
i _______________________________________________________________
Реализация алгоритма
Первоначально отрезок [0, п\ разбивается на четыре сегмента равной длины
[О, л] = [0,л/4] U[л/4,л/2] U [я/ 2,3л/4] U [3л/4,л].
На каждом из этих сегментов в соответствии с выражением (2.11) определяются кусочно-квадратичные пробные функции. Ре шение дифференциального уравнения представляется разложением (2.9). Для построения приближённого решения необходимо решить систему линейных алгебраических уравнений относительно коэффи циентов 7).
В соответствии с приведёнными выражениями и количествен ными значениями физических констант подсчитываются значения интегралов (с учетом определения (2.19) пробных функций):
- коэффициенты и правая часть системы уравнений для сег мента [0, л/4]:
к/Л —^
Mp'oWP <& = ~ = 207,962,
оЗА
*/} •/* о т
р Ж (*Ж/2(*)* = JWifi (*Ж (*)<Ь= - — = -237,671,
ПП
№ /2 (*)f л |
= ^ |
= 475,343, |
|
||
Я/4 |
Л/4 |
|
|
|
|
р-Фо(*Ж(*)^ = /*Ж (*Ж (*)* = — = 29,7089, |
|||||
0 |
0 |
|
|
Зп |
|
|
{x)fdx=1^ |
=207,962, |
|
||
л/4 |
я/4 |
|
|
|
|
р Ж /2(*Ж (*)<& = |
р Ж (jcK/2{x)dx = ~ 3h =-237,671; |
||||
|
sin h |
|
4 cos h |
||
|
~1Г |
|
|
1,03391, |
|
|
|
|
|
||
|
8(l-cos/t) |
4sinA |
= 197,299, |
||
|
|
h1 |
|
h |
|
|
3sinA |
4 |
.1 |
2, |
4 |
|
h |
—:—1 |
COSA--- r- = 94,5598; |
||
|
+ |
|
|
|
' 207,962Го - 237,6717;/2 + 29,70897; = 1,03391,
- - 237,671Г0+ 475,3437;/2- 237,6717; = 197,299,
29,70897’o - 237,67\Ty2+ 207,9627; = 94,5598;
- коэффициенты и правая часть системы уравнений для сегмента [л/4,я/2]:
|
тг/2 |
х( х ) ¥ сЬ с = -1 £ - =207,962, |
|
|
|
||||
|
Ш |
|
|
|
|||||
|
J* |
|
|
за |
|
|
|
|
|
|
л/41 |
|
|
|
|
|
|
|
|
|
л/4 |
я/4 |
|
|
о ) |
|
|
|
|
|
|я.<р;(л> ; /2(х)а!г = |
|Яф'/2(х>р;(х>&= |
=-237,671 |
||||||
|
л /2 |
|
|
|
|
|
|
|
|
|
{^[ф'з/2 |
^ = |
3/7 |
= 475)343> |
|
|
|
||
|
я/4 |
|
|
|
|
|
|
||
|
л/2 |
|
л/4 |
|
|
|
|
|
|
|
JX(p'(х)ф2(x)dx = |
|Хф2(х)ф' (х)Л = — = 29,7089, |
|
||||||
|
f |
|
|
|
|
|
|
|
|
|
л/4 |
|
|
|
|
|
|
|
|
|
ir/2 |
|
|
|
|
|
|
|
|
|
|Х[ф'2(*)]2Л = ^ |
= 207,962, |
|
|
|
||||
|
fJ |
|
|
ЗЛ |
|
|
|
|
|
|
л/4 |
|
|
|
|
|
|
|
|
|
л/2 |
л,/2 |
|
|
о* |
|
|
|
|
|
J ^ ' /2(J:V2(X>/X = |
| ^ - ф2(х) фз/2{х)dx = - — = -237,671; |
|||||||
|
л/4 |
л/4 |
|
|
3А |
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
cos А + 4cos2A = 96,0219, |
|||
|
Ф |
|
|
|
|
|
|
|
|
|
W0 jsinx(p3/2(x)dx = fV0 8(cosA-cos2A) |
4(sin A + sin 2A) |
= 476,323, |
||||||
|
л/4 |
|
h1 |
|
|
A |
|
|
|
Л/2 |
|
|
|
( 4 |
Л |
4 cos A |
|
||
W0 |
(sinA + 3sin2A) |
|
= 134,762; |
||||||
Jsinхф2{x)dx = IV0 ------------------i + |
—r - 1 cos2A- |
|
|||||||
я/4 |
h |
|
U |
) |
|
|
|
207,9627; - 237,671Г3/2 + 29,70897; = 96,0219, -237,6717; +475,3437^2 -237,6717; = 476,323, 29,70897; - 237,671Гу2 + 207,9627; = 134,762;
Зл/2
}Х[<р/2(х)Г ^ = ^ = 207,962,
!'•
|
|
л/2 |
|
|
|
|
|
Зл/2 |
Зл/2 |
= ~ |
=-237,671, |
||
|
}М>2(xV5/2 |
= |
JA.<P;/2 (дг)ф2 |
|||
|
я/2 |
"/2 |
ЗА |
|
||
|
|
|
|
|||
|
|
Зл/2 |
|
|
|
|
|
|
JA.[(P;/2(JC)]2^ = ^ |
= 475,343, |
|
|
|
|
|
л/2 |
|
|
|
|
|
Зл/2 |
|
Зя/2 |
|
|
|
|
jW 2 (л:)<Рз(*>& = |
| A.<PJ(дг)ф2 (x)dx= = 29,7089, |
|
|||
|
л/2 |
|
Л/2 |
зл |
|
|
|
|
|
|
|
||
|
|
Зл/2 |
|
|
|
|
|
|
}^[(рз(д:)]2 C*X = TZ~ 207,962, |
|
|
||
|
|
л/2 |
3/7 |
|
|
|
|
|
|
|
|
|
|
|
Зя/2 |
Зя/2 |
|
|
|
|
|
|я.<р;/2(х>р;(*> &= JA -Фз(*)q>;/2( * > & = - — |
= - 237,671; |
||||
|
Я/2 |
л/2 |
|
|
|
|
Зл/2 |
3sin2/J4-sin3/z ( 4 |
Л .. |
4cos3A |
= 134,762, |
||
0о |
Jsin хф2 (x)ir = W0 |
A |
---- - T - l cos2/7+ |
А |
||
Ф |
U |
J |
|
|||
|
|
|
|
|
||
|
Зл/2 |
|
|
4(sin 2/7+ sin ЗЛ) = 476,323, |
||
^0 |
jsin;c(p5/2(*)<fc= 0'o 8(cos 2ЛcosЗ/7 ) |
|||||
|
л/2 |
|
|
|
|
|
Зл/2 |
|
|
|
4cos2/7 =96,0219; |
||
Jsinjc(p3(jc)i!c= ^0 (sin 2Л+ 3sinЗЛ) ( f 4 |
Л со„эд |
|||||
Ф |
A |
+1,AI _ JC0S |
A5 |
|
||
|
|
|
|
|
207,9627, -237,67175/2 + 29,70897, = 134,762, -237,6717,+475,34375/2 -237,67173 = 476,323, 29,708972 - 237,6717]V2 + 207,96273= 96,0219;
|
|
|
= ^ |
= 207,962, |
|
|
||
|
|
Зя/2 |
|
ЗЛ |
|
|
|
|
|
JX(p;(jc)q>7/2(x)dx = JX(p'7/2(дс>рз |
= ~ |
= -237,671, |
|||||
|
Зл/2 |
|
Зл/2 |
|
|
|
|
|
|
|
№ /2 W f ^ = ^ |
|
= 475,343, |
|
|||
|
|
Зл/2 |
|
3,7 |
|
|
|
|
|
К (jc>p4(•*■)dx = /я,ф' (х)ф'з(х>2г = ^ - = 29,7089, |
|
||||||
|
Зл/2 |
|
Зл/2 |
|
|
|
|
|
|
|
3^2 |
Л = - ^ = 207,962, |
|
|
|||
|
|
|
3/1 |
|
|
|
||
|
) Н /2(*M (x]dx= |
|Лф^(х )ф; 2 (x)dc = ~ |
= -237,671; |
|||||
|
Зл/2 |
|
Зл/2 |
|
|
|
|
|
^0 {зтхфз(зс>& = Ж0 3sin З/г+ sin 4/z |
р --1 |cos3/* + 4cos4Л |
= 94,5598, |
||||||
Зл/2 |
|
|
h |
|
|
|
_ fc5 - |
|
^0 Jsinдхр7/2 (J:)£*C= |
8(cos3A-cos4A) |
4(sin3A + sin4A) |
|
|||||
|
|
|
|
|
= 197,299, |
|||
Зл/2 |
|
|
|
|
|
|
|
|
Wg Jsin |
хф4 (x)t& = W0 |
(sin3A + 3sin4A) |
+ |
4 |
'j |
|
= 1,03391. |
|
|
A |
—5--I cos4A- 4cosЗА |
||||||
3*/2 |
|
|
|
|
|
A* |
|
207,9627’ -237,6717^+29,70897; =94,5598 «-237,6717; + 475,343Г7/2- 237,67Щ = 197,299, 29,70897; - 237,671Г7/2 + 207,962Г4 = 1,03391.
Ансамблирование систем, полученных для каждого сегмента, приводит к формированию общей для всего стержня системы линей ных алгебраических уравнений (табл. 2.4) относительно коэффициен тов Ti приближённого решения (2.9) дифференциального уравнения.
Учёт граничного условия первого рода Т(х]х=о = Т° на левом
торце стержня и граничного условия третьего рода ХТ'(х)|м =
= _а [r(jc)|jr_Jt- 7 ,ee] на правом конце стержня приводят к системе ли
нейных алгебраических уравнений, приведённой в табл. 2.3.
Результаты расчетов
Использование Программы 2.2 даёт решение этой системы уравнений:
Т0=100,0 7^=102,119, 7] =103,407, Г3/2= 103,154, Т2=100,897,
Г5/2=96,4594, Тъ=90,018, 7^=82,0351, Т4=73,2219.
Приближённые решения дифференциального уравнения на отрез ке [0, тс], полученные на основе кусочно-квадратичного представления (2.9), приведены на рис. 2.4 с использованием 4 (а) и 64 (б) сегментов. Полученное приближённое решение даёт возможность определить по токи тепла с торцевых поверхностей стержня. Первое и последнее уравнения
|
|
— Г |
+ — Т = W sin/? |
|
4cosh |
- Q \ |
|
||
|
|
ЗА |
1/2 |
ЗА |
1 0 ~h~ |
|
|
|
|
— T |
- - T |
+ — T =w sin(7i-/?)_f 4 |
Л 4cos(n-h) |
- QL, |
|||||
3/7 |
3/7 |
m~]/1 |
3h |
m |
0 |
|
|
|
|
исключённые ранее из системы, преобразуются к виду |
|
|
|||||||
|
Q " - w . |
|
|
|
- A ( , r . _ Щ г + ^ ), |
|
|||
П1 _ w ^[sin(ji - А)+ А]- 4[l + C O S (TC - |
А)] |
X , |
|
ч |
|||||
V |
"о |
|
|
^2 |
Ж т~' ~ |
+ 1Тт'' |
ГО
ГО
Таблица 2.4
Матрица коэффициентов и правая часть системы линейных алгебраических уравнений для построения приближённого решения дифференциального уравнения с использованием кусочно-квадратичных пробных функций
1.0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
100,0 |
-237,671 |
475,343 |
-237,671 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
197,299 |
29,7089 |
-237,671 |
415,925 |
-237,671 |
29,7089 |
0,0 |
0,0 |
0,0 |
0,0 |
190,582 |
0,0 |
0,0 |
-237,671 |
475,343 |
-237,671 |
0,0 |
0,0 |
0,0 |
0,0 |
476,323 |
0,0 |
0,0 |
29,7089 |
-237,671 |
415,925 |
-237,671 |
29,7089 |
0,0 |
0,0 |
269,523 |
0,0 |
0,0 |
0,0 |
0,0 |
-237,671 |
475,343 |
-237,671 |
0,0 |
0,0 |
476,323 |
0,0 |
0,0 |
0,0 |
0,0 |
29,7089 |
-237,671 |
415,925 |
-237,671 |
29,7089 |
190,582 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
-237,671 |
475,343 |
-237,671 |
197,299 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
29,7089 |
-237,671 |
237,963 |
601,034 |
а
Рис. 2.4. Приближённое решение уравнения A,r'(jt)+ W(x)= 0 с граничными условиями r(jc)|j=0= Т°
и ХТ’( х ) [ т ( х ) х=1t - Г,] при использовании 4 (а) и 64 (б)
сегментов на отрезке [0, л]
Точное решение (1.18) поставленной задачи позволяет опреде лить точные значения потоков на торцах стержня
\(w0-a{T°-T_)) |
Х((У0-а(т°-Т„)} |
Q0=W{x\__0 = WQCOSX + |
= W0+ |
X+an |
X+ an |
X(W0-a{T°-T„ |
Uw0-a{T°-T_)) |
QL =-\T\X\^ = A ^T0 COSJC + - |
= W0- |
X+an |
X+ an |
и оценить погрешность полученных приближённых значений тех же потоков.
Приближённые значения тепловых потоков, полученных при раз личном числе сегментов, приведены в табл. 2.5. Значения потоков, по лученные с использованием точного решения: Q0 =403,34 Вт/м2
и QL =1596,66 Вт/м2
Таблица 2.5
Приближённые значения тепловых потоков на торцах стержня при различном числе т сегментов
т |
<? |
& |
т |
& |
Qt |
2 |
403,341 |
1596,66 |
16 |
403,340 |
1596,66 |
4 |
403,341 |
1596,66 |
32 |
403,341 |
1596,67 |
8 |
403,339 |
1596,66 |
64 |
403,357 |
1596,69 |
В табл. 2.6 приведены значения погрешностей |
|
|
8И = ЦТ ~ТЯI= тах|г(х) - Тя(х) = п^ахШх) - £ |
(х1, |
|
= 1тгп,~тЛ =maxiT2m(х) - |
Tm(х)| = max X 7].(pf (х) - £ |
Г,<р,(х) |
д е [ ° д ; |
х е [ 0 , п \ |
|
определённых с помощью точного решения Т(х) дифференциаль ного уравнения и сравнением двух последовательных приближён ных решений Тт и Т2т дифференциального уравнения при различ ных значениях числа т и 2т сегментов с использованием чебышёвской нормы. Эти же данные представлены в графическом виде на рис. 2.5, б.
Таблица 2.6
Погрешность Ътприближённого решения дифференциального уравнения при различном числе т сегментов
т |
h |
5„, |
& т .2т |
т |
h |
5m |
^ т.2т |
2 |
1,570796 0,316257 |
0,301401 |
16 |
0,196350 |
0,000859 |
0,000813 |
|
4 |
0,785398 0,051170 |
0,049652 |
32 |
0,098175 |
0,000104 |
0,000150 |
|
8 |
0,392699 0,006753 |
0,006613 |
64 |
0,049087 |
0,000013 0,000015 |
Рис. 2.5. Погрешность приближённого решения дифференциально го уравнения, полученного на четырёх сегментах с использованием пробных кусочно-квадратичных функций (я); зависимость по грешностей 8Ш(—о—)и 52т(-А-) приближённого решения
от длин сегментов (б)
С применением формулы (В.1) приближенно определяется поря док погрешности численного решения дифференциального уравне ния разложением (2.9) по системе кусочно-квадратичных пробных функций. Для 8тпорядок погрешности оценивается значением
Ът= (In0,316257 - In0,000013)/(ln 1,570796In0,049087) = 2,91405.
Для S2mэта величина
Ьт2т= (1п 0,30140 - In0,000015)/(1п1,570796In0,049087) = 2,85888.
Погрешность полученного решения дифференциального уравне ния методом Галёркина с использованием системы кусочно-квад ратичных пробных функций в обоих случаях можно приближённо оценить как величину третьего порядка относительно длины сегмен та А, т.е.
I |
т II |
/ \ |
II 2т |
т |
| | / \ |
|
г-х?;.<р, |
- о И , бт.2я= |
|
- Е Ы |
- ° И - |
||
I |
i=0 |
|| |
|
|| 1=0 |
1=0 |
|| |
В силу этого 8т —>0 при И—> 0 или т -> «>. Это позволяет утвер ждать, что последовательность приближённых решений дифференци
ального уравнения, полученных аппроксимацией кусочно-постоянны ми функциями (2.9), сходится равномерно на отрезке [0, к].
Выполненные расчёты (см. табл. 2.5, рис. 2.5, б) показывают, что погрешности 8т и 6^ практически одинаково оценивают погреш ность приближённого уравнения. Это означает, что при отсутствии точного решения, когда погрешность 8топределить невозможно, для контроля погрешности можно применять оценку Ььп-
На рис. 2.6 приведена зависимость времени t выполнения расче тов от числа т сегментов разностной сетки.
Рис. 2.6. Зависимость времени / выполнения расчетов от числа т сегментов разностной сетки
Выводы
1. Процедура метода Галёркина использована для приближённо го решения дифференциального уравнения. Сформирована система линейных алгебраических уравнений для определения коэффициен тов разложения искомого решения по системе кусочно-квадратичных функций.
2.Разработана вычислительная программа определения коэффи циентов разложения решения дифференциального решения по сис теме кусочно-квадратичных функций.
3.С использованием разработанной программы определены ко эффициенты и построены приближённые решения дифференциаль ного уравнения для 2, 4, 8, 16, 32 и 64 сегментов постоянной длины (см. рис. 2.4).
4.Для указанной последовательности разложений определены оценки Ьт и Ььп погрешности приближённого решения (см. табл. 2.5) для различного числа т сегментов.
5.Показано, что с уменьшением длины h сегментов погрешность приближённого решения дифференциального уравнения, определяе мая чебьпиёвской нормой, уменьшается (см. рис. 2.5, б). Установлено, что погрешность аппроксимации имеет третий порядок относительно длины h сегментов (шага интегрирования).
6.Выполненное исследование показывает, что последователь ность приближённых решений дифференциального уравнения, полу ченных методом Галёркина при аппроксимации кусочно-квадра
тичными функциями, сходится равномерно на отрезке [0, 7с].
7. Для численного решения дифференциального уравнения ме тодом Галёркина с использованием кусочно-квадратичных функций на разностной сетке, содержащей 64 сегмента, на компьютере с про цессором Intel® Pentium® 4 (тактовая частота 2,2 ГГц, объем опера тивной памяти 512 Мбайт) требуется 2,5-10-3 с.
2.3.Аппроксимация решения иерархической системой кусочно-непрерывных полиномов 1-й степени
Задание. Методом Галёркина с использованием иерархической системы пробных функций построить приближённое решение од номерного дифференциального уравнения стационарной теплопро водности [ХТ'(х)\ + JT0sinjt = 0 с граничными условиями r(jt)|j=0= f
А,Г'(х)|^ = - a |V ( x ) |^ - J ^ j . Сформировать систему линейных ал
гебраических уравнений относительно коэффициентов разложе ния искомого решения по заданнной системе пробных функций; разработать вычислительную программу для определения коэф фициентов разложения решения дифференциального уравнения по заданнной системе пробных функций для 2, 4, 8, 64 сегментов постоянной длины; при известном точном решении определить погрешность приближённых решений для указанной последователь ности сегментов; исследовать зависимость погрешности искомого решения от длины h сегментов; исследовать сходимость процесса аппроксимации; оценить быстродействие вычислительной програм мы. При выполнении расчётов принять: L = 7Ц А = 70 Вт/м •град, \¥0=\Ш Вт/м3, а = Вт/м2град, Т° = 100й, = 20°
Разрешающие соотношения
Стержень, имеющий длину L = n, разбивается на т сегментов JCJ, / = 1,/и, равной длины h = xi - x i l = к1т. Рассматривается про цесс аппроксимации решения дифференциального уравнения с ис
пользованием двух полиномов первой степени:
Фы М = (х, - x )/h> Ф,(*) = (* - Xi-I)/h ■
Приближённое решение на основе кусочно-линейных пробных функций строится в виде
ТАх)=Ti-\ФмМ + т’/Ф,- W •
Построение системы линейных алгебраических уравнений для определения коэффициентов подробно описано в подразд. 2.1.1. Для сегмента [дг/_\,х,] эта система линейных алгебраических уравнений имеет вид
,sin х —sin Х |
=W01cos*,.,------------------
T, = Wn SinX; - sin X,7-1 - COSJC; \-Vr
Алгоритм решения
Фрагмент программного кода на языке Си, реализующий проце дуру построения приближённого решения дифференциального урав нения с использованием иерархической системы кусочно-непрерыв ных полиномов 1-й степени, представлен Программой 2.1.
Реализация алгоритма
Отрезок [0, к] разбивается на четыре сегмента равной длины
[0,71j = [0,я/4] U [я/4, л/2]U [я/2,Зтх/4] U [3л/4, я].
На каждом из этих сегментов в соответствии с выражением (2.11) определяются кусочно-непрерывные полиномы 1-й степени. Решение дифференциального уравнения представляется разложени ем (2.9). Для построения приближённого решения формируется сис тема линейных алгебраических уравнений (см. табл. 2.1) относитель но искомых коэффициентов Тп / = 0,4.
Результаты расчётов
Результаты приближённого решения заданного дифференциаль ного уравнения с использованием иерархической системы кусочно непрерывных полиномов 1-й степени приведены на рис. 2.1-2.3 и в табл. 2.2-2.3.
Выводы
1.Процедура метода Галёркина использована для приближён ного решения дифференциального уравнения. Сформирована сис тема линейных алгебраических уравнений для определения коэф фициентов разложения искомого решения по системе кусочно линейных функций.
2.Разработана вычислительная программа определения коэф фициентов разложения решения дифференциального решения по системе кусочно-линейных функций.
3.С использованием разработанной программы определены ко эффициенты и построены приближённые решения дифференциаль ного уравнения для 2, 4, 8, 16, 32 и 64 сегментов постоянной длины (см. рис. 2.2).
4.Для указанной последовательности разложений определены
оценки 5т и Ъгт погрешности приближённого решения (см. табл. 2.3) для различного числа т сегментов.
5. Показано, что с уменьшением длины h сегментов погреш ность приближённого решения дифференциального уравнения, оп ределяемая чебышёвской нормой, уменьшается (см. рис. 2.2, б). По грешность аппроксимации имеет первый порядок относительно длин сегментов (шага интегрирования).
6. Выполненное исследование показывает, что последователь ность приближённых решений дифференциального уравнения, полу ченных методом Галёркина при аппроксимации кусочно-линейными функциями, сходится равномерно на отрезке [0, л].
7. Для численного решения дифференциального уравнения ме тодом Галёркина с использованием кусочно-линейных функций на разностной сетке, содержащей 64 сегмента, на компьютере с процес сором Intel® Pentium® 4 (тактовая частота 2,2 ГГц, объем оператив ной памяти 512 Мбайт) требуется 2,5-1 (Г3с.