SciELO - Scientific Electronic Library Online

 
vol.14 issue2OPTIMIZATION OF COVERAGE FOR GEOREFERENCED ZONESHYBRID BSS TECHNIQUES FOR FOETAL ECG EXTRACTION USING FRAMEWORK FOR STRESS-TESTING EXTRACTION ALGORITHMS author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • On index processCited by Google
  • Have no similar articlesSimilars in SciELO
  • On index processSimilars in Google

Share


Iteckne

Print version ISSN 1692-1798

Iteckne vol.14 no.2 Bucaramanga July/Dec. 2017

https://doi.org/10.15332/iteckne.v14i2.1769 

Artículos de investigación e innovación

APLICABILIDAD DEL MÉTODO DE LOS ELEMENTOS FINITOS EN EL ANÁLISIS Y DIMENSIONAMIENTO DE LOSAS JCPC PARA CARRETERAS DE DOS CARRILES

APPLICABILITY OF THE FINITE ELEMENT METHOD IN ANALYSIS AND SIZING OF JCPC SLABS FOR TWO-LANE ROADS

Myriam Rocío Pallares-Muñoz1 

Julián Andrés Pulecio-Díaz2 

1M. Sc. Métodos Numéricos en Ingeniería Universidad Surcolombiana. Neiva, Colombia. myriam.pallares@usco.edu.co.

2M. Sc. Construcción Obras Viales Universidad Surcolombiana. Neiva, Colombia. julian.pulecio@usco.edu.co.


Resumen

Se analizaron y dimensionaron losas planas de concreto con juntas (Jointed-Plain-Concrete-Pavement “JCPC” en inglés) para carreteras de dos carriles por medio del modelador de elementos finitos EverFE 2.25, intentando reproducir de una manera más ajustada a la realidad los esfuerzos del conjunto de losas de pavimento. Para la validación se empleó el método analítico Westergaard 1948. Los resultados de esfuerzos y número de repeticiones de carga permitidos llevaron a concluir que el diseño de losas JCPC para carreteras de dos carriles no requiere el peso propio de la estructura en el análisis y dimensionamiento porque el porcentaje de variación obtenido es poco significativo. Por lo tanto, se recomienda omitir el parámetro densidad en el análisis para generar un coste computacional menor. Este estudio traza una ruta para seguir analizando otros factores que puedan afectar el dimensionamiento de pavimentos rígidos. Se espera en próximas investigaciones integrar cargas de distintos vehículos y los comportamientos producidos por la presencia de dovelas y barras de amarre en juntas transversales y longitudinales.

Palabras clave: Carreteras de dos carriles; EverFE 2.25; losas JCPC; modelado tridimensional de pavimentos; modelado de pavimentos con software libre; Westergaard 1948.

Abstract

Two-lane road JCPC (Jointed-Plain-Concrete-Pavement) slabs were analyzed and dimensioned using the EverFE 2.25 finite element modeler in order to reproduce the efforts of the pavement slab in a more realistic way. The “Westergaard”-1948 analytical method was used for the validation. With the results of efforts and the number of allowed repetitions of load we can conclude that the designs of JCPC slabs for two lanes do not require the analysis and dimensioning with own weight of the structure, because the percentage of variation is little significant. Therefore, we recommend omitting the density parameter in the analysis to generate a lower computational cost. This study shows us a route to continue analyzing other factors that influence the sizing of rigid pavements. We hope in future investigations to integrate loads of different vehicles and the behaviors produced by the presence of segments and tie bars in transverse and longitudinal joints.

Keywords: Two-lane roads; EverFE 2.25; Slabs JCPC; three-dimensional modeling of pavements; freeware pavement modeling; Westergaard 1948.

1. INTRODUCCIÓN

