Pavimentos Elementos Finitos

Tecnura
ISSN: 0123-921X
tecnura@udistrital.edu.co
Universidad Distrital Francisco Jos de Caldas
Colombia
RODRGUEZ CALDERN, WILSON; PALLARES MUOZ, MYRIAM ROCO

Desarrollo de un modelo de elementos finitos para el diseo racional de pavimentos
Tecnura, vol. 9, nm. 17, 2005, pp. 25-37
Universidad Distrital Francisco Jos de Caldas
Bogot, Colombia
Disponible en: http://www.redalyc.org/articulo.oa?id=257021014002
Cmo citar el artculo

Nmero completo
Sistema de Informacin Cientfica
Ms informacin del artculo Red de Revistas Cientficas de Amrica Latina, el Caribe, Espaa y Portugal
Pgina de la revista en redalyc.org Proyecto acadmico sin fines de lucro, desarrollado bajo la iniciativa de acceso abierto
con-ciencias
Desarrollo de un modelo de elementos
finitos para el diseo racional de pavimentos
Development of a finite elements model for rational
design of pavement
WILSON RODRGUEZ CALDERN

Ingeniero Civil de la Universidad Industrial de Santander, Especialista en Gerencia Integral de
Obras de la Escuela de Ingenieros Militares, Magster en Mtodos Numricos para el Clculo y
Diseo en Ingeniera de la Universidad Politcnica de Catalua (Barcelona). Director de los grupos
de investigacin Informtica y Mtodos Matemticos (registrado ante Colciencias, Colombia) y
Simulacin y Control Numrico SICON (reconocido por Colciencias). Docente del programa de
Ingeniera Civil de la Pontificia Universidad Javeriana.
wilroca50@hotmail.com
MYRIAM ROCO PALLARES MUOZ

Ingeniera Civil de la Universidad Industrial de Santander, Magster en Mtodos Numricos para
Ingeniera de la Universidad Politcnica de Catalua.
mrocio03@hotmail.com
Clasificacin del artculo: investigacin

Fecha de recepcin: mayo 02 de 2005 Fecha de aceptacin: diciembre 16 de 2005
Palabras clave: elementos finitos, diseo racional de pavimentos, modelos constitutivos.

Key words: finite elements, rational design of pavements, constitutive models.
RESUMEN
lineal de los modelos constitutivos de los materia-
El modelo desarrollado pretende predecir los es- les que componen la estructura.
fuerzos y deformaciones de la estructura del
pavimento por accin del trfico y el clima, reali- Una vez obtenido el modelo computacional, es po-
zando inicialmente un anlisis estructural lineal y sible proceder en una segunda fase a la calibracin
planteando luego una propuesta no lineal con su de los parmetros mediante medicin directa u op-
respectiva solucin iterativa, teniendo en cuenta el timizacin matemtica. Este trabajo es el comienzo
carcter de las ecuaciones empleadas. El artculo de una investigacin con posibles modificaciones
presenta un modelo de elementos finitos para des- futuras, dadas las mltiples irregularidades de una
cribir de forma real y lgica el comportamiento estructura de pavimento tradicional; no obstante,
fsico-mecnico de una estructura de pavimento, el propsito es avanzar en la metodologa de anli-
teniendo el mayor nmero de variables o parme- sis del estado tenso-deformacional de los
tros permitidos por las formulaciones lineal y no pavimentos ante cargas estndar.

W ILSON R ODRGUEZ C ALDERN / M YRIAM R OCO P ALLARES M UOZ 25
con-ciencias
ABSTRACT structure, keeping the biggest number of variables
or allowed parameters, so much for the lineal for-
The model developed in this investigation tries to mulation as for the non lineal formulation of the
predict the stresses and strains of the pavement materials constitutive models that compose the
structure to traffic and climate. Initially, a lineal struc- structure. This work starts a research with possi-
tural analysis is made; later, a non lineal proposal ble future modifications, considering the high variety
and the iterative solution are suggested, keeping the of irregularities of the traditional pavement structu-
character of the used equations. This article intro- re. We wish to advance in the tension-deformation
duces a finite elements model that tries to describe state of pavements in front of standard charges
the physical-mechanical behavior of a pavement analysis methodology.
***
1. Antecedentes e introduccin caracterizacin de parmetros de entrada; c) hip-

