AutoFEM Analysis De disco con una temperatura constante en la superficie cilíndrica externa se calienta lo largo del eje | |||||||
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. |
Solución analítica |
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