En la actualidad para el dimensionamiento de estructuras de pavimentos rígidos para carreteras de dos carriles se emplean métodos empíricos y empíricos-mecanísticos, definiendo comúnmente para ello: el eje de diseño (eje simple de 8.2T (80kN)) o el espectro de carga, el espesor de las capas y las características de los materiales. De manera que se pueda calcular los máximos esfuerzos principales a tensión en las fibras inferiores o superiores de la losa (fatiga)1-3 y deformaciones verticales en las capas granulares y subrasante (deformación permanente). Orientado a estudios en una a dos dimensiones que conllevan omitir condiciones de frontera más amplias. Este hecho genera la necesidad de realizar un análisis y dimensionamiento mucho más recíproco con la realidad, teniendo en cuenta un modelado en tercera dimensión, con y sin peso propio de la estructura, carga más agresiva y condiciones de liga entre losas y demás capas (base y subrasante). Para este propósito se empleó el modelador tridimensional de elementos finitos EverFE 2.25 y el método analítico Westergaard 1948 para validar los resultados.

2. MÉTODOS

El análisis se abordó por medio del modelador computacional de elementos finitos EverFE 2.25, validando sus resultados con el método analítico Westergaard 1948. La estructura de pavimento rígido está compuesta por una base granular y una capa de rodadura de losas JCPC (Jointed Plain Concrete Pavement) como se muestra en la Fig. 1.

Fuente: los autores.

Figura 1 ESTRUCTURA PAVIMENTO RÍGIDO DE LOSAS JCPC 

Las dimensiones de las capas de la estructura de pavimento modelado se presentan en la Tabla I. La Tabla I presenta los espesores en milímetros (mm) y pulgadas (in), ya que el modelador de elementos finitos EverFE 2.25 y el método analítico de Westergaard 1948 requieren de este tipo de unidades. La condición de frontera considerada para la subrasante es infinita y como es común en el dimensionamiento de este tipo de estructuras se trabajaron las interfaces ligada entre subrasante y base granular y no ligada entre base granular y losas. Las losas son cuadradas de 3.5m*3.5m y las aberturas de las juntas de 5 mm de espesor, ambos son parámetros geométricos que se anticipan a las fallas por contracción y alabeo del concreto hidráulico. Las fallas por contracción ocurren por secado y enfriado del material. Se dan usualmente entre las primeras 12 o 24 horas de ser colocado y colado el concreto hidráulico de cemento tipo I. Las fallas por alabeo se dan en el tiempo Jeong4.

Tabla I VALORES TÍPICOS DE ESPESORES DE CAPA DE UNA ESTRUCTURA DE PAVIMENTO RÍGIDO DE LOSAS JCPC 

Capa Espesor (mm) Espesor (in)
Losas JCPC 228.6 9
Base granular 150 5.905
Subrasante infinita infinita

* sin corregir

Fuente: los autores.

Los módulos de rigidez tomados como referencia para cada una de las capas de acuerdo con las condiciones de trabajo del pavimento se muestran en la Tabla II. El módulo de reacción de la subrasante con denominación en inglés dense liquid subgrade, se determina por medio de un ensayo de placa de carga (PLT) corregido en ocasiones por: a) saturación del suelo cuando es cohesivo, b) asentamiento deficiente de la placa, comportamiento carga deformación no lineal o falla por corte, presentes al inicio del ensayo con un comportamiento variable en la gráfica carga versus deflexión, c) flexión de la placa soportante (cuando las deflexiones registradas mediante LVDTs o diales ubicados perimetralmente en la placa son menores a las del centro) Mallick y El-Korchi5-7. El módulo elástico de la base granular conocido como resiliente se determina por medio de ensayos triaxiales cíclicos Ping et al.8,9, que para el caso es un valor característico del módulo de reacción de la subrasante y el espesor de base granular empleados. El módulo elástico de las losas de concreto hidráulico se establece a través de ensayos estáticos a carga uniaxial Joh et al.10,11.

Tabla II VALORES TÍPICOS DE MÓDULOS DE RIGIDEZ DE LAS CAPAS 

Capa Módulo elástico (MPa) Módulo elástico (lb/pul 2 ) Módulo reacción (MPa/mm) Módulo reacción (lb/pul 3 )
Losas JCPC 28000 4000000 - -
Base granular 82.48 11782.85 - -
Subrasante - - 0.04082* 150*

* sin corregir

Fuente: los autores.

Ya que la variación del módulo de Poisson no tiene incidencia considerable en el comportamiento del pavimento Huang12, se emplearon los valores característicos de la tabla III.

Tabla III VALORES TÍPICOS DE MÓDULO DE POISSON DE LAS CAPAS 

