книги / Прогнозирование теплового состояния изделий при эксплуатации в условиях воздействия солнечного излучения
..pdf6.5. Расчет теплового состояния изделия И5
Изделие И5 представляет собой полый цилиндрический полимерный элемент, покрытый по наружной и внутренней поверхности теплозащитным материалом (резиной) и размещенный среди двух цилиндрических обечаек. Изделие вместе с цилиндрическими обечайками находится в контейнере. Схема размещения конструкционных элементов изделия показана на рис. 6.10. Вся конструкция представляет собой многослойный цилиндр. Геометрические размеры и теплофизические характеристики материалов изделия приведены в табл. 6.7.
Таблица 6 . 7 Геометрические и теплофизические характеристики материала
|
Материал |
Толщина, |
λ, |
с, |
ρ, |
|
мм |
Вт/м·К |
кДж/кг·К |
кг/м3 |
|
1. |
Контейнер (стеклопластик) |
3,5 |
0,124 |
1,382 |
1360 |
2. |
Воздушный зазор |
2,0/1,37 |
0,0175 |
1,009 |
1,252 |
3. |
Корпус (обечайка) |
3,0 |
99,65 |
0,938 |
288 |
4. |
ТЗП (пластмасса) |
1,5 |
0,213 |
1,507 |
1350 |
5. |
Воздушный зазор |
0,55 |
0,0175 |
1,009 |
1,252 |
6. |
ТЗП (резина) |
3,0 |
0,217 |
1,800 |
1400 |
7. |
Полимерный элемент |
37,0 |
0,198 |
1,985 |
1680 |
8. |
ТЗП (резина) |
3,0 |
0,217 |
1,800 |
1400 |
9. |
Воздушный зазор |
1,37 |
0,0175 |
1,009 |
1,252 |
10. ТЗП (пластмасса) |
2,0 |
0,120 |
5,359 |
1260 |
|
11. Корпус (обечайка) |
3,0 |
99,65 |
0,938 |
288 |
Требуется определить распределение температур по своду изделия для 1–3 июля при эксплуатации в естественных климатических условиях при воздействии температуры воздуха и солнечного излучения.
Для численного решения задачи методом конечных элементов выбрано поперечное сечение изделия. Расчетная схема дискретизации представлена на рис. 6.11, при этом параметры дискретизации аналогичны схеме на рис. 6.4, элементы 289–648 относятся к полимерным.
121
122
Рис. 6.10. Схема изделия И5: 1 – контейнер; 2, 5, 9 – воздушный зазор; 3 – корпус; 4 – ТЗП (пластмасса); 6, 8 – ТЗП (резина);
7 – полимерный элемент; 10 – ТЗП (пластмасса); 11 – корпус
Рис. 6.11. Схема дискретизации поперечного сечения изделия И5 (NY = 555; NE = 1008; SH = 17)
В качестве представительного пункта также выбран г. Асуан (северная Африка). Расчеты проводились для контейнера песочного цвета с коэффициентами ε = 0,45 и βc = 0,45. Шаг счета
также принимался равным ∆τсч = 0,1 ч.
Результаты расчетов теплового состояния изделия представлены в табл. 6.8.
123
Таблица 6 . 8
Распределение температуры в изделии И5 по своду и времени (г. Асуан, 1–4 июля, контейнер песочного цвета)
|
Среднесуточнаятемпературавозду- |
Среднесуточнаятемпературавозду- |
|||||
|
ха33,9 °С, (максимальная41,4 °С) |
ха43,6 °С, (максимальная51,1 °С) |
|||||
|
|
Номер узла в конечно |
-элементной сетке |
|
|||
Время, |
|
и место его расположения |
|
|
|||
ч |
472 |
477 |
480 |
472 |
|
477 |
480 |
(середина |
(середина |
|
|||||
|
(поверх- |
(поверх- |
|
(поверх- |
(поверх- |
||
|
свода по- |
ность кор- |
ность кон- |
свода по- |
|
ность кор- |
ность кон- |
|
лимерного |
пуса) |
тейнера) |
лимерного |
|
пуса) |
тейнера) |
|
элемента) |
|
|
элемента) |
|
|
|
1 |
2 |
3 |
4 |
5 |
|
6 |
7 |
0 |
30,0 |
30,0 |
30,0 |
30,0 |
|
30,0 |
30,0 |
2 |
29,54 |
28,54 |
27,17 |
31,49 |
|
33,31 |
33,08 |
4 |
28,59 |
27,14 |
25,40 |
32,42 |
|
33,23 |
33,87 |
6 |
27,62 |
27,53 |
25,12 |
32,90 |
|
34,57 |
33,97 |
8 |
28,49 |
30,15 |
27,22 |
34,87 |
|
37,88 |
36,34 |
10 |
30,86 |
34,58 |
35,52 |
38,12 |
|
42,93 |
45,05 |
12 |
34,91 |
40,77 |
51,79 |
42,88 |
|
49,63 |
61,54 |
14 |
39,91 |
46,56 |
62,56 |
48,45 |
|
55,81 |
72,42 |
16 |
44,36 |
49,83 |
65,44 |
53,35 |
|
59,35 |
75,36 |
18 |
47,18 |
50,49 |
61,85 |
56,51 |
|
60,22 |
72,85 |
20 |
46,53 |
42,98 |
39,16 |
55,97 |
|
53,01 |
49,62 |
22 |
43,24 |
39,81 |
36,07 |
53,13 |
|
50,02 |
46,58 |
24 |
40,52 |
37,03 |
33,09 |
50,61 |
|
47,33 |
43,60 |
26 |
37,83 |
34,24 |
30,19 |
48,04 |
|
44,60 |
40,69 |
28 |
35,20 |
31,72 |
27,85 |
45,50 |
|
42,13 |
38,38 |
30 |
32,90 |
31,15 |
27,09 |
43,28 |
|
41,59 |
37,62 |
32 |
32,65 |
32,91 |
28,73 |
43,02 |
|
43,25 |
39,16 |
34 |
34,09 |
36,67 |
36,46 |
44,43 |
|
47,00 |
46,85 |
36 |
37,39 |
42,36 |
52,40 |
47,73 |
|
52,71 |
62,70 |
38 |
41,81 |
47,77 |
62,99 |
52,17 |
|
58,15 |
73,25 |
40 |
45,83 |
50,76 |
65,77 |
56,20 |
|
61,14 |
75,99 |
42 |
48,31 |
51,20 |
62,11 |
58,70 |
|
61,59 |
72,34 |
44 |
47,22 |
43,56 |
39,44 |
57,66 |
|
54,12 |
50,13 |
46 |
43,93 |
40,29 |
36,30 |
54,46 |
|
50,90 |
46,99 |
48 |
41,07 |
37,40 |
33,28 |
51,65 |
|
48,02 |
43,93 |
50 |
38,27 |
34,54 |
30,34 |
48,87 |
|
45,15 |
40,95 |
52 |
35,55 |
31,96 |
27,98 |
46,15 |
|
42,57 |
38,59 |
54 |
33,18 |
31,35 |
27,20 |
43,79 |
|
41,93 |
37,79 |
56 |
32,87 |
33,05 |
28,80 |
43,42 |
|
43,51 |
39,28 |
58 |
34,26 |
36,78 |
36,50 |
44,74 |
|
47,20 |
46,93 |
60 |
37,52 |
42,44 |
52,43 |
47,97 |
|
52,86 |
62,76 |
62 |
41,92 |
47,84 |
63,02 |
52,35 |
|
58,27 |
73,29 |
64 |
45,91 |
50,81 |
65,79 |
56,34 |
|
61,23 |
76,02 |
124 |
|
|
|
|
|
|
|
Окончание табл. 6 . 8
1 |
2 |
3 |
4 |
5 |
6 |
7 |
66 |
48,37 |
51,24 |
62,12 |
58,80 |
61,66 |
72,36 |
68 |
47,27 |
43,59 |
39,46 |
57,74 |
54,18 |
50,15 |
70 |
43,97 |
40,31 |
36,32 |
54,53 |
50,94 |
47,01 |
72 |
41,10 |
37,42 |
33,29 |
51,70 |
48,06 |
43,95 |
74 |
38,29 |
34,56 |
30,35 |
48,91 |
45,18 |
40,97 |
76 |
35,57 |
31,98 |
27,98 |
46,18 |
42,59 |
38,60 |
78 |
33,20 |
31,36 |
27,21 |
43,82 |
41,95 |
37,79 |
80 |
32,88 |
33,06 |
28,80 |
43,44 |
43,52 |
39,29 |
82 |
34,27 |
36,78 |
36,49 |
44,75 |
47,20 |
46,93 |
84 |
37,53 |
42,45 |
52,43 |
47,98 |
52,87 |
62,76 |
86 |
41,92 |
47,84 |
63,03 |
52,36 |
58,27 |
73,30 |
88 |
45,91 |
50,81 |
65,80 |
56,35 |
61,24 |
76,03 |
90 |
48,38 |
51,25 |
62,13 |
58,81 |
61,67 |
72,37 |
92 |
47,28 |
43,60 |
36,46 |
57,75 |
54,18 |
50,15 |
94 |
43,97 |
40,31 |
36,32 |
54,53 |
50,95 |
47,01 |
96 |
41,11 |
37,43 |
33,29 |
51,71 |
48,06 |
43,95 |
По данным табл. 6.8 на рис. 6.12 и 6.13 представлены зависимости температуры от времени.
Лимитирующими срок эксплуатационной пригодности изделия И5 могут быть как процессы, протекающие в продукте, так
ипроцессы массообмена в зоне «продукт–покрытие». Для решения задачи обоснования допустимых сроков эксплуатации указанного изделия более подробно рассмотрены эквивалентные температуры эксплуатации самых теплонапряженных участков полимерного элемента, зоны «полимерный элемент – ТЗП» и «ТЗП». Рассчитанные исходные данные по годовому распределению температур
ипродолжительности их действия для определения соответствующих эквивалентных температур эксплуатации изделия приведены в табл. 6.9. На основании данных табл. 6.9 в табл. 6.10 представлены зависимости значений эквивалентной температуры от величины температурного коэффициента Е, лимитирующего процесс эксплуатации. Для всех трех зон получены достаточно близкие значения эквивалентной температуры.
Сравнение данных, представленных на рис. 6.9 и 6.14, показывает, что в прочих равных условиях благодаря более мощной теплозащите теплонапряженность зоны контакта «полимерный элемент–теплозащитное покрытие» для изделия И5 существенно ниже, чем для изделия А5.
125
126
127
Таблица 6 . 9
Годовое распределение температур и продолжительность их воздействия в т. 473 (поверхностный слой полимерного элемента), т. 474 (зона контакта «полимерный элемент– теплозащитное покрытие») и т. 475 (наружная поверхность теплозащитного покрытия), г. Асуан, изделие И5
Интервал |
Продолжительность действия температуры, ч |
||
температуры, °С |
|
|
|
т. 473 |
т. 474 |
т. 475 |
|
8–10 |
– |
14 |
28 |
10–12 |
135 |
159 |
188 |
12–14 |
362 |
302 |
312 |
14–16 |
250 |
306 |
237 |
16–18 |
344 |
418 |
394 |
18–20 |
479 |
413 |
457 |
20–22 |
422 |
395 |
410 |
22–24 |
469 |
449 |
504 |
24–26 |
552 |
530 |
482 |
26–28 |
526 |
546 |
474 |
28–30 |
576 |
590 |
683 |
30–32 |
667 |
728 |
705 |
32–34 |
815 |
717 |
686 |
34–36 |
592 |
573 |
613 |
36–38 |
516 |
521 |
523 |
38–40 |
361 |
458 |
358 |
40–42 |
515 |
383 |
391 |
42–44 |
352 |
363 |
422 |
44–46 |
350 |
306 |
268 |
46–48 |
254 |
254 |
295 |
48–50 |
223 |
304 |
244 |
50–52 |
– |
31 |
86 |
52–54 |
– |
– |
– |
128
Таблица 6 . 1 0
ЗависимостьзначенийэквивалентнойтемпературыТэкв отвеличинытемпературногокоэффициентаЕ, лимитирующего
процессдлязон: поверхностныйслойполимерногоэлемента(т. 473), контакт«полимерныйэлемент–ТЗП» (т. 474), наружная поверхностьТЗП(т. 475), изделиеИ5, г. Асуан
Е, |
|
Тэкв, К |
|
|
кДж/моль |
т. 473 |
т. 474 |
т. 475 |
|
(ккал/моль) |
||||
|
|
|
||
20,93 (5,0) |
304,29 |
304,38 |
304,39 |
|
41,87 (10,0) |
305,52 |
305,67 |
305,71 |
|
62,80 (15,0) |
306,67 |
306,88 |
306,94 |
|
83,74 (20,0) |
307,73 |
307,99 |
308,08 |
|
104,67 (25,0) |
308,69 |
309,00 |
309,11 |
|
125,60 (30,0) |
309,55 |
309,91 |
310,05 |
|
146,54 (35,0) |
310,33 |
310,74 |
310,89 |
|
167,47 (40,0) |
311,03 |
311,47 |
311,65 |
|
188,41 (45,0) |
311,66 |
312,14 |
312,3 |
|
209,34 (50,0) |
312,23 |
312,73 |
312,94 |
На рис. 6.14 зависимость Тэкв от температурного коэффициента Е для зоны контакта «полимерный элемент–ТЗП» (т. 474) показана графически.
Рис. 6.14. Зависимость Тэкв от температурного коэффициента Е для зоны контакта «полимерный элемент–ТЗП», изделие И5, контейнер песочного цвета
129
6.6.Обсуждение результатов расчета
Сприменением программы TеSunRad для двух изделий различных геометрических размеров и с различными теплофизическими характеристиками и компоновкой составляющих конструкционных элементов рассчитаны температуры, возникающие
вэлементах изделий в процессе их непрерывной эксплуатации
виюле в условиях повышенных температур окружающего воздуха и воздействия прямой солнечной радиации; рассчитаны зависимости для определения эквивалентных температур эксплуатации элементов и узлов изделий.
Анализ результатов расчетов показывает, что реализованная
впрограмме расчетная схема качественно правильно отражает влияние формирующих индивидуальное тепловое состояние изделия факторов: геометрических размеров, параметров теплообмена с окружающей средой, теплофизических характеристик конструкционных элементов, схемы компоновки изделия.
Что касается количественной оценки результатов расчетов, то следует полагать, что величина возможной расчетной ошибки определяется в основном тремя факторами:
–качеством расчетной схемы, реализующей решение системы уравнений нестационарной теплопроводности, моделирующей тепловое состояние изделия;
–точностью определения исходных характеристик элементов изделия;
–надежностью и обоснованностью назначения граничных условий решения задачи.
Первый фактор достаточно подробно рассмотрен в главах 4 и 5 настоящей работы, где определены условия устойчивости расчетной схемы и показано, что даже при относительно больших шагах по времени пространственной координаты расхождение между результатами расчетов тестовых задач по программе TeSunRad и апробированными методами не превышает 2 %.
Теплофизические характеристики входящих в состав изделий материалов определены принятыми в отрасли аттестованными методами либо взяты из справочной литературы и, как
130