En el diseo de pavimentos es viable suponer solu- tesis sobre dominio y condiciones de contorno del
ciones a partir de las que, mediante un proceso de problema; d) compatibilidad existente entre capas;
ajustes sucesivos, puede lograrse una condicin e) solicitaciones y condiciones ambientales para el
favorable en la estructura del pavimento que ga- funcionamiento. En la medida en que el modelo se
rantice unas condiciones de equilibrio durante adecue a las condiciones reales, ser mayor la
determinado perodo de tiempo. Este procedimien- aproximacin con la cual se caracterice la respuesta
to para el dimensionamiento de las estructuras de estructural y los esfuerzos y deformaciones crti-
pavimento se conoce como mtodos analticos o cas para el diseo de espesores.
racionales; ellos se fundamentan en el anlisis del
2. Marco terico
estado de esfuerzos y deformaciones que les indu-
cen diferentes solicitaciones como trnsito, clima y Dentro de los mtodos ms difundidos para el dise-
caractersticas de los materiales de cada una de las o racional de pavimentos se encuentran aquellos
capas y del suelo de fundacin o subrasante. que representan la estructura como un conjunto de
capas de espesor constante y longitud infinita hori-
Para el anlisis de estructuras de pavimento exis- zontal, apoyadas sobre un dominio semi-infinito o
ten dos grandes vertientes o grupos de mtodos: subrasante. Las capas se suponen compuestas por
los de carcter emprico y los fundamentados en materiales homogneos, istropos y gobernados por
criterios racionales, comnmente llamados meca- un modelo constitutivo lineal elstico. Sin embargo,
nicistas. Para la implementacin adecuada de un debe notarse que en la definicin del mdulo elstico
modelo racional se requiere de un computador que de la capa asfltica se tienen en cuenta la tempera-
desarrolle los procesos matemticos y del criterio tura y frecuencia de aplicacin de la carga, con lo
ingenieril adquirido por experiencia directa en ca- cual, implcitamente, se considera el comportamien-
racterizacin de los materiales, determinacin de to visco-elstico real de las mezclas asflticas. Bajo
patologas causantes de dao en estructuras del esta perspectiva se han desarrollado importantes
pavimento y familiaridad con rdenes de magnitud aplicaciones computacionales, como ALIZE III,
de parmetros de entrada y salida del problema. ELSYM5 y DEPAV1 ; paquetes de software que
trabajan con las restricciones propias de su formu-
Los mtodos racionales pueden presentar diferen-
cias en la prediccin del comportamiento de la
estructura del pavimento de acuerdo con: a) rigu-
1
rosidad de modelos constitutivos de material; b) Software comercial para el diseo de pavimentos.
26 Tecnura ao 9 No.17 segundo semestre de 2005

con-ciencias
lacin y estn restringidos en capacidad en cuanto b) Introduccin de la no linealidad en las deforma-
al nmero de capas que pueden representar. ciones (large strain) a travs de la activacin
de la no linealidad geomtrica en el modelo
Las tcnicas ms avanzadas para el anlisis y dise- elastoplstico.
o racional de pavimentos consideran el uso de
c) Inclusin de la cohesin y del ngulo de fric-
elementos finitos como la va de mayor flexibilidad
cin como parmetros caractersticos de los
y versatilidad para la representacin de los mate-
geomateriales; esto permite, a travs de un
riales, las condiciones de compatibilidad y la
modelo constitutivo de Drucker-Prager, carac-
representacin de n capas; en general se destaca
terizar el material no linealmente (aproximacin
la gran capacidad para incluir complejidades al pro-
elastoplstica al modelo de Mohr-Coulomb).
blema, generando mayor certeza de los resultados
alcanzados. A este grupo pertenecen las aplicacio- d) Curvas de esfuerzo-deformacin (modelos reo-
nes multifsicas que originalmente no estn diseadas lgicos).
para abordar el problema particular del comporta- e) Condicin de compatibilidad ligada o no ligada
miento tenso-deformacional de la estructura del entre capas.
pavimento; sin embargo, cuentan con todas las he- f) Solucin iterativa e incremental para la aplica-
rramientas que aprovechadas por un usuario cin de la carga.
experto pueden conducir a modelos satisfactorios,
g) Determinacin de zonas plsticas.
ms centrados en la realidad del comportamiento
del pavimento. En este artculo se plantea el uso de h) Inexistencia de restriccin en el nmero de ca-
ANSYS como modelador de alto rango; no obstan- pas constituyentes del modelo geomtrico.
te, pueden emplearse otros modeladores similares, i) Posibilidad futura de tener en cuenta modelos
por ejemplo, ABAQUS o NASTRAN2 , teniendo constitutivos de material y anlisis ms avanza-
en cuenta que en cada uno de ellos es necesario dos, dado el carcter abierto del mtodo.
realizar una adaptacin de acuerdo con las posibili-
dades de modelacin numrica que ofrecen. 2.2 Clasificacin y fundamentacin del anlisis
En este trabajo se presenta un contraste entre los Los modelos fundamentan su anlisis en curvas de
modelos elsticos y elastoplsticos, generando gran- esfuerzo-deformacin, teora elstica y elasto-plas-
des expectativas sobre la validez de las dos teoras ticidad perfecta en geomateriales; sus clculos se
en la representacin del comportamiento de los di- basan en una representacin de las condiciones de
ferentes materiales de un pavimento. campo en trminos del estado de tensiones y de-
formaciones inducidas por diferentes solicitaciones
(cargas, clima, etc.) a los materiales constitutivos
2.1 Ventajas del mtodo de elementos finitos en
el anlisis de la estructura del pavimento.
El mtodo de elementos finitos cuenta con innume- En el modelo lineal es importante sealar que la
rables ventajas para el desarrollo de un modelo imposicin o no de condicin de liga entre capas
conceptual y numrico de anlisis racional de pavi- marca notablemente el anlisis del comporta-
mentos. Entre ellas pueden destacarse: miento del pavimento. Igualmente, en el anlisis
estructural no lineal deben destacarse los siguien-
a) Inclusin de la temperatura como parmetro en
tes efectos:
la capa asfltica. Ello se logra gracias a la posi-
bilidad de incorporar curvas maestras que Grandes deformaciones (large strain): la de-
caracterizan el mdulo dinmico de la mezcla formacin no est definida en forma lineal sino
asfltica en el intervalo de 5 a 40C, u otro in-
tervalo de trabajo. 2
Software comercial de elementos finitos.