Capa Módulo poisson
Losas JCPC 0.15
Base granular 0.40
Subrasante -

Fuente: los autores.

La determinación de los parámetros de diseño de pavimentos rígidos tiene principio en el diseño empirico-mecanístico, que consiste en modelar la estructura a partir de la definición de los espesores y las rigideces (módulos elásticos, reacción y de Poisson) de cada capa para calcular las tensiones y compresiones provocadas por una carga tipo e identificar: a) el máximo esfuerzo principal a tensión en las fibras inferiores o superiores de la losa (σt) capaz de producir rotura en las losas; este esfuerzo se compara con el límite admisible del concreto hidráulico (σtadm) que es función del tránsito, esfuerzo a tensión aplicado y módulo de rotura, siendo este último parámetro determinado mediante el ensayo de viga de cuatro puntos; b) máximas deformaciones a compresión (εc) en las capas granulares y subrasante que se compara con la admisible (εcadm) y que depende del nivel de tránsito LCPC13.

2.1 Método analítico Westergaard 1948

Se calcularon los máximos esfuerzos principales a tensión en el interior, esquina y borde de la losa, a partir de las expresiones (1), (2), (3), (4) y (5) de Westergaard 1948, introduciendo una presión de inflado distribuida uniformemente sobre un área circular de radio 0.2286 m (9 in) bajo una solicitación de carga de 40.1 kN (9000 lb), que representa la mitad de un eje simple de rueda sencilla. El análisis permite identificar los esfuerzos característicos de control del diseño estructural del pavimento rígido, tales como los esfuerzos a tensión en las fibras inferiores o superiores de la losa (fatiga) Papagiannakis y Masad14 Reddy15.

(1)

(2)

(3)

(4)

(5)

Donde E es el módulo elástico del concreto hidráulico, d espesor de la losa de concreto hidráulico, μ módulo de poisson del concreto hidráulico, K módulo de reacción de la subrasante, P solicitación de carga, a radio del área circular, radio de rigidez relativa, σI máximo esfuerzo principal a tensión en el interior de la losa, σE máximo esfuerzo principal a tensión en la esquina de la losa y σB máximo esfuerzo principal a tensión en el borde de la losa.

2.2 Modelos de elementos finitos EverFE 2.25

Los modelos de elementos finitos con el modelador tridimensional elástico EverFE 2.25 fueron validados con el método analítico Westergaard 1948. Las Figs. 2, 3, 4 y 5 muestran los elementos finitos que emplea el software Davids et al.16. El EverFE 2.25 trabaja con mallados de distintos grados de refinamiento y se compone de seis interfaces que realizan diferentes funciones: a) Geometry, contiene una herramienta CAD que permite construir un modelo geométrico e introducir los datos característicos de las capas y definir las fronteras del dominio del modelo. b) Material, almacena los parámetros de rigidez de cada material, densidades, coeficiente de expansión térmica del concreto hidráulico y las condiciones de liga entre capas. c) Loading, almacena los parámetros de carga, tipo de eje, lado y ancho del contacto de rueda (área de contacto rectangular), ubicación sobre el plano x-y de la(s) losa(s) y carga por temperatura. d) Dowel, contiene los parámetros de rigidez y geometría de las dovelas y barras de amarre, además, condiciones de liga entre dovelas y el concreto hidráulico de la(s) losa(s). e) Interlock, almacena abertura de las juntas con su respectiva rigidez, acompañado de la opción de modelo lineal a no lineal. f) Meshing, determina los parámetros de mallado. g) Solver, ejecuta el solucionador de elementos finitos. h) Results, permite obtener los resultados del análisis. Los pasos típicos del procedimiento de análisis de un problema con EverFE 2.25, son: ingresar los datos de dimensionamiento, número de capas, condiciones de frontera, espesores y parámetros de rigidez (módulos elásticos, reacción y de Poisson), otros factores inherentes al material (densidades y coeficiente de expansión térmica del concreto hidráulico), carga, ancho y largo del área de contacto (el programa trabaja con área de contacto rectangular) espaciamiento de las llantas, carga por temperatura, dovelas y barras de amarre con sus propiedades geométricas, de rigidez y liga con el concreto hidráulico, rigidez y abertura de las juntas, elegir combinación de mallado, ejecutar el solucionador y obtener resultados Davids17-19.

