Métodos de Análisis Esfuerzos Deformaciones Desplazamientos en Pavimentos Flexibles

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 63

Universidad Nacional de Colombia

Sede Manizales
Facultad de Ingeniera y Arquitectura
Departamento de Ingeniera Civil
Pavimentos. Profesor: Luis Ricardo Vsquez Varela, M.Sc.

Mtodos de Anlisis:

Esfuerzos, deformaciones y desplazamientos


en pavimentos flexibles.
Referencias:

Yoder & Witczak (1975). Principles of Pavement Design.


Huang (2004). Pavement Analysis and Design.
Papagiannakis & Masad (2008). Pavement Design and Materials.
TM 5-822-13/AFJMAN 32-1018 Pavement Design for Roads, Streets and Open Storage Areas, Elastic
Layered Method (Joint Departments of the Army and Air Force, 1994).

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

Desarrollo histrico.
Los modelos matemticos son las herramientas para la
solucin de los problemas de la ingeniera.
Bases de una solucin:
Tener en cuenta los requerimientos fsicos de la estructura para
soportar las cargas externas, las deformaciones y los esfuerzos en los
elementos.

Definir el comportamiento mecnico de los materiales de acuerdo con


las leyes bsicas de la mecnica que gobiernan el movimiento y las
fuerzas.
14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

Semiespacio continuo homogneo,


istropo y elstico lineal de Boussinesq.

La teora original de Boussinesq (1885) se basa en la aplicacin de una carga


puntual concentrada sobre un semiespacio infinito en planta y profundidad.

El semiespacio es continuo, homogneo, istropo y elstico lineal.

Los esfuerzos, deformaciones y desplazamientos bajo un rea uniformemente


cargada se determinan mediante integracin.

Los neumticos de los vehculos se pueden modelar como reas circulares.

El modelo de Boussinesq es aplicable en pavimentos cuya estructura y subrasante


tengan mdulos de elasticidad muy similares. Es un caso muy poco representativo.

Una aplicacin ms realista es la evaluacin del mdulo movilizado de la


subrasante a partir de medidas de deflexin (desplazamiento vertical) bajo una
carga de configuracin conocida (intensidad, forma).

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

Grficos de Foster y Ahlvin (1954) para el


modelo de Boussinesq.
Solucin de un sistema en coordenadas
cilndricas (r, , z) con origen en el centro de
un rea circular superficial de radio (a).

El rea cargada aplica una presin uniforme


(q), es decir, no tiene rigidez propia.
El material se caracteriza por el mdulo
elstico (E) y la relacin de Poisson ().
Se puede obtener:

14/05/2015

Esfuerzo vertical (z).


Esfuerzo radial (r).
Esfuerzo tangencial ().
Esfuerzo cortante (zr).
Desplazamiento o deflexin vertical (uz).

LUIS RICARDO VSQUEZ VARELA, M.Sc.

= 0.5

Los nmeros en las curvas


indican (r / a)

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

Los nmeros en las curvas


indican (r / a)

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

Los nmeros en las curvas


indican (r / a)

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

Los nmeros en las curvas


indican (r / a)

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

Los nmeros en las curvas


indican (r / a).

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

Ejemplo:
Sobre una subrasante nivelada y compactada se aplica una carga igual
a 40 kN sobre una placa flexible de 0.15 metros de radio.
Calcule los esfuerzos producidos por la carga:
Bajo el centro del rea cargada y en la superficie del terreno.
A 0.75 metros del centro del rea cargada y a 0.3 metros de profundidad.

Considere un mdulo de elasticidad de 100,000 kPa y calcule el


desplazamiento vertical en la superficie del terreno bajo el centro del
rea cargada.

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

10

Solucin:
El rea de la superficie cargada es igual a:

= 0.15

La presin uniforme aplicada es igual a:

= 40

= 0.071

0.071

= 566

Se calculan las profundidades (z) y distancias radiales (r) de los puntos


de inters en trminos del radio del rea cargada (a = 0.15 m):

14/05/2015

Punto

z (m)

z/a

r (m)

r/a

0.0

0.0

0.0

0.0

0.30

2.0

0.75

5.0

LUIS RICARDO VSQUEZ VARELA, M.Sc.

11

Se obtienen los factores de influencia de los grficos apropiados para


cada respuesta estructural en cada punto:
Punto

z/a

r/a

z /q (%)

r /q (%)

/q (%)

zr /q (%)

0.0

0.0

100%

100%

100%

0%

1.5

2.0

5.0

0.3%

1.6%

0%

0.68%

0.17

Los esfuerzos y el desplazamiento vertical en cada punto son:

14/05/2015

Punto

z (kPa)

r (kPa)

(kPa)

zr (kPa)

uz (mm)

0.0

0.0

566

566

566

1.27

0.3

0.75

1.7

9.1

3.8

0.14

LUIS RICARDO VSQUEZ VARELA, M.Sc.