con-ciencias
en trminos de una funcin logartmica (e = Ln Figura 1. Estructura tpica de un pavimento
(DL/L)) que permite minimizar el error en los
clculos ante grandes deformaciones. Es de
notar que esta expresin es muy conveniente aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa
aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa
en el anlisis de geomateriales, segn la litera- aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa
aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa
tura tcnica sobre este tema. aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa
aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa
Variacin en la forma de aplicacin de la carga

a medida que la estructura se deforma: las car- aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa

gas concentradas y el peso mantienen su

direccin despus de la deformacin, mientras
que la presin de inflado siempre se ejerce de
manera perpendicular a la superficie, ajustn- Las dimensiones de las capas constitutivas de la
dose a la estructura deformada. estructura del pavimento modelado se presentan en
la tabla 1.
Variacin de la respuesta estructural de acuer-
do con el nivel de esfuerzos: cada punto del Tabla 1. Espesores de capa de la estructura del pavimento
modelo puede trabajar de manera elstica o
inelstica de acuerdo con el nivel de esfuerzo
solicitado, representando la no linealidad del
problema.
3. Planteamiento del problema
3.1 Modelacin de la estructura del pavimento La subrasante es semi-infinita; sin embargo, en los
En el anlisis de la estructura de un pavimento nue- modelos de elementos finitos las fronteras se ubi-
vo se requiere de un predimensionamiento, esto es, can de tal forma que puedan establecerse lneas en
de los espesores de las capas que lo componen donde no existen desplazamientos importantes, de-
obtenidos de diseos preliminares. En estructuras bido a la lejana que presentan respecto a la carga
existentes los espesores de capa son determinados aplicada.
mediante inspeccin directa a travs de sondeos,
3.2 Mdulos de elasticidad de los materiales
apiques o extraccin de ncleos.
Los mdulos de elasticidad del suelo y de las capas
En la figura 1 se esquematiza una estructura de granulares se determinan por medio de ensayos
pavimento tpica. En un modelo de este tipo desa- triaxiales dinmicos para suelos. El mdulo elstico
rrollado con ANSYS pueden aadirse capas y variar de la capa asfltica se establece a travs de ensa-
la dimensin de los espesores y tipos de material, si yos de mdulos dinmicos a la temperatura y a la
el analista as lo requiere. El modelo desarrollado frecuencia de diseo del pavimento, o por medio de
en este trabajo representa slo la mitad derecha pistas de prueba o tramos experimentales sobre los
utilizando una formulacin axisimtrica bidimensio- cuales se realizan ensayos de viga Benkelman do-
nal del problema de esfuerzo-deformacin. La lnea ble o utilizando deflectmetros.
roja representa el eje de simetra para el dominio y
la carga circular de magnitud q y radio a. En cuan- El ensayo de mdulos dinmicos se realiza emplean-
to a la condicin de unin de las interfaces puede do una mquina de traccin indirecta, como la
optarse por compatibilidad ligada o no ligada, de mquina dinmica de ensayos Nottingham Asphalt
acuerdo con los materiales componentes y la expe- Tester NAT. El ensayo es una prueba de traccin
riencia del analista en proyectos similares. indirecta: una carga vertical es aplicada con de-

con-ciencias
terminada frecuencia a una probeta de mezcla as- Los mdulos de elasticidad tomados como referen-
fltica midiendo instantneamente el desplazamiento cia se muestran en la tabla 2.
horizontal de la muestra mediante sensores (LVDT)
Tabla 2. Mdulos de las capas constitutivas
calibrados para capturar desplazamientos del orden
micromtrico. La frecuencia de carga es controlada
a travs del parmetro rise-time (tiempo pico), que
determina el tiempo de ascenso de la carga y permi-
te tener en cuenta el comportamiento viscoelstico
real de las mezclas asflticas. Generalmente, el en-
sayo se realiza a tres temperaturas (5, 25 y 40C)
para caracterizar lo que comnmente se denomina 3.3 Relaciones de Poisson
curva maestra de la mezcla asfltica o del ncleo El mdulo de Poisson se define como la relacin
extrado in situ. Las figuras 2, 3 y la grfica 1 ilus- entre las deformaciones horizontal y vertical indu-
tran ciertos detalles sobre el ensayo de mdulos. cidas en la estructura del pavimento por accin de
la carga. La variacin de este parmetro no tiene
Figura 2. Montaje de la prueba de mdulos dinmicos NAT
una incidencia considerable en su comportamiento;
por esta razn, se emplean valores caractersticos
como los mostrados en la tabla 3.
Tabla 3. Relaciones de Poisson de las capas constitutivas de

