Визначення стаціонарного температурного поля у пластині з кільцевим перерізом при заданих граничних умовах першого роду.
Société Française des Mécaniciens, Commission Validation des progiciels de calcul de Structures, groupe de travail Thermique (2D et 3D) et thermoélasticité, Paris, 1989.
Визначити радіальний розподіл температур у тілі пластини, а також значення температур у характерних точках (вузлах) по радіусу.
Розрахункова модель являє собою пластину перерізу порожнистого циліндра у формі кільця, із ізотропного матеріалу з коефіцієнтом теплопровідності.
До граней окружності циліндра на внутрішньому та зовнішньому радіусі прикладені граничні умови першого роду та тепловий потік по внутрішньому контуру.
|
а |
б |
Внутрішній радіус Ri = 0,3 м;
Зовнішній радіус Re = 0,35 м;
Довжина циліндра l = 1 м.
Коефіцієнт теплопровідності λ = 1 Вт/м°C.
По внутрішньому радіусу:
Задана температура Тi = 100 °С.
По зовнішньому радіусу:
Задана температура Тe = 20 °С.
Задача розв’язується в просторовій схемі загального виду з 1 ступенем вільності у вузлі – температура t (ознака схеми 15).
Для побудови схеми використано СЕ 1509 – чотирикутні СЕ теплопровідності.
Сітка скінченних елементів має 5 елементів вздовж радіусу і 145 вздовж окружності.
Кількість вузлів: 870. Кількість елементів: 725.
|
Мозаїка температур T, °C |
Мозаїка теплового потоку, Дж/с |


φ = Φ/S
де Φ - тепловий потік, (1Вт = 1 Дж/c);
φ – щільність теплового потоку, (Вт/м2 = 1 Дж/(c*м2 )).
| Шукана величина | Радіус R, м | Аналітичне рішення | LIRA-FEM | Похибка, % |
| T, °C | 0,3 | 100 | 100 | 0 |
| 0,31 | 82,98 | 82,983 | 0,004 | |
| 0,32 | 66,51 | 66,507 | 0,005 | |
| 0,33 | 50,54 | 50,537 | 0,006 | |
| 0,34 | 35,04 | 35,044 | 0,011 | |
| 0,35 | 20 | 20 | 0 | |
| Φ/l, Вт/м | [0,3/0,35] | 3260,8 | 3261,573 | 0,024 |
Помилка в тексті? Виділіть її та натисніть Ctrl + Enter, щоб повідомити нам.
Коментарі