12

Ejemplo 2:
Se realiz un ensayo con deflectmetro de
impacto liviano (ASTM E2583 & IAN73 - UK) sobre
un suelo de subrasante.
El equipo aplica la carga a travs de una placa de 0.15
m de radio.
El desplazamiento vertical se mide mediante un
acelermetro o un gefono.
La carga se mide mediante una celda de presin.

Se obtuvo una carga pico de 7.0 kN y un


desplazamiento vertical de 0.5 mm.
Calcule el mdulo elstico movilizado por la
subrasante.
http://www.dynatest.com/equipment/structural/lwd

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

13

Solucin:
Considerando la ecuacin del grfico de deflexin:

. .

. .
=

La deflexin medida por este equipo es superficial y bajo el centro del rea
cargada, de forma que (z/a) = 0 y (r/a) = 0.
Se puede proponer para el mdulo que:

Se obtiene un parmetro F = 1.5.


El mdulo es:
=

14/05/2015

7
0.15

. 0.15 . 1.5

0.5
1,000

= 44,536 /

LUIS RICARDO VSQUEZ VARELA, M.Sc.

14

Teora de capas elsticas.

http://www.columbia.edu/cu/civileng

La primera solucin para un sistema


generalizado multicapa elstico fue
presentada por Donald M. Burmister
en la Universidad de Columbia
(1943).
Sistemas elsticos de N capas.

/ling/burmister/burmister.html

a
q

Capa superior

1, E1, h1
2, E2, h2

Soluciones especficas para sistemas de


dos y tres capas.
Cargas uniformes, aplicadas de forma
normal sobre un rea circular.

Cada capa se compone de materiales


que son istropos, homogneos,
elsticos lineales y carentes de peso
propio.

Capas intermedias

i, Ei, hi
n-1, En-1, hn-1

Capa inferior

n, En, hn =
z

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

15

El sistema de capas elsticas es


compuesto, es decir, puede existir
continuidad de los esfuerzos o
deformaciones a travs de las
interfases entre capas.

La mayor parte de las soluciones para


pavimentos
asumen
que
los
materiales son elsticos lineales.

Esto no es un problema si se estima


adecuadamente el mdulo elstico
de cada material para el nivel de
esfuerzo (subrasante y capas
granulares) o para la temperatura y
frecuencia de carga (capas asflticas).

14/05/2015

Schiffman extendi los trabajos de


Burmister para formas generalizadas
de carga asimtrica, incluyendo
esfuerzos cortantes en la superficie.

Para la aplicacin prctica del sistema


de capas elsticas se han tabulado
coeficientes y se han elaborado
bacos por autores como Jones,
Peattie, Huang y el mismo Burmister.

En la actualidad las soluciones


tabulares o grficas han sido
sustituidas por software.

LUIS RICARDO VSQUEZ VARELA, M.Sc.

16

Solucin de un slido de revolucin con deformacin axialmente simtrica


alrededor de z.
X

Ecuaciones de equilibrio de esfuerzos:


r rz r

0
r
z
r

d
dz

zr z zr

0
r
z
r

Ecuaciones de deformacin:

ur
u
r
r
r
u u
rz r z
z
r

14/05/2015

u z
z

r
rz

LUIS RICARDO VSQUEZ VARELA, M.Sc.

17

Funciones de esfuerzo en
coordenadas cilndricas.

Funciones de desplazamiento
en coordenadas cilndricas.

2
2 2 1

r 2

r
dr r r r

1
ur
E

2 2
2
r
r

2
1 2 1 2 1

2 2

z
r
r r r r r

1
u
E

1 2


2
r
r z


2
2
z 2 2
z
z

1
uz
E

2
2
21 z 2

1 2 2 2


2
r z r r r 2 z 2

1
z
r

rz

Criterio de solucin.

4 0 4 0

2 2
2
1 z 2 r 2


2 1 2
2

r
z 2 r z
2

14/05/2015

2 1
1 2
2

r 2 r r r 2 2 z 2

LUIS RICARDO VSQUEZ VARELA, M.Sc.

18

Soluciones para sistemas de dos capas.

Burmister y Huang proponen bacos para


sistemas de dos capas a partir de la solucin
general del sistema multicapa.

Se resuelve un sistema en coordenadas


cilndricas (r, , z) con origen en el centro del
rea cargada.

Se aplica una presin uniforme (q) sobre un rea


circular de radio (a).

Los bacos permiten obtener respuestas de


inters en el anlisis de pavimentos flexibles:

Distribucin del esfuerzo vertical (z) para h1 = a.

Esfuerzo vertical en la interfase entre la


subrasante y el pavimento bajo el centro de la
carga (c).

Desplazamiento vertical en la superficie bajo el


centro de la carga (uz=0).

Deformacin de traccin en el fondo de la capa 1


bajo el centro de la carga (r).

14/05/2015

Capa 1

1, E1

h1