Fuente: 15.

Figura 2 ELEMENTO FINITO BRICK CON 20 NODOS 

Fuente: 15.

Figura 3 ELEMENTO FINITO INTERFACE CON 16 NODOS 

Fuente: 15.

Figura 4 ELEMENTO FINITO DENSE LIQUID CON 8 NODOS 

Fuente: 15.

Figura 5 ELEMENTO FINITO BEAM CON 2 NODOS 

2.2.1 Área de contacto

Como el método analítico Westergaard 1948 maneja en su formulación un área de contacto circular, diferente al presentado por el modelador tridimensional elástico EverFE 2.25, se calcula en las expresiones (6) y (7) la relación entre las áreas de contacto circular y rectangular. La Fig. 6 justifica la procedencia del nuevo dato de cálculo. Un estudio de los efectos de otros tipos de distribuciones de carga se trata en la referencia de los autores20.

(6)

(7)

Donde L es el ancho del área de contacto rectangular, A el área de contacto y a el radio del área circular.

Fuente: los autores.

Figura 6 ÁREAS DE CONTACTO CIRCULAR Y RECTANGULAR 

2.3 Funciones de transferencia PCA

Se aplicaron las funciones de transferencia de la Portland Cement Association - PCA Packard y Tayabji21, calculando el número de repeticiones de carga permitidas con las expresiones (8), (9) y (10).

(8)

(9)

(10)

Donde Nf número de repeticiones de carga permitidos, σ esfuerzo de tensión aplicado, Sc módulo de rotura del concreto y σ/Sc Relación de esfuerzo.

3. RESULTADOS

En la Tabla IV se presentan los resultados del método analítico Westergaard 1948.

Tabla IV RESULTADOS DE CONTROL DE DISEÑO ESTRUCTURAL POR FATIGA DE UN PAVIMENTO RÍGIDO DE LOSAS JCPC 

Ubicación carga Parámetro Analítico
Interior σI (Mpa) 1.05
Esquina (Corner) σE (Mpa) 1.35
Borde (Edge) σB (Mpa) 1.99

Fuente: los autores.

Se observa que el máximo esfuerzo principal a tensión está ubicado en el borde de la losa. Lo cual conlleva tomar la decisión de dimensionar las losas JCPC para pavimentos con solicitación de carga en el borde.

La tabla V contiene los resultados con solicitación de carga en el borde del método analítico Westergaard 1948 - utilizado como validador - y el modelo de elementos finitos EverFE 2.25.

Tabla V VALIDACIÓN DE LOS RESULTADOS DE CONTROL DE DISEÑO ESTRUCTURAL DE UN PAVIMENTO RÍGIDO DE LOSAS JCPC 

Ubicación carga Parámetro Analítico Elementos Finitos 3D %Variación
Borde (Edge) σB (Mpa) 1.99 1.88 5.5 %

Fuente: los autores.

Se observa que los porcentajes de variación del modelo numérico frente al analítico son inferiores al 5.5% lo cual permite concluir que las decisiones tomadas en la modelación con EverFE 2.25 son adecuadas para implementar los modelos que introduce base granular y un modelo más robusto de losas JCPC. En la Fig. 7 se presenta el resultado del modelo de elementos finitos EverFE 2.25.

Fuente: los autores.

Figura 7 MÁXIMOS ESFUERZOS PRINCIPALES EN EL PLANO X-Y 

La Fig. 7 muestra el resultado gráfico de esfuerzos, correspondiente al parámetro de diseño empírico-mecanístico por medio de isocontornos. En la Fig. 7 se observa que los máximos esfuerzos se concentran bajo la llanta y en las fibras inferiores de la losa tal como lo sugiere el método analítico Westergaard 1948.

3.1 Modelo de losas JCPC para carreteras de dos carriles

Se desarrollaron dos modelos para analizar la respuesta de esfuerzos de la estructura de pavimento, conformado por tres losas cuadradas por carril en serie, siendo cada una de 3.5m*3.5m de dimensión en planta y abertura de juntas de 5 mm de espesor, dimensiones características para anticipar las fallas por contracción y alabeo del concreto hidráulico. Se evaluó con y sin peso propio la estructura, para contemplar factores que lleven a tener análisis más realistas.