la estructura del pavimento
Figura 3. Aplicaciones de carga en pruebas de mdulos dinmicos

* Los datos fueron tomados del documento Caracterizacin
de materiales viales. Su aplicacin al diseo estructural de
pavimentos de la Universidad del Cauca, en el cual se recopilan
importantes trabajos de investigacin de diferentes pases
incluido Colombia.
En el modelo que aqu se propone se incorpor la

variacin del mdulo de Poisson respecto a la tem-
Grfica 1. Curvas maestras de mdulos dinmicos

3
El valor del mdulo dinmico fue tomado de ensayos sobre
muestras obtenidas en el Laboratorio de Asfaltos del Instituto
Colombiano del Petrleo, en el marco de la investigacin de los
autores, titulada Aditivos mejoradores de adherencia para asfaltos.
Las conclusiones ms importantes de este proyecto se publicaron
en la revista de Corasfaltos Asfaltos y Pavimentos (septiembre de
2000). La granulometra de la mezcla corresponde al huso medio de
la especificacin MDC-2 del Invas-96 para capa de rodadura de
mezcla densa en caliente. El asfalto utilizado fue una muestra de
Apiay-1999. La dosificacin se efectu con el mtodo Marshall y el
porcentaje ptimo de asfalto obtenido fue de 6,3%. El valor del
mdulo dinmico empleado en el modelo se tom a una temperatura
de 25C y una frecuencia de 5 Hz.

con-ciencias
peratura en la capa asfltica utilizando la siguiente Grfica 3. Curva esfuerzo-deformacin de un geomaterial
funcin, propuesta por la Universidad de Notting- (modelo Drucker-Prager)
ham:
Para T< 10 C = 0,25
T> 10 C = 0,01*T + 0,15 (1)
Suponiendo condiciones normales de temperatura

(25C) en la capa asfltica, 0,40. La grfica 2
ilustra el comportamiento de Poisson en la capa
asfltica respecto a la temperatura, como lo indica
la expresin (1).
Grfica 2. Curva de la Universidad de Nottingham para la Para el modelo no lineal se establecieron valores
relacin de Poisson en mezclas asflticas
de prueba obtenidos de la literatura tcnica en el
tema. Sin embargo, es importante sealar que este
modelo es una propuesta cuyo xito depender de
la calibracin de los parmetros por medio de en-
sayos de laboratorio y de campo que permitan
establecer su magnitud real. En cualquier caso, el
modelo Drucker-Prager ha sido probado con xito
en problemas geotcnicos de cierta utilidad prcti-
ca en ingeniera. La tabla 4 muestra los valores
adoptados para probar el modelo no lineal de ele-
mentos finitos en este trabajo.
3.4 ngulo de friccin interna y cohesin Tabla 4. Cohesin y ngulo de friccin para el modelo no
lineal propuesto
Estos dos datos son imprescindibles para el modelo
elastoplstico no lineal propuesto en este trabajo.
Empleando estos dos parmetros es posible definir
el esfuerzo de fluencia incorporando materiales
constitutivos del tipo elasto-plstico perfecto (gr-
fica 3), que son una aproximacin bastante cercana Se conoce que en la capa asfltica el ligante aporta
del modelo constitutivo de Mohr-Coulomb para geo- cohesin a la mezcla asfltica y el agregado mine-
materiales. ral la componente friccionante (SUPERPAVE). En
las dems capas la cohesin y el ngulo de friccin
Para geomateriales, el modelo de Drucker-Prager
son parmetros conocidos y utilizados ampliamen-
establece que el esfuerzo de fluencia puede ser
te en problemas de esfuerzos bajo zapatas y
calculado por medio de:
estabilidad de taludes, referenciados con suficien-
(2) cia en la literatura tcnica sobre el tema.
donde: y Esfuerzo de fluencia 3.5 Presin de inflado y radio de carga

Angulo de friccin interna La carga aplicada est representada por una pre-
C Cohesin > 0 sin vertical distribuida uniformemente sobre un rea