Capa 2

LUIS RICARDO VSQUEZ VARELA, M.Sc.

1, E2

h2=

19

z / q

Capa 1

Interfase 1-2
z/a

Capa 2

Capa de subrasante

Note que el espesor de


pavimento es igual al
radio del rea cargada.

Capa de
base del
pavimento

z / q

Distribucin del esfuerzo vertical en un sistema bicapa (Burmister, 1958)

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

20

Si el rea cargada es rgida:


=0 =

1.5
2
2

=0 =

1.18
2
2

Deflexin vertical en la superficie de sistemas bicapa (Burmister, 1943)

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

21

Esfuerzo vertical en la interfase de sistemas bicapa (Huang, 1969)

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

22

Factor de deformacin, Fe

Deformacin horizontal de traccin en el fondo de la capa 1 (Huang, 1973)

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

23

Ejemplo:
Un neumtico flexible aplica una presin de 700 kPa en un rea circular de 0.1
metros de radio sobre una capa de concreto asfltico de 0.2 metros de
espesor construida sobre una subrasante de espesor infinito.
Los mdulos de elasticidad de la capa asfltica y la subrasante son 1,400 MPa
y 140 MPa, respectivamente.
Calcule, bajo el centro del rea cargada:
La deflexin superficial (uz=0).
El esfuerzo vertical de la subrasante (c).
La deformacin de traccin en el fondo de la capa asfltica (r).

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

24

Solucin:
Para resolver el problema se requieren los parmetros:
h1 / a
E1 / E2

= 0.20 m / 0.10 m
= 1,400 MPa / 140 MPa

=2
= 10

Deflexin superficial bajo el centro del rea cargada (uz=0):


F2 = 0.3 del grfico de deflexin.

=0 =

14/05/2015

1.5. . . 2
1.5 700 0.1 0.3
=
= 0.24
2
140,000

LUIS RICARDO VSQUEZ VARELA, M.Sc.

25

Solucin:
Esfuerzo vertical en la subrasante (c) bajo el centro del rea cargada:
a / h1
E1 / E2
Se obtiene del grfico

= 0.10 m / 0.20 m
= 1,400 MPa / 140 MPa
: c / q

= 0.5
= 10
= 0.1

= 0.1 . = 0.1 700 = 70

Deformacin horizontal de traccin en el fondo de la capa asfltica.


Factor de deformacin, Fe = 0.35

.
700 0.35
=
=
= 175 106
1
1 400,000

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

26

Soluciones para sistemas de tres capas.

Corresponde al sistema de coordenadas cilndricas(r,


, z) con origen en el centro del rea superficial
cargada con presin uniforme (q) y de radio (a).

Se obtienen los esfuerzos en las interfases 1 y 2.

Interfase 1:

Interfase 1

Anlisis de fatiga en el fondo de las capas asflticas


(r).

Interfase 2:

Anlisis de fatiga en el fondo de las capas [base]


asflticas (r).

Anlisis de deformacin permanente en la parte


superior de la subrasante (z).

En este modelo la traccin es negativa y la


compresin positiva.

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

r3

1
+

1
+

Interfase 2

27

Llanta con carga y


presin especficas.

Capa cementada con asfalto.


(Espesor finito).
Capa de base.
(Espesor finito).
Suelo de subrasante.
(Se asume de espesor infinito).

1.
2.
3.
4.

Ubicacin

Deflexin en la superficie del pavimento.


Deformacin horizontal de tensin en el fondo de la capa asfltica.
Deformacin vertical de compresin en la parte superior de la base.
Deformacin vertical de compresin en la parte superior de la subrasante.

Respuesta

Razn de su empleo

Deflexin (uz)

Empleada para imponer restricciones de


carga durante el descongelamiento de
primavera y en el diseo de sobrecapas.

(2) Fondo de la capa de mezcla


asfltica en caliente

Deformacin horizontal
de traccin (x y)

Empleada para predecir la falla por fatiga


en la mezcla asfltica en caliente

(3) Parte superior de una capa


intermedia (base o subbase)

Deformacin vertical de
compresin (z)

Empleada para predecir la falla por


ahuellamiento en la base o la subbase

(4) Parte superior de la


subrasante

Deformacin vertical de
compresin (z)

Empleada para predecir la falla por


ahuellamiento en la subrasante

(1) Superficie del pavimento

16/06/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

http://training.ce.washington.edu/wsdot/Mod
ules/06_structural_design/06-4_body.htm

Ejemplos (Asphalt Institute)

25.64
=

1
= 0.414

3.291

4.196

0.854

= 1.365 10

4.477

= 1.365 10

4.477

28

Tablas de Jones
Los grficos y tablas completas estn en la referencia de Yoder & Witczak (1975).

Se definen cuatro relaciones adimensionales para el


sistema:

Se obtienen los valores de influencia para el clculo de los


esfuerzos verticales y radiales en las interfases 1 y 2:

1 =
=

1
2

2 =

2
3

1
2

1 = 1
2 = 2
1 1 = 1 1
2 2 = 2 2
2 3 = 2 3

Con los esfuerzos obtenidos y las ecuaciones constitutivas


de la elasticidad lineal se obtienen las deformaciones.

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

1
+

1
= +

29

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

30

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

31

Ejemplo:
Un pavimento se compone de tres capas con las caractersticas que se presentan a
continuacin:
Capa

Espesor (m)

Rel. Poisson

Mdulo de elasticidad (MPa)

0.10

0.5

1,400

0.20

0.5

700

Infinito

0.5

35

Se aplica una carga uniforme (q) de 600 kPa sobre un rea circular de 0.16 m de radio
(a).
Calcule:
Los esfuerzos en las interfases.
La deformacin horizontal de traccin en el fondo de la capa asfltica.
La deformacin vertical de compresin en la parte superior de la subrasante.

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

32

Solucin:
Se obtienen los parmetros adimensionales del sistema:

k1 =
k2 =
H=
A=

E1 / E2
E2 / E3
h1 / h2
a / h2

=
=
=
=

1,400 MPa / 700 MPa


700 MPa / 35 MPa
0.10 m / 0.20 m
0.16 m / 0.20 m

=2
= 20
= 0.5
= 0.8

De la tabla de factores de Jones se obtiene:

14/05/2015

ZZ1
ZZ2
ZZ1 RR1
ZZ2 RR2
ZZ2 RR3

=
=
=
=
=

0.69098
0.06476
0.86191
0.91168
0.04558

LUIS RICARDO VSQUEZ VARELA, M.Sc.

33

Esfuerzos en las interfases:


1 = 1 = 600 0.69098 = 414.6
1 1 = 1 1 = 600 0.86191 = 517.1
1 = 1 517.1 = 414.6 517.1 = 102.5
2 = 2 = 600 0.06476 = 38.86
2 2 = 2 2 = 600 0.91168 = 547
2 = 2 547 = 38.86 547 = 508.1

3 = 2 = 38.86
2 3 = 2 3 = 600 0.04558 = 27.35
3 = 2 27.35 = 38.86 27.35 = 11.51

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

34

Deformacin horizontal de traccin en el fondo de la capa asfltica.


1 =

1
1 + 1
1
1
1 400,000

1 =

102.5 0.5 414.6 102.5 = 184 106

Deformacin vertical de compresin en la parte superior de la subrasante:


3 =

1
3 + 3
3

3 =

14/05/2015

1
38.86 0.5 11.51 + 11.51 = 781 106
35,000

LUIS RICARDO VSQUEZ VARELA, M.Sc.

35

Transformacin de capas.
Las soluciones grficas y tabulares estn limitadas hasta tres capas.
Para un nmero mayor de capas, se pueden unir materiales mediante la
conversin de sus espesores a un valor equivalente de un material de
referencia con igual rigidez a la flexin.
Considere dos vigas de ancho unitario, compuestas por los materiales 1 y 2,
con la misma relacin de Poisson (1 = 2):
h1

E1

1 1 = 1

13
12

h2

b1 = 1.0

E2

2 2 = 2

23
12

b2 = 1.0

Al igualar la rigidez a la flexin de las dos vigas, se obtiene la relacin que


deben cumplir sus alturas para que su capacidad estructural sea equivalente.
14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

36

1 1 = 2 2

13
23
1
= 2
12
12

2 =

1
1
2

A partir de esta relacin, Odemark (1949) propuso una ecuacin de la forma:

1 13 = 2 23

1
= 0.9
2

Donde:
heq:

Espesor equivalente de un material con mdulo E1 al transformarlo en el


material con mdulo E2.

h:

Espesor original del material con mdulo E1.

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

37

Ejemplo:
Convierta el sistema de tres capas del ejemplo anterior (ver Tabla) en un sistema equivalente
de dos capas y obtenga el esfuerzo vertical en la subrasante mediante el baco propuesto por
Huang (1969).
Capa

Espesor (m)

Rel. Poisson

Mdulo de elasticidad (MPa)

0.10

0.5

1,400

0.20

0.5

700

Infinito

0.5

35

Haga el anlisis transformando:

Las capas 1 y 2 en el material de la capa 1 (E1 = 1,400 MPa).

Las capas 1 y 2 en el material de la capa 2 (E2 = 700 MPa).

Se aplica una carga uniforme (q) de 600 kPa sobre un rea circular de 0.16 m de radio (a).

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

38

Solucin:
Transformando las capas 1 y 2 en el material de la capa 1 (E1 = 1,400 MPa), el
espesor equivalente del pavimento es:
1 = 0.10 + 0.90

700
0.20 = 0.2429
1,400

Transformado las capas 1 y 2 en el material de la capa 2 (E2 = 700 MPa), el


espesor equivalente del pavimento es:
3