3.1.1 Sin peso propio de la estructura

La Fig. 8 presenta el desarrollo del esfuerzo en el plano X-Y de las fibras inferiores de las losas JCPC para pavimento de carretera de dos carriles sin peso propio. Bajo la llanta se concentra el máximo esfuerzo principal σB = 1.89Mpa.

Fuente: los autores.

Figura 8 ESFUERZOS EN EL PLANO X-Y DEL PAVIMENTO DE CARRETERA DE DOS CARRILES SIN PESO PROPIO 

3.1.2 Con peso propio de la estructura

La Fig. 9 presenta el desarrollo del esfuerzo en el plano X-Y de las fibras inferiores de las losas JCPC para pavimento de carretera de dos carriles con peso propio. Peso que se activa mediante el parámetro gravimétrico densidad, siendo para el concreto hidráulico 2200 kg/m3 y la base granular 2100 kg/m3. Bajo la llanta se concentra el máximo esfuerzo principal σB = 1.86Mpa.

Fuente: los autores.

Figura 9 ESFUERZOS EN EL PLANO X-Y DEL PAVIMENTO DE CARRETERA DE DOS CARRILES CON PESO PROPIO 

3.1.3 Resumen de resultados con y sin peso propio de la estructura

La Tabla VI contiene el resumen de resultados con y sin peso propio de la estructura, determinando porcentajes de variación poco significativos.

Tabla VI VARIACIÓN DE LOS RESULTADOS DE CONTROL DE DISEÑO ESTRUCTURAL DE UN PAVIMENTO RÍGIDO DE LOSAS JCPC CON Y SIN PESO PROPIO DE LA ESTRUCTURA 

Ubicación carga Parámetro Sin peso propio Con peso propio %Variación
Borde (Edge) σB (Mpa) 1.89 1.86 1.6 %

Fuente: los autores.

3.1.4 Número de repeticiones de carga permitidos

Para la estructura de pavimento rígido de losas JCPC aplicado para carreteras de dos carriles estudiado hasta el momento, podría soportar 25.253,697 repeticiones de carga, valor calculado con un esfuerzo de tensión de 1.89 Mpa, módulo de rotura de 4 Mpa y relación de esfuerzo de 0.46.

4. CONCLUSIONES

Los resultados obtenidos del modelo de elementos finitos desarrollado con el software libre Ever-FE 2.25 fueron satisfactoriamente validados con el método analítico Westergaard 1948, lo cual permite analizar y dimensionar losas JCPC para carreteras de dos carriles.

Se evidenció que no es necesario contemplar el análisis y dimensionamiento con peso propio de la estructura, porque el porcentaje de variación es poco significativo. Por lo tanto, se recomienda omitir el parámetro densidad sobre el análisis y así generar un coste computacional menor.

El método de los elementos finitos ha demostrado ser muy eficaz y en el caso particular, permitió analizar y dimensionar losas JCPC para carreteras de dos carriles, por lo tanto, marca un camino por seguir para estudiar la sensibilidad de distintos factores que pueden afectar el diseño de pavimentos rígidos.

Generado el interés de profundizar en la aplicabilidad del método de los elementos finitos en el análisis y dimensionamiento de losas JCPC para carreteras de dos carriles, en próximos estudios se analizarán condiciones de carga de distintos vehículos y los comportamientos producidos por la presencia de dovelas y barras de amarre en juntas transversales y longitudinales. De tal manera que se pueda generar análisis y dimensionamiento más realista.

REFERENCIAS

1. A. M. Ioannides, Pavement Fatigue Concepts: A Historical View. Proc., 6th International Conference on Concrete Pavement Design and Materials for High Performance, vol. 3, Purdue University, West Lafayette, Ind., 1997. [ Links ]

2. P. A. Okamoto, “Report on Review of Concrete Fatigue Models”, PCA R&D Serial No. 2213. Portland Cement Association, Skokie, Ill, 1999. [ Links ]

3. L. Titus-Glover, J. Mallela, M. Darter, G. Voigt and S. Waalkes, Enhanced Portland Cement Concrete Fatigue Model for StreetPave. Transportation Research Record: Journal of the Transportation Research Board, vol. 1919, pp. 29-37, Jan, 2005. [ Links ]