con-ciencias
circular. Se ha demostrado que la presin transmi- tiendan a cero. Las capas se consideran ligadas o no
tida equivale a la presin de inflado de las llantas. ligadas, segn el caso, y la geometra utilizada busca
Como el modelo utilizado en este trabajo es bidi- aprovechar la caracterstica cilndrica del modelo real,
mensional, en el programa es necesario especificar permitiendo un ahorro considerable en elementos y
que la carga es aplicada en un rea circular, lo cual tiempo computacional.
determina la naturaleza axisimtrica del modelo. El
radio de carga empleado es de 10,8 cm, bajo una 3.7 Parmetros de control del diseo
presin de inflado de 5,6 kg/cm2.
En la estructura del pavimento, la deformacin por
traccin en la fibra inferior de la capa asfltica
3.6 Condiciones de contorno aplicadas al
(et) y el esfuerzo vertical por compresin aplica-
modelo
do sobre la subrasante (sv) se establecen como
El modelo est confinado lateralmente y el eje de la parmetros de control. Se destaca que el esfuerzo
carga circular est restringido en sus movimientos de traccin en la fibra inferior de la capa asfltica
horizontales debido a la simetra de la carga respecto (st) es un parmetro que puede establecer diferen-
al eje central. La base de la estructura est restringida cias entre un modelo elstico y uno elastoplstico.
en su movimiento vertical y su ubicacin se establece En la figura 4 se indica la ubicacin de estos par-
en la medida que los desplazamientos de dicha base metros en la estructura del pavimento.
Figura 4. Ubicacin de los parmetros de control del diseo de la estructura del pavimento
3.8 Organizacin del programa ANSYS 3.9 Elemento finito empleado
El programa ANSYS est organizado en dos niveles El elemento PLANE424 es bidimensional y definido
bsicos: entrada y procesamiento (o de rutinas). El por cuatro nodos, con dos grados de libertad por nodo
nivel de entrada acta como una puerta dentro y fue- (UX,UY) y capacidad para representar estados de
ra del programa; es usado para tener un control global tensin plana, deformacin plana y axisimetra. El
del programa. En el nivel de procesamiento estn dis- elemento posee propiedades plsticas, creep, rigidez
ponibles varios procesadores, cada uno constituido por al esfuerzo, grandes deflexiones y capacidad para
un conjunto de funciones que realizan tareas especfi- grandes deformaciones; tambin posee propiedades
cas del anlisis: en el preprocesador general (/PREP7) para materiales ortotrpicos5. En la tabla 5 se ilus-
se construye el modelo paramtricamente (geome- tran la geometra, localizacin de nodos y sistema de
tra, materiales, acoplamiento y restricciones a los coordenadas para este elemento, con el que se cons-
grados de libertad), en el procesador de solucin truyeron todas las capas del modelo.
(/SOLU) se aplican las cargas y se obtiene la solu-
4
cin del problema, y en el posprocesador general Disponible en la librera de elementos del programa comercial de
(POST1) se evalan los resultados sobre todo el elementos finitos ANSYS.
5
modelo en puntos especficos. Las direcciones ortotrpicas del material corresponden a las
direcciones de las coordenadas del elemento.

con-ciencias
Tabla 5. Datos de entrada del elemento para anlisis estructural 2D, PLANE42
4. Resultados y discusin mecnica de los geomateriales y pueden ser

Con el objeto de observar los efectos de los mode- determinados en el laboratorio.
los constitutivos de los materiales, las condiciones Los resultados de esfuerzo y deformacin de los
de liga o no liga entre capas y la verificacin de tres modelos en consideracin se presentan por
parmetros de diseo (esfuerzos y deformaciones medio de isocontornos.
de traccin en la interfase capa asfltica-base
granular y los esfuerzos verticales en la subrasan- 4.1 Modelo elstico lineal de capas ligadas
te) en la estructura de un pavimento, se realizaron
tres modelos6 : En la figura 5 se ilustra el modelo geomtrico para
las tres alternativas planteadas anteriormente. De
1) Modelo elstico lineal de capas ligadas: la hip- abajo hacia arriba, la primera capa es la subrasan-
tesis de compatibilidad en todas las capas es de te; la siguiente, la sub-base; la tercera, la base y la
tipo ligada, lo que significa compatibilidad total cuarta, la capa asfltica.
en los desplazamientos horizontales y verticales.
2) Modelo elstico lineal de capas no ligadas: la
hiptesis de compatibilidad en todas las capas
es de tipo no ligada, que simboliza compatibili- 6
Trabajo de investigacin realizado por uno de los autores
dad slo en los desplazamientos verticales. (Rodrguez W.) en el Instituto Colombiano del Petrleo.
7
3) Modelo elastoplstico no lineal de capas liga- El anlisis de comparacin se realiza entre un modelo lineal y un
elastoplstico (no se hace con un viscoelstico), dado que existen
das7 : la hiptesis de compatibilidad en todas las
serias restricciones para obtener parmetros viscoelsticos de
capas es de tipo ligada, con compatibilidad total mezclas asflticas que puedan emplearse en un modelo de
en los desplazamientos horizontales y vertica- elementos finitos. Hasta el momento, en Colombia no existen
les. Este modelo es una propuesta que debe ser mquinas para realizar ensayos que puedan obtener estos
parmetros va experimental. El modelo elastoplstico sugerido
calibrada y verificada en estudios futuros; sin
tiene la ventaja de poseer parmetros obtenibles en el laboratorio y,
embargo, es factible numrica y fsicamente, si por tanto, es un modelo factible que va ms all de las tcnicas
se tiene en cuenta que los parmetros de cohe- convencionales y que verifica un problema reportado en la literatura:
sin y ngulo de friccin son comunes en la la plastificacin de la base.

