Energy
education

сайт для тех, кто хочет изучать энергетику

Контрольные задания №1. Задание 5

Пластина толщиной $\delta = 60$ мм и коэффициентом теплопроводности $\lambda = 65$ Вт/м·К с равномерно распределенными внутренними источниками теплоты $q_v = 10^6$ Вт/м3 охлаждается с обеих сторон. С одной стороны на пластину наложена изоляция толщиной $\delta_{из} = 4$ мм и коэффициентом теплопроводности $\lambda_{из} = 12$ Вт/м·К . Температура охлаждающей жидкости $t_ж = 150$ °C и коэффициент теплоотдачи $\alpha = 10^4$ Вт/м2·К. Условия с неизолированной стороны - теплоотдача

Определить максимальную температуру, температуры в центре, на поверхностях пластины и изоляции.

Решение:

Для упрощения задачи путем учета изоляции в качестве дополнительного термического сопротивления теплоотдачи, рассчитаем приведенный коэффициент теплоотдачи, Вт/м2·К:

$$\alpha_{пр} = \frac{1}{\frac{1}{\alpha}+\frac{\delta_{из}}{\lambda_{из}}} = \frac{1}{\frac{1}{10^4}+\frac{0.004}{12}} = 2300.$$

При указанных условиях температура пластины будет изменяться только вдоль оси x, тогда дифференциальное уравнение принимает вид:

$$\frac{\partial^2t}{\partial x^2}+\frac{q_v}{\lambda} = 0.$$

После интегрирования получим:

$$\frac{\partial t}{\partial x} = -\frac{q_v·x}{\lambda} + C_1;$$ $$t = -\frac{q_v·x^2}{2·\lambda}+C_1·x+C_2.$$

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

при $x = 0$

$$\lambda·\frac{\partial t}{\partial x} = \alpha·(t-t_ж).$$

при $x = \delta$

$$-\lambda·\frac{\partial t}{\partial x} = \alpha_{пр}·(t-t_ж).$$

Тогда:

$$\lambda·\left(-\frac{q_v·0}{\lambda}+C_1\right)=\alpha·\left(\left(-\frac{q_v·0^2}{2·\lambda}+C_1·0+C_2\right)-t_ж\right);$$ $$-\lambda·\left(-\frac{q_v·\delta}{\lambda}+C_1\right)=\alpha_{пр}·\left(\left(-\frac{q_v·\delta^2}{2·\lambda}+C_1·\delta+C_2\right)-t_ж\right).$$

Или

$$65·\left(-\frac{10^6·0}{65}+C_1\right)=10^4·\left(\left(-\frac{10^6·0^2}{2·65}+C_1·0+C_2\right)-150\right);$$ $$-65·\left(-\frac{10^6·0.06}{65}+C_1\right)=2300·\left(\left(-\frac{10^6·0.06^2}{2·65}+C_1·0.06+C_2\right)-150\right).$$

Решая систему находим: $C_1 = 567.2$, $C_2 = 154$.

Задаем распределение температуры внутри пластины:

$$t(x) = -\frac{q_v·x^2}{2·\lambda}+C_1·x+C_2= -\frac{10^6·x^2}{2·65}+567.2·x+154.$$

Находим температуры в середине и на краях пластины:

при $x = 0$ м температура равна $t(0) = 153$ °C;

при $x = \delta/2 = 0.03$ м температура равна $t(0.03) = 164$ °C;

при $x = \delta = 0.06$ м температура равна $t(0.06) = 160$ °C.

Находим максимальную температуру используя свойства экстремумом:

при $\partial t/\partial x = 0$ или $0 = -10^6·x_{max}/65+567.2$ находим $x_{max} = 0.037$ м;

при $x = x_{max} = 0.037$ м температура равна $t_{max}(0.037) = 164.1$ °C.

Находим температуру поверхности изоляции:

$$\frac{\lambda_{из}}{\delta_{из}}·(t(\delta)-t_{из})=\alpha·(t_{из}-t_ж);$$ $$\frac{12}{0.004}·(164-t_{из})=10^4·(t_{из}-150).$$

Откуда температура изоляции равна $t_{из} = 152$ °С.

Пример решения в MathCAD:

Пример решения в ANSYS Steady-State Thermal:

Создаем модель

Создаем сетку

Задаем внутренние источники тепла

Задаем граничные условия конвективного теплообмена с обоих сторон

Получаем распределение температур