De disco con una  temperatura  constante en la superficie cilíndrica externa se calienta lo largo del eje

Vamos a considerar un problema de dos dimensiones de una distribución  temperatura  de estado estacionario en la sección transversal-cara de extremo de un disco. Un calentador de espesor cero (una línea) con una potencia de P = 100 W está situado en el eje del disco, Y en la periferia una  temperatura  constante de 273.15 oK se lleva a cabo.

Parámetros de un disco: disco de metal con un espesor d = 5 mm, radio R = 100 mm y conductividad térmica  K=50 W/(m • oK) [aleación de acero (SS)].

(ver la figura).

La solución de este problema en el que se considera como la fuente distribuida se puede conseguir por resolver el problema de una fuente puntual de calor. Sea el origen del sistema de coordenadas coincide con el centro de la superficie circular de un disco. La ecuación diferencial que tiene que ser resuelto para una fuente puntual tiene la forma:

Donde ρ – densidad de una fuente de calor distribuida. Para nuestro caso:  ρ = P / d. Solución de esta ecuación es la función G del verde (función de fuente de calor)

En esta solución el punto cero corresponde al centro del disco. El eje Z del sistema de coordenadas, en el que se busca la solución, se dirige a lo largo del eje del disco.

El modelo de elementos finitos con condiciones de frontera aplicados

Comparemos solución analítica con la solución obtenida a partir de AutoFEM. Después de llevar a cabo el siguiente cálculo se obtiene los siguientes resultados:

Tabla 1. Parámetros de malla de elementos finitos

Tipo de elementos finitos.

Número de Nodos

Número de elementos finitos.

Tetraedro lineal

1263

3842

Tabla 2. Resultado " temperatura "

Distancia del centro, mm

Solución numérica.
temperatura T*, K

Solución analítica
temperatura T, К

Error δ = 100%* |T* - T| / |T|

R 30

350.7618

349.7972

0.27

R 40

331.7566

331.4828

0.08

R 50

317.1303

317.2771

0.05

R 60

305.7533

305.6701

0.03

 

Conclusiones:

Error relativo de la solución numérica. en comparación con la analítica solución no exceder de 0,3% cuando se utiliza elemento lineal. La exactitud del cálculo crece cuando se aleja de la fuente de calor debido a la singularidad de la analítica Solución en el punto de aplicación de fuente de calor (en el centro de un disco).

La trama de la dependencia de la  temperatura  en el radio muestra que las soluciones analíticas numéricas prácticamente coinciden. Esto implica que las distribuciones de máximos y mínimos   temperatura  son la misma, y por lo del tanto, el cálculo de los flujos de calor Térmico y de energía se lleva a cabo con el mismo grado de exactitud.

 

Los resultados de pruebas numéricas dependen de la malla de elementos finitos y pueden diferir ligeramente de los que figuran en la tabla.

 

Lea más acerca de AutoFEM Térmico Analysis

autofem.com

Volver al índice