con-ciencias
Figura 5. Dominio del modelo geomtrico Grfica 5. Isocontornos de deformacin total horizontal
En las grficas 4 y 5 puede observarse la concentra-

cin de esfuerzos y deformaciones horizontales en Grfica 6. Isocontornos de esfuerzo vertical
la interfase constituida por la capa asfltica y la base,
como lo predicen otras teoras (por ejemplo, la els-
tica lineal multicapa, empleada por programas
comerciales como el ALIZE III y el DEPAV). La
grfica 6 revela que el bulbo de esfuerzos reduce su
influencia a medida que se acerca a la subrasante,
como se espera para garantizar la capacidad portan-
te del suelo. Finalmente, la grfica 7 permite observar
una concentracin de esfuerzos de corte en la capa
asfltica, que comnmente no es empleada como pa-
rmetro de control, pero que sera importante evaluar
por los posibles efectos de dao que le produciran,
ya que de ser as debera tenerse en cuenta en el
diseo racional de pavimentos.
Grfica 7. Isocontornos de esfuerzo cortante
Grfica 4. Isocontornos de esfuerzo horizontal

con-ciencias
4.2. Modelo elstico lineal de capas no ligadas Grfica 10. Isocontornos de esfuerzo vertical
Las grficas 8 y 9 ponen de manifiesto la influencia

de la compatibilidad no ligada en el comportamien-
to estructural del pavimento, demostrando que los
esfuerzos y las deformaciones horizontales son ab-
sorbidos en gran parte por la capa asfltica. La grfica
10 demuestra que el bulbo de esfuerzos vertical se
corta y solo permanece en la capa asfltica, deno-
tando de cierta manera la discontinuidad provocada
por el carcter no ligado de las capas. La grfica 11
muestra la persistencia de la concentracin de
esfuerzos cortantes en la capa asfltica.

Grfica 9. Isocontornos de deformacin total horizontal 4.3. Modelo elastoplstico no lineal de capas
ligadas

con-ciencias
Grfica 13. Isocontornos de deformacin total horizontal Grfica 16. Isocontornos de deformacin plstica horizontal
Grfica 14. Isocontornos de esfuerzo vertical En las grficas 12 y 13 se observa la concentracin

de esfuerzos y deformaciones horizontales en la
interfase capa asfltica-base, tal como lo predice
el modelo elstico lineal de capas ligadas; sin em-
bargo, los valores difieren un poco por el cambio
mismo del modelo constitutivo. La grfica 14 mues-
tra un bulbo de esfuerzos verticales similar al del
modelo elstico de capas ligadas. La grfica 15
expone la misma configuracin en la concentracin
de los esfuerzos de corte de la capa asfltica del
modelo lineal elstico de capas ligadas, con dife-
rencias en los valores. Finalmente, la grfica 16
revela los fenmenos plsticos en la capa de base,
los cuales se sugiere analizar con mucho cuidado,
si se tiene en cuenta que ellos estn controlados
por los parmetros de cohesin y ngulo de fric-
cin del material granular constitutivo.
5. Comparacin de resultados
En la tabla 6 se presentan comparaciones entre
esfuerzos y deformaciones para cada uno de los
modelos. Sx (kg/cm2) representa el esfuerzo hori-
zontal en la interfase capa asfltica y base granular;
ex, la deformacin horizontal en la interfase capa
asfltica y base granular; Sxy (kg/cm2) es el es-
fuerzo cortante mximo en la capa asfltica.

con-ciencias
Tabla 6. Comparacin de parmetros de respuesta en los tres modelos
El esfuerzo horizontal Sx es mayor en el modelo no costosas tcnicas de retroclculo, debido a la gran

