Placa Isotrópica bajo flujo de calor y convección

Consideremos una placa de acero bxa = 100x200 mm de espesor D = 5 mm con la conductividad térmica de K = 50 W / (m oK) [aleación de acero (SS)]. Vamos a mostrar que el problema de la transferencia de calor en estado estacionario también puede resolverse para el caso de que en el límite de la  temperatura  no se mantiene constante (para el modelado de la  temperatura  correcta debe ser conocida de antemano). Vamos a considerar dos tipos de cargas: flujo de calor (igual a cero para algunos límites) de transferencia de calor y por convección. En la superficie de la placa (de ambos lados) prescribimos el intercambio de calor con el entorno ambiental con el coeficiente de transferencia de calor C=40 W/(m2oK) y la  temperatura  del medio ambiente igual a cero. En el cuarto lado (de 100 mm de longitud) prescribimos una transferencia de calor por convección con el coeficiente H=20 W/(m2oK) y la  temperatura  del entorno ambiental es igual a u0=283.15 oK. (ver la figura).

Vamos a localizar el origen del sistema de coordenadas en la esquina de la placa, dejar que el eje O y se dirige a lo largo del lado más largo de la placa, y convección entonces se aplica al borde y = b. A continuación vamos a inferir que dos de cada tres bordes, en los que se prescribe un flujo de calor cero, corresponden a los valores x = 0yx = a. Esto significa que la solución que se obtiene no tiene que depender de la coordenada x. Solución tendrá la forma:

A continuación les mostramos las líneas isotérmica para la solución obtenida mediante el software de arce.

La convección es aplicado para la coordenada Y=b.

líneas isotérmicas para la solución obtenida mediante el software de arce. La convección se aplica al borde y = b. Cuando y = 0 el flujo de calor térmico es cero

Vamos a localizar los sensores de la  temperatura  como se muestra en la figura y crea un complot para ellos en y = 170; 180; 190; 200 mm. En los puntos dados compararemos la Solución numérica. obtenido utilizando el Análisis AutoFEM con la solución analítica.

El modelo de elementos finitos con condiciones de frontera aplicados

Después de llevar a cabo el cálculo se obtienen 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

462

1200

Tabla 2. Resultado " temperatura "

Coordinate y, mm

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

Solución analítica
temperatura T, К

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

170

273.2773

273.2783

3.55E-04

180

273.3029

273.3033

1.43E-04

190

273.3325

273.3332

2.67E-04

200

273.3669

273.3690

7.86E-04

 

Conclusiones:

El error relativo de la solución numérica. en comparación con la solución analítica no exceda 0,0008% para los elementos lineales.

El problema se resolvió casi exactamente con el mínimo coste computacional.

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