DEPARTAMENTO DE INGENIERÍA CIVIL Y AMBIENTAL

Diseño Avanzado de Pavimentos – ICYA4602

2024-20

Laura Sofía Méndez - 202014539

María Paula Perilla – 202013117

Tarea 4
1. Introducción
La evaluación de los efectos de las condiciones climáticas sobre la superficie de un pavimento
es crucial porque determina el comportamiento de las capas, especialmente las asfálticas. Esto
se debe a que las propiedades de estos materiales, como la resistencia y durabilidad, se ven
afectadas por la temperatura y la radiación solar. Con el propósito de observar la evolución de
los componentes del flujo de calor y de la temperatura a lo largo del tiempo dentro de la
estructura del pavimento, se ha seleccionado Barrancabermeja como lugar de estudio y el 15 de
junio como día de análisis, bajo condiciones climáticas soleadas y con cielo despejado. La
evaluación se llevará a cabo utilizando datos iniciales de temperatura del aire (Ta), humedad
relativa (Uw) y velocidad del viento (Va) para cada hora del día seleccionado. Se analizarán los
resultados en diferentes momentos del día (6h, 12h, 18h y 24h) para obtener una comprensión
detallada de las variaciones térmicas. Los espesores por capa del caso de análisis se muestran en
la figura 1. Además, se desarrollará un mapa de temperatura que muestre estas variaciones a
diferentes profundidades y horas del día. Esta información es clave para el diseño y
mantenimiento de los pavimentos, asegurando un desempeño óptimo y aumentando su vida útil.

Figura 1. Espesores de las capas de pavimento por analizar

2. Procedimiento/Resultados/Análisis
Para poder solucionar el problema numéricamente se realizó el método de diferencias finitas
teniendo en cuenta los parámetros iniciales y siguiendo la serie de pasos descrita en figura 2.
 Figura 2. Diagrama de flujo

La estructura de pavimentos a analizar tiene 4 capas de las cuales sus características se describen
en la tabla 1. Cabe mencionar que la densidad del material específica para cada capa se obtuvo
según una subbase clase 1 y una base clase 1A (HOLCIM, s. f.), según una capa asfáltica
convencional (Ortiz Hernández et al., 2019) y según un estudio de suelos realizado la ciudad de
Barrancabermeja (TOPOGRAFIA Y LABORATORIO DE SUELOS, 2022).
Tabla 1. Parámetros iniciales de las capas del pavimento

Con térmica Cap

Espesor Cap calorífica Densidad
Capa KH calorífica CH
(m) CH (J/(m3K)) material (kg/m3)
(W/(m°C)) (J/(kgK))
Asfáltica 0.2 1.500 2.00E+06 2375 1000*0.842
Base Granular 0.2 2.159 2.58E+06 2382 1000*1.082
Subbase
0.3 1.793 2.71E+06 2365 1000*1.144
Granular
Subrasante Infinita 1.506 3.05E+06 1830 1000*1.668
Por otro lado, la emisividad usada va a ser igual a la absortividad que tenia por entrada un valor
de 0.85 tal como lo muestra la ecuación 1. Es importante mencionar que la constante de Stefan-
Boltzmann (𝜎) es igual a 5.67E-8 W/(m2 K4) y la constante solar es igual a 1368 W/m2 .
𝛼 = 𝜖 = 0.85
Ecuación 1. Absortividad y emisividad

Las variables ambientales de temperatura del aire (Ta), humedad relativa (Uw) y velocidad del
viento (Va) usadas para el ejemplo por cada hora del día de análisis (15 de junio) se muestran en
la tabla 2 (Weather Spark, 2023). De cada serie de datos se muestra el promedio, la desviación
estándar y su coeficiente de radiación. Es importante mencionar que para el método se partió de
una temperatura inicial para la ciudad de 23°C. Lo anterior dado que es una temperatura cercana
a la que se tiene en la ubicación en el día de análisis y se desea que el método converja rápidamente
a las temperaturas del aire presentadas. Adicionalmente, la ciudad de Barrancabermeja tiene una
latitud de 7° 3’ 55’’ (Date and Time, 2024).
Tabla 2. Parámetros atmosféricos iniciales

Ta Uw Va Va
Hora (°C) (%) (km/h) (m/s)
1 25.50 88 2.15 0.597
2 25.00 88 2.20 0.611
3 24.50 90 2.25 0.625
4 24.50 91 2.30 0.639
5 24.50 92 2.34 0.650
6 24.50 92 2.20 0.611
7 24.00 89 2.00 0.556
8 25.00 80 2.05 0.569
9 26.00 75 2.50 0.694
10 27.00 70 2.90 0.806
11 27.50 64 3.40 0.944
12 28.00 60 3.50 0.972
13 28.00 58 3.60 1.000
14 28.50 55 3.60 1.000
15 29.00 54 3.20 0.889
16 28.50 55 2.60 0.722
17 28.50 61 1.90 0.528
18 28.00 70 1.70 0.472
19 27.50 76 2.20 0.611
20 27.00 78 2.10 0.583
21 26.50 81 2.05 0.569
22 26.00 82 2.00 0.556
23 25.50 85 2.00 0.556
24 27.00 86 2.00 0.556
Promedio 26.50 75.83 2.45 0.68
 Desv 1.57 13.25 0.59 0.16