lineal y se observa una diferencia significativa en- informacin que requieren para obtener ajustes me-
tre los resultados obtenidos con condicin ligada y dianamente satisfactorios, y al consumo de tiempos
no ligada. En las deformaciones horizontales (x) considerables para llegar a obtener resultados con
ocurre lo mismo que con los esfuerzos horizonta- mrgenes de confiabilidad aceptables. En este sen-
les, Sx. En los esfuerzos de corte Sxy se observan tido, se considera que la propuesta no lineal
valores similares entre el modelo elstico lineal li- presentada en este trabajo es viable, si se tiene en
gado y el modelo no lineal, mientras que el modelo cuenta que los parmetros de cohesin, friccin y
lineal no ligado presenta esfuerzos menores. La con- mdulo pueden ser determinados con pruebas de
dicin no ligada es difcil de garantizar en campo, laboratorio reconocidas.
dado que en las interfases en verdad existe un fe-
nmeno de contacto complejo y, en algunos casos, 6. Conclusiones y recomendaciones
la condicin de no liga puede ser inconveniente para El mtodo de elementos finitos demuestra su
la estabilidad de la estructura. No obstante, los mo- superioridad respecto a otras metodologas y
delos evidencian que el efecto estructural de las marca un camino a seguir para realizar diseos
cargas es menor cuando la condicin de compatibi- de pavimentos ms ajustados al comportamiento
lidad de las capas es no ligada. real. En una segunda fase se propone realizar
En general, el modelo no lineal es ms crtico que el las actividades de calibracin pertinentes y la
lineal. Esta situacin genera la pregunta: es ne- verificacin de los modelos para determinar su
cesario avanzar hacia modelos no lineales de grado de validez. Sin embargo, debe destacar-
prediccin de esfuerzos y deformaciones en pavi- se que el desarrollo de los modelos numricos
en s mismos constituyen una importante herra-
mentos? Esta pregunta no puede ser respondida
mienta para futuras investigaciones en el campo
categricamente, descalificando las aproximacio-
del anlisis y diseo de pavimentos.
nes hechas por los modelos lineales; no obstante,
los resultados obtenidos del anlisis que aqu se plan- El esfuerzo de corte es un dato no empleado
tea advierten que vale la pena realizar esfuerzos normalmente como parmetro de control en el
por ir ms all de la teora de la elasticidad, ya que, diseo de pavimentos; sin embargo, los mode-
efectivamente, en los pavimentos suceden fenme- los de elementos finitos advierten la existencia
nos no lineales como deformaciones plsticas, de concentraciones de esfuerzos de corte que
fluencia y fatiga, entre otros, que slo pueden abor- probablemente induzcan fallas importantes a la
darse considerando otros modelos constitutivos de estructura. Mediante la implementacin de un
mayor alcance. modelo numrico adecuado podra predecirse
de manera ms clara su influencia y la relevan-
Sin embargo, deber evitarse emplear modelos ex- cia de tenerlo en cuenta en el diseo. Para pistas
cesivamente complejos, pues en muchas ocasiones de aeropuertos, el clculo de esfuerzos de cor-
no es factible obtener sus parmetros por medio te tiene gran importancia si se tiene en cuenta
del laboratorio y no queda otra salida que acudir a que el impacto de los aviones durante el proce-

con-ciencias
so de aterrizaje produce esfuerzos de corte de y econmica de emplear este tipo de materiales
magnitud considerable. en estructuras.
La condicin de liga entre capas no es muy Se recomienda usar los resultados de los mode-
ventajosa respecto de la de no liga desde el punto los numricos de elementos finitos con buen
de vista de los esfuerzos y deformaciones ob- juicio y rigor matemtico, fsico e ingenieril para
tenidos del anlisis. No obstante, en la realidad tomar decisiones basadas en las restricciones e
suele existir un fenmeno complejo de contac- hiptesis propias de los modelos; de otra mane-
to entre capas intermedio entre las condiciones ra, puede llegarse a concepciones erradas o poco
de liga y no liga. En la topografa colombiana ajustadas a la realidad. No debe olvidarse que
(zonas montaosas en gran parte del territorio) todo modelo matemtico o fsico es una aproxi-
no convienen los estados de no liga, que pue- macin a la realidad, y que es imprescindible
den afectar la estabilidad de las estructuras. reconocer el alcance de las teoras que susten-
Algunas veces la condicin de liga entre capas tan las predicciones halladas.
puede verse afectada por la inclusin de geoma- El modelo elastoplstico desarrollado en este
llas en las interfases (por ejemplo, geomallas trabajo predice deformaciones plsticas en la
en la interfase capa asfltica-base); por esto, base granular. Esta situacin es reportada por
se recomienda evaluar la conveniencia tcnica la literatura en casos de campo.
Referencias bibliogrficas
[1] CHEN, W. F. Y BALADI, G. Y. (1985). Soil Plasticity. [9] OWEN, D. R. J. Y HINTON, E. (1980). Finite ele-
Theory and implementation. Elsevier. ments in plasticity. Pineridge Press.
[2] CHRISTIAN, J. F. Y DESAI, C. S. (1977). Constituti- [10] PERZYNA, P. (1966). Fundamental problems in vis-
ve laws for geologic media. Numerical Methods in coplasticity. Advanced Applied Mechanics, vol. 9, pp.
Geotechnical Engineering. Desai, C. S. & Christian, J. 243-377.
T. (eds.), McGraw-Hill. [11] SCOTT, R. F. (1985). Plasticity and constitutive rela-
[3] DAFALIAS, H. F. Y HERRMANN. L. R. (1982). tions in soil mechanics. J. Geotech. Engrg. ASCE, vol.
Bounding surface formulation of soil plasticity. Soil 111, pp. 559-605.
Mechanics Transient and cyclic loads. Pande, G.N. & [12] CHEN, W. F. Y MIZUNO, E. (1990). Nonlinear analy-
Zienkiewicz, O.C. (eds.) Wiley. sis in Soil Mechanical. Theory and implementation.
[4] DAVID M. POTTS, LIDIJA ZDRAVKOVIC (1999). [13] SWANSON ANALYSIS SYSTEMS, Inc. Houston.
Finit Element Analysis in Geotechnical Engineering - Ansys User Manual/Revision 5.0. vol. I (Procedu-
Theory. Thomas Telford, London. res), vol. II (Comands), vol. III (Elements), vol. IV
[5] DESAI, C. S. Y SIRIWARDANE, H. J. (1984). Cons- (Theory).
titutive laws for engineering materials. Prentice Hall, [14] WOOD, D. M. (1984). Choice of models for geotech-
Englewood Cliffs, New Jersey. nical problems. Mechanics of engineering materials.
[6] DRUCKER, O. C. Y PRAGER, N. (1962). Soil me- Desai, C. S. & Gallagher, R. H. (eds.), Wiley.
chanics and plastic analysis of limit design. Quarterly [15] WOOD, D. M. (1990). Soil Behaviour and Critical
Applied Math, vol. 10, pp. 157-164. State Soil Mechanics. Cambridge University Press.
[7] GUDEHUS, G. (1985). Requirements for constitutive [16] ZIENKIEWICZ, O. C. y CORMEAU, I. (1974).
relations for soils. Mechanics of Geomaterials. Rocks, Viscoplasticity-plasticity and creep in elastic solids -
Concrete, Soils. Bazant, Z. P. (ed.) Wiley, cap. 4. a unified numerical solution approach. Internatio-
[8] MARQUES, J. M. M. C. (1984). Stress computation nal Journal for Numerical Methods in Engineering.
in elastoplasticity. Engineering Computations, vol. 1, vol. 8, pp.
pp. 42-51.