4. J.H. Jeong and D. Zollinger, Development of Test Methodology and Model for Evaluation of Curing Effectiveness in Concrete Pavement Construction. Transportation Research Record: Journal of the Transportation Research Board, vol. 1861, pp. 17-25, 2003. [ Links ]

5. R. Mallick and El-Korchi, Pavement engineering: principles and practice. United States of American, CRC Press Taylor & Francis Group., 2013, pp. 48-49. [ Links ]

6. R. Mallick and El-Korchi. Pavement engineering: principles and practice. United States of American, CRC Press Taylor & Francis Group., 2009, pp. 40-41. [ Links ]

7. ASTM, Standard Test Method for Nonrepetitive Static Plate Load Tests of Soils and Flexible Pavement Components, for Use in Evaluation and Design of Airport and Highway Pavements, American Society for Testing and Materials, United States of American, 2016. [ Links ]

8. Ping, Z. Yang, and Z. Gao, “Field and Laboratory Determination of Granular Subgrade Moduli”. ASCE Journal of Performance of Constructed Facilities, vol. 16, No. 4, pp. 149-159, 2002. [ Links ]

9. AASHTO, Standard Method of Test for Determining the Resilient Modulus of Soils and Aggregate Materials. American Association of State Highway and Transportation Officials, United States of American, 2012. [ Links ]

10. S. H. Joh, M. R. Cho, T. H. Kang, S. Kwon and J. Nam, “Pavement Integrity Scanner to Characterize Modulus Contrast Between Near-Surface and Deeper Material in Concrete Pavements”. Transportation Research Record: Journal of the Transportation Research Board , vol. 2070, pp. 39-48, 2008. [ Links ]

11. ASTM, Standard Test Method for Static Modulus of Elasticity and Poisson’s Ratio of Concrete in Compression. American Society for Testing and Materials, United States of American, 2014. [ Links ]

12. Y. Huang, Pavement analysis and desing. Kentucky United States of American, Prentice hall, 2004. [ Links ]

13. LCPC, Conception et dimensionnement des structures de chaussée, Guide technique. LCPC and SETRA, Paris, 1994, pp. 163. [ Links ]

14. A. Papagiannakis and E. Masad, Pavement design and materials. United States of American, John Wiley & Sons, 2008, pp. 221-225. [ Links ]

15. J. N. Reddy, Mechanics of Laminated Composite Plates and Shells Theory and Analysis, 2nd ed. CRC Press, Boca Raton, Fla, 2004. [ Links ]

16. W. G. Davids, Z. M. Wang, G. Turkiyyah, J. Mahoney and D. Bush, “Finite Element Analysis of Jointed Plain Concrete Pavement with EVERFE 2.2.” Transportation Research Record 1853: Journal of the Transportation Research Board, TRB, National Research Council, United States of American, pp. 92-99, 2003. [ Links ]

17. W. G. Davids, Foundation Modeling for Jointed Concrete Pavements. Transportation Research Record: Journal of the Transportation Research Board , vol. 1730, pp. 34-42, 2000. [ Links ]

18. W. G. Davids, EverFE Theory Manual. United States of American, University of Maine Dept. of Civil and Environmental Engineering, 2003. [ Links ]

19. W. G. Davids, G. M. Turkiyyah and J. P. Mahoney, “Ever-FE: Rigid Pavement Three-Dimensional Finite Element Analysis Tool”. In Transportation Research Record 1629, TRB, National Research Council, United States of American, 1998, pp. 41-49. [ Links ]

20. M. R. Pallares and J.A. Pulecio, “Effect of a dual tire pressure on the design parameters of thick asphalt pavements using finite element freeware,” Dyna, No. 196, pp. 194-203, Apr. 2016. [ Links ]

21. R.G. Packard and S.D. Tayabji. “New PCA thickness design procedure for concrete highway and street pavements”, The Third International Conference on Concrete Pavement Design and Rehabilitation, Purdue University, West Lafayette, United States of American, 1985. [ Links ]

Recibido: 23 de Enero de 2017; Aprobado: 15 de Mayo de 2017

Creative Commons License Este es un artículo publicado en acceso abierto bajo una licencia Creative Commons