De disco que tiene una temperatura constante en la superficie cilíndrica externa está siendo calentado a 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 un temperatura constante de 273.15 oK se lleva a cabo.

es celebrada. 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 cálculo de los siguientes resultados se obtienen:

Tabla 1. Parámetros de malla de elementos finitos

Tipo de elementos finitos.

Número de Nodos

Número de elementos finitos.

Tetraedro lineal

3389

6592

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

349.9406

349.7972

0.04

R 40

331.5473

331.4828

0.02

R 50

317.2970

317.2771

0.01

R 60

305.6963

305.6701

0.01

 

Conclusiones:

Error relativo de la solución numérica. en comparación con la analítica Solución no exceder de 0,04% cuando se usa 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 numérica.s analíticas y solución prácticamente coinciden. Esto implica que la distribución de los máximos y mínimos temperatura son los mismos, y por lo tanto, el cálculo de los flujos de calor termino y el poder se llevarán a cabo con el mismo grado de precisión.

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