1 = 0.90

14/05/2015

1,400
0.10 + 0.20 = 0.3134
700

LUIS RICARDO VSQUEZ VARELA, M.Sc.

39

En el primer caso, el esfuerzo vertical en la subrasante es:


a / h1
E1 / E2
Relacin c / q

= 0.16 m / 0.2429 m
= 1,400 MPa / 35 MPa

= 0.66
= 40
= 0.07

= 0.07 . = 0.07 600 = 42

En el segundo caso, el esfuerzo vertical en la subrasante es:


a / h1
E1 / E2
Relacin c / q

= 0.16 m / 0.3134 m
= 700 MPa / 35 MPa

= 0.51
= 20
= 0.07

= 0.07 . = 0.07 600 = 42


Las respuestas se ajustan razonablemente con el valor de 38.86 kPa obtenido
con las tablas de coeficientes de Jones para sistemas de tres capas.
14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

40

Software disponible.
La solucin del sistema descrito en esta presentacin se conoce como
Anlisis de Capas Elsticas o LEA:
Layered Elastic Analysis.

Se ha publicado un nmero importante de programas para computadora


que emplean este sistema para el anlisis y diseo de pavimentos
flexibles.
La disponibilidad se ha reducido en los ltimos aos debido a problemas
de compatibilidad con los sistemas operativos modernos de 64 bits.
Se sugiere emplear el programa WESLEA para familiarizarse con la
aplicacin del anlisis de capas elsticas en pavimentos flexibles.

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

41

Programa

Autor

Licencia

Empleo actual

CHEVRON

Warren y Dieckman, USA

Recodificado en otros programas

BISAR

Shell Petroleum, UK

Comercial

Problemas de compatibilidad con


Windows de 64 bits

ELSYM5

Federal Highway
Administration & UC Berkeley

Comercial

Limitado (DOS)

KENPAVE

Yang H. Huang

Comercial

Extendido

JULEA

Jacob Uzan (Technion, Israel)

Comercial

Hace parte de la nueva MEPDG

LEAF

Gordon Hayhoe (FAA, USA)

Pblico

En aumento.

Extendido

WESLEA

F. J. Van Cauwelaert

Variable

Varias licencias
PerROAD
WESLEA
EVERSTRESS

Aliz 3

Variable

Limitado (DOS)

ALIZE-LCPC

Laboratoire Central des Ponts


et Chausss (Francia)

Comercial

Extendido

CIRCLY

Mincad Systems (Australia)

Comercial

Extendido

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

42

Solucin del ejemplo de tres capas en WESLEA.

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

43

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

44

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

45

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

46

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

47

Comparacin de resultados:
bacos de Jones.
Punto

Capa

Z (m)

z (kPa)

r (kPa)

r (m/m)

z (m/m)

0.10

414.6

-102.5

- 184.7

--

0.30

38.86

11.51

--

781.4

Punto

Capa

Z (m)

z (kPa)

x = y (kPa)

x = y (m/m)

z (m/m)

0.10

414.3

-103.1

-184.8

--

0.30

38.85

11.51

--

781.4

WESLEA.

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

48

Ejemplo de asignacin de propiedades


mecnicas en un pavimento flexible.

Considere una estructura de pavimento flexible compuesta por concreto asfltico,


base granular y subbase granular sobre una fundacin de suelo fino.

Para cada material se presentan las ecuaciones constitutivas y el modelo de


comportamiento correspondiente.

Se deben asignar las propiedades de cada material y, mediante un proceso


iterativo, verificar que el mdulo de elasticidad movilizados se ajusta a las
condiciones de trabajo (esfuerzo) de cada capa:
Subrasante: El mdulo elstico es funcin del nivel de esfuerzo. El punto de control queda en
la parte superior de la capa y se analiza la deformacin permanente.
Capas granulares: El mdulo elstico es funcin del nivel de esfuerzo. El punto de control
queda en la parte superior de cada capa granular.
Capa asfltica: El mdulo elstico es funcin de la temperatura de la mezcla asfltica y la
frecuencia esperada de las cargas de trnsito. El punto de control queda la parte inferior de la
capa y se analiza la fatiga de la misma.

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

49

Subrasante:
El mdulo resiliente de la subrasante se define mediante la ecuacin:

= 235.2

0.429

Donde:
Mr:
Pa:
d:

Mdulo resiliente del suelo fino de subrasante (kPa).


Presin atmosfrica (100 kPa).
Esfuerzo desviador en el punto de control de comportamiento (kPa).

El modelo de comportamiento frente a la deformacin permanente se define


mediante la ecuacin:
13,626
=

3.870

Donde:
Nd:
z:

Nmero de repeticiones admisibles de carga para desarrollar un ahuellamiento de 13


mm de profundidad.
Deformacin vertical de compresin en el punto de control de comportamiento (m/m).

El punto de control para el anlisis de subrasante se encuentra en la parte


