Рішення

На першому етапі рішення задачі застосовується до функції Т(х,t) пряме перетворення Фур'є по перемінній х

T(x,t) ®F(s,t)= (4)

де T(x,t) – оригінал (шукана функція); F(s,t) зображення для функції
T(-¥ < s< + ¥ - деякий параметр).

Метою даного етапу є перебування зображення F(s,t) за допомогою перетвореної системи рівнянь (1) – (2) відповідно до формули (4). Використовуючи формули, аналогічні (4), одержуємо для лівої і правої частин рівняння (1), а також початкової умови (2) наступні перетворення Фур'є

, (5)

(6)

Застосовуючи до інтегралу, що розміщений в правій частині (6), двічі метод інтегрування по частинам і враховуючи граничні умови (3), одержуємо, що

(7)

Тоді

(8)

Для початкової умови (2) маємо

(9)

Таким чином, враховуючи (5) - (9) для знаходження зображення F(s,t) одержуємо звичайне диференціальне рівняння 1-го порядку з розподілювальними перемінними (замість диференціального рівняння в частинних похідних 2-го порядку для функції Т(х,t)

(10)

з початковою умовою

(11)

Рішення рівняння (10) з урахуванням (11) одержуємо прямим інтегруванням лівої і правої його частин

(12)

На другому етапі рішення задачі для знаходження шуканого рішення Т(х,t) використовується зворотне перетворення Фур'є

(13)

Підставляючи (12) у (13), одержуємо

(14)

З урахуванням (9) для визначення T(x,t) одержуємо наступний подвійний інтеграл

(15)

I

 


Для обчислення внутрішнього інтеграла I скористаємося відомою формулою Ейлера

(16)

Одержуємо

(17)

I'

 

 


Враховуючи, що I¢ = 0 (інтеграл від непарної функції із симетричними межами), маємо

(18)

Вираз для інтегралів типу (18) можна знайти в спеціальних довідниках[1]

(19)

Підставляючи (19) у (15), одержуємо наступне рішення вихідної задачі

(20)

Як приклад використання формули (20) розглянемо окремий випадок із заданою функцією j(х) у вигляді деякого прямокутного імпульсу

(21)

де х – у м, j(х)- у °С.

Підставляючи (21) у (20), маємо

(22)

Зведемо інтеграл, що розміщений в правій частині (18), до спеціальної функції (інтеграл помилок, інтеграл Лапласа, інтеграл імовірностей)

erf(z) = , (23)

значення якої для різних значень аргументу z можна знайти в довідковій літературі. Для цього зробимо заміну перемінних в інтегралі (23)

(24)

Одержуємо

(25)

З урахуванням (23) одержуємо наступну формулу для розподілу температури в розглянутому середовищі в будь-який момент часу

(26)

де Т - у °С; t- у с.

Задача 2 Процес нагрівання виробу рухливим джерелом теплового впливу описується наступним рівнянням теплопровідності (рис.3, рухлива система координат).

Рисунок 3 – Схема теплової дії на виріб ( - вузький тонкоплівочний елемент, тонкий, довгий стрижень та інше)

 

(27)

початковою умовою

(28)

і граничними умовами

при х®±¥, (29)

де V - швидкість руху джерела; Fv(x,t) - об'ємна густина теплового потоку;
Т0 = const - початкова температура виробу. Потрібно, використовуючи метод інтегральних перетворень Фур'є, знайти загальне рішення задачі T(x,t).

Рішення.

Шляхом відповідної заміни перемінних перетворимо систему рівнянь
(27) – (29) до стандартного (фундаментального) виду, а потім для рішення отриманих рівнянь використаємо метод інтегральних перетворень Фур'є. Спочатку зробимо заміну

(30)

Одержуємо

(31)

(32)

при х ®±¥ (33)

Далі виключимо з рівняння (31) член , зробивши заміну перемінних

(34)

Підставляючи (34) у (31) - (33), одержуємо для функції уже фундаментальну систему рівнянь

(36)
(35)

(38)
при х ® ±¥, (37)

Як і при рішенні першої задачі, для знаходження функції на 1-муетапі використовуємо пряме перетворення Фур'є по координаті х

(39)

Також як і при рішенні попередньої задачі, знайдемо F( , перетворивши (35) - (37) відповідно до (39)

(40)

(41)

(42)

Використовуючи (40)-(42), для знаходження F(s,t) одержуємо звичайне лінійне диференціальне рівняння 1-го порядку, яке уже можна вирішити одним із стандартних методів (наприклад, методами Бернуллі, Ейлера та ін.)

(43)

(44)

Знайдемо рішення (43) - (44) методом Бернуллі

(45)

де m(t) і n(t) – дві невідомі функції, що знаходяться в процесі рішення
(43) – (44). Підставляючи (45) у (43), одержуємо

(46)

Функцію n(t) знаходимо, спростивши (46)

(47)

Інтегруючи (47), (константу в рішенні опускаємо) одержуємо

(48)

З урахуванням (46) і (48) знаходимо

(49)

Звідси

, (50)

де с – деяка постійна. Підставляючи (48) і (50) у (45), одержуємо

(51)

З урахуванням початкової умови (44), остаточно, для F(s,t) одержуємо наступний вираз

(52)

На другому етапі рішення задачі для знаходження оригіналу скористаємося зворотним перетворенням Фур'є для функції F(s,t)

(53)

Підставляючи (42) і (52) у (53), одержуємо

(54)

Згрупувавши інтеграли в (54) і розбивши їх на внутрішній і зовнішній, знаходимо

(55)

I'

 

 


Також як і в попередній задачі, для знаходження внутрішнього інтеграла використовується (16) і довідкова література

(56)

Підставляючи (56) у (55), одержуємо наступне рішення для функції

(57)

Використовуючи (30), (34) і (57), остаточно, для рішення вихідної задачі маємо

(58)

Розглянемо джерело теплового впливу з нормально розподіленою потужністю (гаусового типу)

, (59)

де Р(t) – потужність джерела в центрі впливу; k(t) коефіцієнт зосередженості джерела; Н, В – відповідно товщина і ширина елемента. Підставляючи (59) у (58) і враховуючи що Р(t) =Р0 і k(t) =k00 і k0 - постійні величини), одержуємо

(60)

Формула (60) дозволяє проводити розрахунки розподілів температури в різних тонкоплівкових виробах машинобудування, точного приладобудування, космічної й авіаційної техніки та ін. при їхній обробці рухливими концентрованими джерелами тепла (різець, фреза, електронний промінь, лазерне випромінювання та ін).









Дата добавления: 2015-08-21; просмотров: 532;


Поиск по сайту:

При помощи поиска вы сможете найти нужную вам информацию.

Поделитесь с друзьями:

Если вам перенёс пользу информационный материал, или помог в учебе – поделитесь этим сайтом с друзьями и знакомыми.
helpiks.org - Хелпикс.Орг - 2014-2023 год. Материал сайта представляется для ознакомительного и учебного использования. | Поддержка
Генерация страницы за: 0.02 сек.