Pavimentos Elementos Finitos

Cargado por

Copyright:

Formatos disponibles

Pavimentos Elementos Finitos

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Pavimentos Elementos Finitos

Cargado por

Copyright:

Formatos disponibles

Tecnura

RODRGUEZ CALDERN, WILSON; PALLARES MUOZ, MYRIAM ROCO

Disponible en: http://www.redalyc.org/articulo.oa?id=257021014002

Cmo citar el artculo

WILSON RODRGUEZ CALDERN

MYRIAM ROCO PALLARES MUOZ

Clasificacin del artculo: investigacin

Palabras clave: elementos finitos, diseo racional de pavimentos, modelos constitutivos.

Desarrollo de un modelo de elementos finitos para el diseo racional de pavimentos

1. Antecedentes e introduccin caracterizacin de parmetros de entrada; c) hip-

26 Tecnura ao 9 No.17 segundo semestre de 2005

Desarrollo de un modelo de elementos finitos para el diseo racional de pavimentos

Variacin en la forma de aplicacin de la carga

a medida que la estructura se deforma: las car- aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa

gas concentradas y el peso mantienen su

3. Planteamiento del problema

28 Tecnura ao 9 No.17 segundo semestre de 2005

Tabla 3. Relaciones de Poisson de las capas constitutivas de

Figura 3. Aplicaciones de carga en pruebas de mdulos dinmicos

En el modelo que aqu se propone se incorpor la

Grfica 1. Curvas maestras de mdulos dinmicos

Desarrollo de un modelo de elementos finitos para el diseo racional de pavimentos

T> 10 C = 0,01*T + 0,15 (1)

Suponiendo condiciones normales de temperatura

donde: y Esfuerzo de fluencia 3.5 Presin de inflado y radio de carga

30 Tecnura ao 9 No.17 segundo semestre de 2005

3.8 Organizacin del programa ANSYS 3.9 Elemento finito empleado

Desarrollo de un modelo de elementos finitos para el diseo racional de pavimentos

4. Resultados y discusin mecnica de los geomateriales y pueden ser

32 Tecnura ao 9 No.17 segundo semestre de 2005

En las grficas 4 y 5 puede observarse la concentra-

Desarrollo de un modelo de elementos finitos para el diseo racional de pavimentos

Las grficas 8 y 9 ponen de manifiesto la influencia

Grfica 8. Isocontornos de esfuerzo horizontal

Grfica 12. Isocontornos de esfuerzo horizontal

34 Tecnura ao 9 No.17 segundo semestre de 2005

Grfica 14. Isocontornos de esfuerzo vertical En las grficas 12 y 13 se observa la concentracin

Desarrollo de un modelo de elementos finitos para el diseo racional de pavimentos

El esfuerzo horizontal Sx es mayor en el modelo no costosas tcnicas de retroclculo, debido a la gran

36 Tecnura ao 9 No.17 segundo semestre de 2005

Desarrollo de un modelo de elementos finitos para el diseo racional de pavimentos

También podría gustarte