superior de la capa, bajo el centro del eje del rea cargada.

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

50

Capas granulares:
Los mdulos resilientes de las capas de base granular y subbase granular
se definen mediante la ecuacin:

= 306.1

1.415

0.467

Donde:
Mr:
Pa:
:

Mdulo resiliente del material granular (kPa).


Presin atmosfrica (100 kPa).
Primera invariante de esfuerzos normales en el punto de control (kPa).

= + +
oct:

Esfuerzo cortante octadrico en el punto de control (kPa).

1 2

El punto de control para el anlisis de las capas granulares se encuentra en


la parte superior de cada capa, bajo el centro del eje del rea cargada.

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

51

50

Concreto asfltico:

50 mm

45

100 mm
400 mm
600 mm
35

30

T HMA (C)

La relacin entre la temperatura del aire,


el espesor y temperatura de trabajo de la
mezcla asfltica se obtendr del grfico
publicado por SHELL.

200 mm

40

25

20

El mdulo resiliente del concreto asfltico


se obtiene a partir de los valores
reportados por el USACE.

15

10

Considerando una temperatura media


mensual ponderada del aire de 20C y una
frecuencia de carga de 8 Hz se proponen
los siguientes valores de mdulo elstico:

14/05/2015

Espesor
(cm)

Temperatura
mezcla (C)

E (MPa) para
f = 8 Hz.

31.0

1,775

10

30.0

2,000

20

28.6

2,325

LUIS RICARDO VSQUEZ VARELA, M.Sc.

10

15

20

25

30

MAAT or w-MAAT (C)

52

50

50 mm

45

100 mm
200 mm

40

400 mm
600 mm

35

T HMA (C)

30

25

20

15

10

0
0

10

15

20

25

30

MAAT or w-MAAT (C)

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

53

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

54

El modelo de comportamiento frente a la


fatiga se define con la ecuacin calibrada del
Asphalt Institute (1982):

3.291

0.414
3.291

0.854

Donde:
N f:
t:
EHMA:

Nmero de repeticiones admisibles


de carga.
Deformacin horizontal de
traccin en el punto de control (m/m).
Mdulo elstico del concreto asfltico
(kPa).

Para los espesores de capa asfltica y sus


mdulos elsticos se obtienen los valores de
f1 en la tabla adjunta.
El punto de control para el anlisis del
concreto asfltico se encuentra en la parte
inferior de la capa, bajo el centro del eje del
rea cargada.

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

hHMA (cm)
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20

THMA (C)
30.9
30.8
30.6
30.4
30.3
30.1
30.0
29.8
29.7
29.5
29.3
29.2
29.0
28.9
28.7
28.6

EHMA (MPa)
1,766
1,830
1,885
1,935
1,980
2,021
2,059
2,094
2,127
2,158
2,187
2,215
2,241
2,266
2,290
2,313

f1
1.915E-06
1.858E-06
1.811E-06
1.771E-06
1.736E-06
1.706E-06
1.679E-06
1.655E-06
1.633E-06
1.613E-06
1.595E-06
1.578E-06
1.562E-06
1.547E-06
1.534E-06
1.521E-06

55

Se proponen una serie de estructuras de pavimento con las siguientes


caractersticas:
Capa

Espesor (cm)

()

E (kPa)

Peso unitario (kN/m)

Concreto asfltico
(w-MAAT = 20C, f = 8 Hz)

K0

5
()
20

0.35

1766,000
()
2313,000

24

Base granular

15

0.35

Por estimar

21

3.0

Subbase granular

45

0.35

Por estimar

21

3.0

Subrasante

Semi-infinito

0.45

Por estimar

18

0.5

El anlisis no lineal de las capas granulares y la subrasante requiere considerar el


esfuerzo geoesttico inducido por el peso propio de los materiales.

El esfuerzo geoesttico horizontal se estima mediante el coeficiente K0, el cual


vara entre 0.4 para materiales inalterados y 3.0 para materiales compactados
(Harichandran & Baladi, 2000).

El incremento de esfuerzo debido a la carga de trnsito se representa con una


huella circular de 15 cm de radio, sobre la cual se aplica una presin uniforme de
566 kPa.

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

56

Estructura con 5 centmetros de concreto asfltico.


Capa
HMA
BG
SBG
SR

Espesor
(m)

Mdulo
(MPa)

R. Poisson
()

Peso
unitario
(kN/m)

0,05
0,15
0,45

1.766
Ver tabla
Ver tabla
Ver tabla

0,35
0,35
0,35
0,45

24
21
21
18

K0
= 306.1
3,0
3,0
0,5

= 235.2

1.415

0.467

Pa = 100 kPa

0.429

Iteracin 0
zg (kPa)
1,20
4,35
13,80

xg (kPa) z (kPa) x (kPa)


3,60
13,05
6,90

z (kPa)
1,20
4,35
13,80

x (kPa)