COV 5.94% 17.48% 23.96% 23.96%

Teniendo en cuenta que el día a evaluar es el 15 de junio, el intervalo escogido de análisis para
que la temperatura converja con el paso de los días es desde el día 150 al 166 del año, siendo
este último el 15 de junio. Ahora bien, para la discretización en el tiempo y espacio se usa un
delta de tiempo de 120 segundos y un delta de profundidad de 0.02 metros.

A partir de los parámetros de entrada descritos anteriormente, se proceden a calcular los

parámetros ambientales del ejemplo a lo largo del tiempo. Las ecuaciones 2 a 5 muestran cómo
se calcularon los parámetros ambientales teniendo en cuenta las variables de entrada siendo 𝑇0
la temperatura absoluta de la superficie en grados Kelvin.

𝑚𝑇𝑎 17.62𝑇𝑎
𝑇𝑛 [ + ln(𝑈𝑤 )] 243.12 [ + ln(𝑈w )]
𝑇𝑛 + 𝑇𝑎 243.12 + 𝑇𝑎
𝑇𝑑 = =
𝑚𝑇 17.62𝑇𝑎
𝑚 − [𝑇 + 𝑎𝑇 + ln(𝑈𝑤 )] 17.62 − [ )]
𝑛 𝑎 243.12 + 𝑇𝑎 + ln(𝑈w
Ecuación 2. Temperatura de rocío

𝑇sky = (𝑇𝑎 + 273.15)(0.77 + 0.0038𝑇𝑑 )0.25 [K]

Ecuación 3. Temperatura atmosférica superficial

ℎ𝑐 = 5.8 + 4.1𝑉𝑎 [W/m2 K]

Ecuación 4. Coeficiente de convección

2
ℎrad = ϵσ(𝑇0 + 𝑇sky )(𝑇02 + 𝑇sky ) [W/m2 K]
Ecuación 5. Coeficiente de radiación

De esta manera, los parámetros atmosféricos a lo largo del día se presentan en la siguiente figura.
Es relevante mencionar que para la humedad relativa y la velocidad del viento a lo largo del día
se usó el promedio para la discretización en términos de tiempo y espacio.
 Figura 3. Evolución de las variables atmosféricas a lo largo del día

Respecto a los resultados obtenidos de la evolución de las variables atmosféricas a lo largo del
día, es posible observar que la temperatura del aire alcanza un valor máximo de
aproximadamente 29°C alrededor de las 3 p.m., mientras que la humedad relativa alcanza un
valor mínimo de aproximadamente 55% en ese punto. Por otro lado, la velocidad del viento varía
entre 2 y 3.7 km/h, alcanzando el mínimo alrededor de las 6 p.m. y el máximo alrededor de la 1
p.m. Según estos datos iniciales, se puede evidenciar que la temperatura del rocío varía entre 18
y 24°C, y la temperatura atmosférica superficial varía entre 276 y 285 K. Para estas dos variables
anteriores, se alcanza el punto máximo a primeras horas de la mañana y el mínimo alrededor de
las 2 p.m. Finalmente, para los coeficientes de convección y radiación se observa una variación
de entre 14 y 20 W/m2 K para el primero y un valor constante (alrededor de 4 W/m2 K) para el
segundo. Cabe mencionar que el valor máximo del coeficiente de radiación se alcanza alrededor
de la 1 p.m.
Para encontrar los valores de temperatura a lo largo del asfalto durante cierto periodo de tiempo
dado se utiliza la solución explicita que permite discretizar el problema en deltas de tiempo y de
profundidad. Dependiendo de la ubicación del nodo a evaluar se usa ciertas ecuaciones, pues las
condiciones de borde afectan en el análisis, teniendo en cuenta los parámetros de conductividad
térmica y capacidad calórica definidos en la tabla 1. Es importante mencionar que todos los
parámetros en el tiempo “t” descritos en las siguientes ecuaciones son conocidos, pues se
estableció previamente las condiciones ambientales en cada hora del día, el delta de tiempo
evaluado fue dos minutos. Para la distribución de calor en la capa superficial del pavimento,
punto amarillo mostrado en la figura 3, se deben tener en cuenta que la temperatura es controlada
por las condiciones climáticas.
2Δ𝑡
𝑇1𝑡+Δ𝑡 = [Q ] + 𝑇1𝑡
𝑐𝐻 Δ𝑧 Hsurf
Ecuación 6. Temperatura en tiempo t+dt para nodo superior (Caicedo, 2021)

Donde el flujo de calor de la superficie es

4 𝑇𝑖𝑡 − 𝑇𝑖
𝑄𝐻surf = α𝐼 + ℎ𝑐𝑡 (𝑇𝑎𝑡 − 𝑇𝑖𝑡 ) + ϵσ [(𝑇𝑠𝑘𝑦
𝑡
) − (𝑇𝑖𝑡 )4 ] − 𝑘𝐻
Δ𝑧
Ecuación 7. Flujo de calor para el nodo superior (Caicedo, 2021)

Figura 4. Ubicación nodo superior

Partiendo de este dato de temperatura inicial longitudinalmente es posible calcular la

temperatura entre nodos internos, punto amarillo mostrado en la figura de a continuación, entre
cada capa por medio de la siguiente ecuación
𝑘𝐻 Δ𝑡
𝑇𝑖𝑡+Δ𝑡 = (𝑇 𝑡 + 𝑇𝑖+1
𝑡
− 2𝑇𝑖𝑡 ) + 𝑇𝑖𝑡
𝑐𝐻𝑣 Δ𝑧 2 𝑖−1
Ecuación 8. Temperatura en tiempo t+dt para nodo entre capas (Caicedo, 2021)

Figura 5. Ubicación nodo entre capas

Es importante destacar que, para los límites entre capas, punto amarillo mostrado en la figura
siguiente, esta ecuación varía teniendo en cuenta que 𝑘𝐻 y 𝐶𝐻𝑣 son diferentes entre capa.

𝑇𝑖𝑡+Δ𝑡 = 𝐴int 𝑇𝑖−1

𝑡
+ 𝐵int 𝑇𝑖𝑡 + 𝐶int 𝑇𝑖+1
𝑡

Ecuación 9. Temperatura en tiempo t+dt para nodo en límite entre capas (Caicedo, 2021)

Los valores A, B y C son los coeficientes por multiplicar por la temperatura entre los nodos
vecinos.
𝑖−1/2
2Δ𝑡 𝑘𝐻
𝐴int = 2 ( 𝑖−1/2 )
Δ𝑧 𝑐 + 𝑐𝐻𝑖
𝐻

Ecuación 10. Coeficiente A para encontrar la temperatura del nodo límite

𝑖−1/2 𝑖+1/2
2Δ𝑡 𝑘𝐻 𝑘𝐻
𝐵int = 1 − 2 ( 𝑖−1/2 + 𝑖+1/2 ) = 1 − (𝐴int + 𝐶int )
Δ𝑧 𝑐 + 𝑐𝐻𝑖 𝑐 + 𝑐𝐻𝑖
𝐻 𝐻

Ecuación 11. Coeficiente B para encontrar la temperatura del nodo límite

Note que B, el cual corresponde al nodo justo que se encuentra en el límite de las capas, es una
especie de promedio entre los valores para el nodo que lo antecede y el que lo precede.
𝑖+1/2
2Δ𝑡 𝑘𝐻
𝐶int = 2 ( 𝑖+1/2 )
Δ𝑧 𝑐 + 𝑐𝐻𝑖
𝐻

Ecuación 12. Coeficiente C para encontrar la temperatura del nodo límite

Figura 6. Ubicación nodo límite

Finalmente, para el último nodo evaluado la ecuación varía de la siguiente manera.

𝑘𝐻 Δ𝑡 𝑡
𝑇𝑛𝑡+Δ𝑡 = (2𝑇𝑛−1 − 2𝑇𝑛𝑡 ) + 𝑇𝑛𝑡
𝑐𝐻 Δ𝑧 2
Ecuación 13 Temperatura en tiempo t+dt para último nodo (Caicedo, 2021)
 Figura 7. Ubicación nodo inferior

Partiendo del procedimiento descrito anteriormente, se evalúa el flujo de calor en la superficie

y para ello, se involucraron los componentes que lo describen. Estos son el calor de convección,
radiación, emisión térmica y el flujo de calor entrando a la estructura de pavimento. En la figura
7 se exponen los valores de estas variables a lo largo del día. Se puede notar que, en todos los
casos los valores máximos de calor se alcanzan sobre 11 a.m. y la 1 p.m. del día. Además, el
pico más alto se encuentra en el calor de radiación con un valor de aproximadamente 800, para
el calor de convección el pico se encuentra en 270 y para los demás se alcanzan los máximos en
350. Así, se debe tener en cuenta que el módulo cambia con la temperatura dependiendo de la
radiación.

Figura 8. Evolución a lo largo del día de los componentes del flujo de calor

Sumado a lo anterior, en la siguiente figura se muestra la variación de la temperatura a lo largo del

día para la superficie de cada una de las capas del pavimento. Tal como se puede observar, la
superficie (capa asfáltica) alcanza las temperaturas más altas (de alrededor de 48°C) mientras que
la subrasante conserva las temperaturas más bajas. Igualmente, se puede evidenciar que entre más
cerca de la superficie de la vía se encuentre la capa, se enfrentará a mayores temperaturas. Cabe
recalcar que un factor que influencia la sensación de temperaturas más extremas en primera capa
es que los materiales asfalticos son más sensibles a los cambios de temperatura y, por ende, a fallas
en condiciones extremas.

Figura 9. Temperatura de las capas de pavimento a lo largo del día

Por otro lado, la variación de la temperatura se describe en la figura de a continuación en ciertas

horas particulares (6h, 12h, 18h y 24h) a lo largo de la profundidad. Se puede evidenciar que las
mayores temperaturas se encuentran a menor profundidad. Adicionalmente, se nota que las
menores temperaturas se tienen a las 6 y 24 horas mientras que las mayores son a las 12 y 18 horas.
Lo anterior sucede debido a que es después de las 6 primeras horas donde la ciudad alcanza las
mayores temperaturas que se ven reflejadas tiempo después en el pavimento. Así, las siguientes
horas después de alcanzar el máximo, las temperaturas del ambiente bajan y, por ende, tiempo
después las de la estructura.
 Figura 10. Temperatura en las horas de análisis a lo largo de la profundidad

Finalmente, se muestra en la siguiente figura el mapa de calor por toda la superficie del pavimento
a lo largo del día de análisis (15 de junio) después de que la temperatura haya convergido después
de los días seleccionados. Como se puede notar, las mayores temperaturas se encuentran más
cercanas a la superficie entre las 12 p.m. y 3 p.m.
 Figura 11. Mapa de calor para el día de análisis

3. Conclusiones
Es de gran importancia determinar correctamente los datos ambientales, pues en un proyecto
grande no sólo se requerirían valores de temperatura del ambiente, humedad relativa y velocidad
del viento para un día normal, por el contrario, se debería considerar el día con las condiciones
climáticas más críticas con el fin de estudiar los efectos sobre el asfalto y el pavimento en
general.
En el modelo se utilizó una humedad y velocidad del viento promedio en el día evaluado debido
a que en Barrancabermeja no se presentó una variación significativa en el rango de estos datos,
sin embargo, en caso de que varíen considerablemente se deben generar nuevos ciclos
evaluativos en los cuales se establezcan las condiciones ambientales para cada hora.
Para lograr un diseño óptimo de mezcla asfáltica en Barracabermeja se deberían tener en cuenta
las temperaturas más altas del ambiente, que se dan aproximadamente entre 10 y 12 pm, pero
que se ven reflejadas en la estructura hasta casi las 3 pm.
En general, se pudo observar que el material granular disipa en gran medida la temperatura
absorbida por el asfalto, pues en la figura de temperatura según profundidad se puede observar
que de casi 50°C que tiene la superficie a los 2 metros de profundidad logra alcanzar a ser casi
la mitad.
4. Referencias
Caicedo, B. (2021). Geotechnics of Roads: Advanced Analysis and Modeling.
https://www.amazon.com/Geotechnics-Roads-Advanced-Analysis-
Modeling/dp/113860058X
Date and Time. (2024). Las coordenadas geográficas de Barrancabermeja. La latitud, la
longitud y la altitud sobre el nivel del mar de Barrancabermeja, Colombia.
dateandtime.info.
https://dateandtime.info/es/citycoordinates.php?id=3689169#google_vignette
HOLCIM. (s. f.). AGREGADOS Material granular para vías. En HOLCIM.
Ortiz Hernández, E., Macías, L., & Ortiz Moncayo, E. (2019). COMPARATIVO DE LAS
PROPIEDADES DE UN DISEÑO DE MEZCLA ASFÁLTICA EN CALIENTE
CONVENCIONAL y EL USO DE POLÍMEROS EN LA CARRETERA TOSAGUA.
UNIVERSIDAD, CIENCIA y TECNOLOGÍA.
TOPOGRAFIA Y LABORATORIO DE SUELOS. (2022). Informe de Estudio de Suelos para
la Construcción de Vías Urbanas. En INVERSIONES & SOLUCIONES HD (NIT:
1120740031-1).
Weather Spark. (2023). El tiempo en Barrancabermeja en junio, temperatura promedio
(Colombia) - Weather Spark. https://es.weatherspark.com/m/24382/6/Tiempo-
promedio-en-junio-en-Barrancabermeja-Colombia