zg (kPa)
1,20
4,35
13,80

xg (kPa) z (kPa) x (kPa) z (kPa)


3,60
199,8
78,84 201,00
13,05
142,97
28,56 147,32
6,90
48,48
6,92
62,28

x (kPa)

Capa

zg (kPa)

xg (kPa) z (kPa) x (kPa)

z (kPa)

x (kPa)

(kPa)

oct (kPa)

BG
SBG
SR

1,20
4,35
13,80

498,50
99,75
36,88

247,99
12,59
8,57

994,48
124,93

118,09
41,09

Capa
BG
SBG
SR

3,60
13,05
6,90

(kPa)
8,40
30,45

oct (kPa)
1,13
4,10

d (kPa)

Mr (MPa)

3,45

7,5
25,3
99,7

d (kPa)

Mr (MPa)

Error

24,23

251,8
138,2
43,2

3275,2%
446,4%
-56,7%

d (kPa)

Mr (MPa)

Error

14,16

689,0
63,5
54,4

0,0%
0,7%
-0,5%

Iteracin 1
Capa
BG
SBG
SR

82,44
41,61
13,82

(kPa)
365,88
230,54

oct (kPa)
55,89
49,83

Iteracin 9

14/05/2015

3,60
13,05
6,90

497,3
95,4
23,08

244,39
-0,46
1,67

LUIS RICARDO VSQUEZ VARELA, M.Sc.

57

Estructura con 10 centmetros de concreto asfltico.


Capa
HMA
BG
SBG
SR

Espesor
(m)

Mdulo
(MPa)

R. Poisson
()

Peso
unitario
(kN/m)

0.1
0.15
0.45

2,021
Ver tabla
Ver tabla
Ver tabla

0.35
0.35
0.35
0.45

24
21
21
18

K0
= 306.1
3.0
3.0
0.5

= 235.2

1.415

0.467

Pa = 100 kPa

0.429

Iteracin 0
zg (kPa)
2.40
5.55
15.00

xg (kPa) z (kPa) x (kPa)


7.20
16.65
7.50

z (kPa)
2.40
5.55
15.00

x (kPa)

zg (kPa)
2.40
5.55
15.00

xg (kPa) z (kPa) x (kPa) z (kPa)


7.20
93.68
30.29
96.08
16.65
68.45
13.67
74.00
7.50
28.27
4.24
43.27

x (kPa)

Capa

zg (kPa)

xg (kPa) z (kPa) x (kPa)

z (kPa)

x (kPa)

(kPa)

oct (kPa)

BG

2.40

7.20

212.2

17.17

214.60

24.37

263.34

SBG
SR

5.55
15.00

16.65
7.50

91.61
22.89

-2.67
0.86

97.16
37.89

13.98
8.36

125.12

Capa
BG
SBG
SR

7.20
16.65
7.50

(kPa)
16.80
38.85

oct (kPa)
2.26
5.23

d (kPa)

Mr (MPa)

3.75

14.4
31.9
54.7

d (kPa)

Mr (MPa)

Error

15.77

119.3
97.5
52.0

729.2%
206.2%
-5.0%

d (kPa)

Mr (MPa)

Error

89.68

126.8

-0.3%

39.21

65.1
53.4

1.0%
-0.9%

Iteracin 1
Capa
BG
SBG
SR

37.49
30.32
11.74

(kPa)
171.06
134.64

oct (kPa)
27.62
20.59

Iteracin 6

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

14.77

58

Estructura con 18 centmetros de concreto asfltico.


Capa
HMA
BG
SBG
SR

Espesor
(m)

Mdulo
(MPa)

R. Poisson
()

Peso
unitario
(kN/m)

0.18
0.15
0.45

2,266
Ver tabla
Ver tabla
Ver tabla

0.35
0.35
0.35
0.45

24
21
21
18

K0
= 306.1
3.0
3.0
0.5

= 235.2

1.415

0.467

Pa = 100 kPa

0.429

Iteracin 0
Capa
BG
SBG
SR

zg (kPa)
4.32
7.47
16.92

xg (kPa) z (kPa) x (kPa)


12.96
22.41
8.46

z (kPa)
4.32
7.47
16.92

x (kPa)

zg (kPa)
4.32
7.47
16.92

xg (kPa) z (kPa) x (kPa) z (kPa)


12.96
45.5
10.95
49.82
22.41
33.33
5.66
40.80
8.46
16.11
2.72
33.03

x (kPa)

zg (kPa)
4.32
7.47
16.92

xg (kPa) z (kPa) x (kPa) z (kPa)


12.96
73.88
7.59
78.20
22.41
44.82
1.85
52.29
8.46
15.43
1.22
32.35

x (kPa)

12.96
22.41
8.46

(kPa)
30.24
52.29

oct (kPa)
4.07
7.04

d (kPa)

Mr (MPa)

4.23

25.1
42.2
54.7

d (kPa)

Mr (MPa)

Error

10.93

79.0
109.0
54.7

214.5%
158.1%
0.0%

d (kPa)

Mr (MPa)

Error

11.34

72.2
79.7
54.7

-0.4%
-0.4%
0.0%

Iteracin 1
Capa
BG
SBG
SR

23.91
28.07
11.18

(kPa)
97.64
96.94

oct (kPa)
12.21
6.00

Iteracin 6
Capa
BG
SBG
SR

14/05/2015

20.55
24.26
9.68

(kPa)
119.30
100.81

LUIS RICARDO VSQUEZ VARELA, M.Sc.

oct (kPa)
27.18
13.21

59

Las mdulos de Young movilizados


para las estructuras finalmente
obtenidas son:
hHMA (cm)

EHMA (MPa)

EBG (MPa)

ESBG (MPa)

ESR (MPa)

5
6
7
7.5
8
9
10
11
12
13
14
15
16
18
20

1,766
1,830
1,885
1,911
1,935
1,980
2,021
2,059
2,094
2,127
2,158
2,187
2,215
2,266
2,313

689.0
363.3
245.8
213.1
187.2
151.5
126.8
111.3
99.9
90.8
84.4
79.7
76.3
72.2
71.0

63.5
62.8
63.0
63.4
63.6
64.5
65.1
66.1
67.2
68.3
69.9
71.7
74.0
79.7
87.5

54.4
53.2
52.8
52.7
52.8
53.0
53.4
54.1
54.7
54.7
54.7
54.7
54.7
54.7
54.7

14/05/2015

Para cada estructura se obtienen las


mximas deformaciones unitarias de
traccin (x o y) en la capa asfltica y
de compresin (z) en la subrasante,
las cuales permiten analizar los
comportamientos de fatiga y
deformacin permanente.

Debe recordarse que la estimacin


del mdulo de las capas granulares y
de la subrasante se hizo slo en los
puntos de control seleccionados.

Un mtodo como los elementos


finitos (p. e. MICHPAVE) permitira
representar la variacin del mdulo
de elasticidad en diferentes puntos
de las capas de comportamiento no
lineal.

LUIS RICARDO VSQUEZ VARELA, M.Sc.

60

La estructura de 5.0 centmetros de capa


asfltica moviliza el valor mximo definido
del mdulo de la base granular (689 MPa), lo
cual slo podra satisfacerse con materiales
de alta calidad.

En general, los mdulos movilizados de la


capa de subbase granular y de la subrasante
presentan una menor variacin que el
mdulo movilizado de la base granular.

De hecho, el mdulo movilizado de la


subrasante se mantiene en un valor cercano
al mximo (54.7 MPa) asociado con el
esfuerzo desviador mnimo de 14 kPa.

Para un espesor de capa asfltica mayor que


17 centmetros, el mdulo movilizado de la
subbase granular supera a aquel movilizado
por la base granular.

Debe recordarse que ambas capas tienen la


misma ecuacin constitutiva no lineal en este
ejercicio.

Se esperara que un material de base granular


triturado tuviese coeficientes (k1, k2, k3)
diferentes a un material sin trituracin.

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

61

Los comportamientos por fatiga y ahuellamiento son:

La estructura de peor comportamiento (mxima deformacin y mnimas repeticiones


admisibles) frente a la fatiga es aquella con 9.5 centmetros de capa asfltica.
La estructura de peor comportamiento frente al ahuellamiento es aquella con 7.0
centmetros de capa asfltica.

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

62

Conclusiones:
En el ejercicio anterior debe resaltarse que las propiedades mecnicas de las capas del
pavimento y de la subrasante variaron en diferente medida para cada uno de los espesores
propuestos de capa asfltica entre 5.0 y 20.0 centmetros.
Los factores que definen dicha variabilidad son: la temperatura, el espesor y la frecuencia de
carga en la capa asfltica y el nivel de esfuerzo en las capas granulares y la subrasante.
Es evidente que el anlisis de capas elsticas (LEA) para el diseo de pavimentos no puede
restringirse a valores constantes de los mdulos de Young de los materiales considerados en
un amplio espectro de espesores propuestos.
Asimismo, se tendr una mayor complejidad en el problema si se consideran las diferentes
magnitudes de las cargas aplicadas y la variacin de la temperatura del aire. Estas cantidades
fueron constantes en el ejercicio anterior (eje de 80 kN, frecuencia de 8 Hz y w-MAAT = 20C).
El buen diseo de pavimentos por metodologas mecanicistas
representar, dentro del mtodo de anlisis escogido, numerosas
factores de diseo considerados para establecer cmo se consumir
estructural del pavimento en el tiempo (fatiga y ahuellamiento)
serviciabilidad.

14/05/2015

LUIS RICARDO VSQUEZ VARELA, M.Sc.

empricas requiere
combinaciones de los
la oferta de capacidad
y su efecto sobre la

63

También podría gustarte