XXVI Reunin Nacional de Mecnica de Suelos

Sociedad Mexicana de e Ingeniera Geotcnica

Ingeniera Geotcnica, A.C. Noviembre 14 a 16, 2012 Cancn, Quintana Roo

Incremento de esfuerzos producido por un terrapln

Stress increment produced by an embankment load

1
Miguel MANICA

1
Estudiante de posgrado, Universidad Nacional Autnoma de Mxico.

RESUMEN: La determinacin del incremento de esfuerzo vertical en una masa de suelo debido a cargas superficiales es
una labor fundamental en el clculo de asentamientos. En el caso de cargas producidas por la construccin de un
terrapln, las soluciones disponibles consideran nicamente una porcin del mismo, y es necesaria la adecuada adicin
o sustraccin de dichas porciones para tomar en cuenta el efecto del terrapln completo. El presente artculo describe la
obtencin de las expresiones fcilmente programables que permiten la determinacin del incremento de esfuerzos
(vertical, horizontal y cortante) producido por un terrapln completo, en cualquier punto del medio. Se presenta un
ejemplo para ilustrar la aplicacin de las expresiones obtenidas, as como los resultados de un programa de
computadora generado a partir de estas expresiones.

ABSTRACT: The vertical stress increment computation within a soil mass due to surface loads is a fundamental task for
settlements calculations. In the case of embankment loads, the available solutions only consider a portion of the
embankment and a suitable addition or subtraction of such portions is needed to take account the effect of the whole
structure. This paper describes the obtainment of the easily programmable equations allowing the stress increment
computation (vertical, horizontal and shear) due to the whole embankment, at any point of the half-space. An example is
presented to illustrate the use of the equations as well as the results of a computer program built from such equations.

1 INTRODUCCIN al centro del terrapln, basta con considerar una

1.1 Antecedentes porcin que represente la mitad del mismo, y
multiplicar el esfuerzo obtenido por dos para tomar
Debido a las condiciones de los proyectos para el en cuenta el efecto del terrapln completo. Para
desarrollo de infraestructura, la construccin de puntos diferentes del centro, o incluso fuera del
terraplenes sobre suelos blandos es cada vez ms terrapln, Osterberg (1957) utiliza el principio de
comn. En este tipo de suelo, generalmente los superposicin, y mediante la correcta adicin y
estados lmite de servicio rigen el diseo debido a su sustraccin de estas porciones, obtiene el
gran compresibilidad, por lo que la estimacin de los incremento de esfuerzos verticales en cualquier
asentamientos es una tarea fundamental en la punto del medio. Adems, para la determinacin de
proyeccin de estas estructuras. No es suficiente manera expedita del aporte de cada una de las
calcular los asentamientos en puntos caractersticos porciones, presenta una carta de influencia (Figura
como en el centro y en las orillas del terrapln; es 1). El procedimiento de Osterberg, as como su carta
necesario determinarlos a lo largo de toda su de influencia, han sido extensamente utilizados,
seccin trasversal a fin de poder evaluar su incluso en la literatura moderna (Tan, 2000; Jones et
comportamiento y estar en condiciones de tomar al, 2008). Desafortunadamente, este procedimiento
decisiones respecto a procedimientos constructivos tiene el inconveniente de volverse tardado y
o viabilidad de los proyectos. laborioso si se requiere el anlisis en numerosos
En el clculo de asentamientos, es necesaria la puntos.
determinacin del incremento de esfuerzos
verticales en la masa de suelo. Diversos autores han
1.2 Objetivos
recopilado soluciones obtenidas a partir de la teora
de la elasticidad, tales como Jrgenson (1934), Gray El objetivo de este trabajo es el proporcionar las
(1936), Newmark (1940), Poulos y Davis (1974). En expresiones que puedan ser fcilmente
el caso de terraplenes, una de estas soluciones es la implementadas en una hoja de clculo, de tal forma
obtenida por Carothers, la cual considera una que, una vez introducidas las dimensiones del
porcin del terrapln limitada en uno de sus terrapln y el peso especfico del material que lo
extremos por un plano vertical. Analizando un punto conforma, nicamente sea necesario ir cambiando

SOCIEDAD MEXICANA DE INGENIERA GEOTCNICA, A.C.

 2 Distribucin de esfuerzos producida por un terrapln sobre un medio elstico

las coordenadas del punto bajo anlisis para obtener

p x z 2 z r1
el incremento de esfuerzos (vertical, horizontal y x = + + 2 (x b ) + ln (2)
cortante) producido por el terrapln completo en a r2 a r0
dicha posicin.
p z z 2
xz = (3)
a r22
donde: p = h; y = peso especfico del material que
conforma el terrapln.

Figura 1. Carta de influencia para el clculo de esfuerzos

verticales debido a una porcin de terrapln (Osterberg,
1957). Figura 2. Solucin de Carothers para la distribucin de
esfuerzos producida por una porcin de terrapln.
1.3 Alcances
Es necesario poner los argumentos angulares y
Se describe el desarrollo que permiti obtener las
en funcin de las coordenadas del punto bajo
expresiones antes mencionadas, se resuelve un
anlisis y de la geometra del terrapln. Se considera
ejemplo para ilustrar la aplicacin de stas, y se
el tringulo que forman los puntos ADB de la Figura
presentan los resultados de un programa de
2, tal como se muestra en la Figura 3, donde:
computadora desarrollado en el lenguaje de
programacin FORTRAN, con el fin de demostrar la bx
versatilidad de las soluciones obtenidas. = tan 1 (4)
z
Adicionalmente se proporciona un enlace para
descargar una hoja de clculo en EXCEL, donde se x a
encuentran implementadas las ecuaciones = tan 1 (5)
presentadas. z

2 DETERMINACIN DE LAS EXPRESIONES

Se parte de las ecuaciones obtenidas por Carothers

(Gray, 1936) para el clculo del incremento de
esfuerzos (vertical, horizontal y cortante) en el caso
de deformacin plana mostrado en la Figura 2, en
donde es analizada una porcin de un terrapln
homogneo, sobre un medio semi-infinito,
homogneo, istropo y linealmente elstico
(ecuaciones 1, 2 y 3).

p x z
z = + 2 (x b ) (1)
a r2

Figura 3. Determinacin del argumento angular , y de las

distancias r1 y r2.

SOCIEDAD MEXICANA DE INGENIERA GEOTCNICA, A.C.

 MANICA M. 3

Por lo tanto: xB = 2 b x (20)

bx
' = 90 = 90 tan 1 (6)
z

x a
' = 90 = 90 tan 1 (7)
z

+ ' +' = 180 (8)

bx x a
+ 90 tan 1 + 90 tan 1 = 180 (9)
z z

bx x a
= tan 1 + tan 1 (10)
z z
Se ponen tambin las distancias r2 y r1 en funcin Figura 4. Ensamble de las dos porciones de un terrapln
de las coordenadas del punto bajo anlisis y de la
geometra del terrapln. Por conveniencia r2 De lo anterior, se obtienen los valores de , , r0 y
permanece elevado al cuadrado. r1 de ambos sistemas en funcin de las variables de
inters (ecuaciones 21 a la 28). En la Figura 4 puede
r22 = (b x )2 + z 2 (11) observarse que el valor de r2 es el mismo en ambos
sistemas, por lo que no se hace distincin de este
1/ 2 parmetro para cada uno de ellos (Ecuacin 29).
[
r1 = (x a)2 + z 2 ] (12) Los argumentos angulares de las ecuaciones (21) a
la (24) se obtienen directamente en radianes.
De forma similar se determina y r0 analizando el
tringulo formado por BDC de la Figura 2, con lo que bx x a
A = tan 1 + tan 1 (21)
se obtiene: z z 180
ax x
= tan 1 + tan 1 (13) x b 2 b x a
z z B = tan 1 + tan 1 (22)
z z 180
(
r0 = x 2 + z 2 ) 1/ 2
(14)
ax x

Ahora todas las variables de las expresiones (1), A = tan 1 + tan 1 (23)
z z 180
(2) y (3) estn en funcin de la geometra del
terrapln (a, b, h), del peso especfico del material
a 2b + x 2 b x
que lo conforma () y de la posicin del punto bajo B = tan 1 + tan 1 (24)
anlisis (x, z), pero todava se analiza nicamente z z 180
una porcin del terrapln.
Considerando un terrapln con simetra, en la
Figura 4 se ensamblan dos sistemas (sistema A y B)
(
r0 A = x 2 + z 2 ) 1/ 2
(25)
como el mostrado en la Figura 2.
1/ 2
Por simetra:
[
r0 B = (2 b x )2 + z 2 ] (26)
bA = bB = b (15)
1/ 2
aA = aB = a (16) [
r1A = (x a )2 + z 2 ] (27)

1/ 2
hA = hB = h (17)
[
r1B = (2b x a)2 + z 2 ] (28)

zA = zB = z (18)
r22 = (b x )2 + z 2 (29)
Para tener ambos sistemas en funcin de las
mismas variables se considera lo siguiente: Finalmente, basta con sumar el aporte de cada
uno de los sistemas a partir de las expresiones (1),
xA = x (19) (2) y (3) tal como se hace en las ecuaciones (30),
(31) y (32), las cuales estn referidas a la Figura 5.

SOCIEDAD MEXICANA DE INGENIERA GEOTCNICA, A.C.

 4 Distribucin de esfuerzos producida por un terrapln sobre un medio elstico

x A z
A + 2 (x b ) +
h a r2
z = (30)
(2 b x ) B z
B + 2 (b x )

a r2

x A z 2 z r1 A
A + + 2 (x b ) + ln +
a r2 a r0 A
h (2 b x ) B

x = B + (31)
a Figura 6. Geometra del terrapln y ubicacin del punto
z 2 z r1B
bajo anlisis.
+ 2 (b x ) + ln

r2 a r0 B
1 9.15 12.2
tan
h z 2 z 2 B = 6 . 1 = 0.519 rad
xz = ( A + B ) 2 (32)
a r2 1 2 (12.2 ) 9.15 6.1 180
+ tan
6.1

6.1 9.15
tan 1
6 .1 = 0.519 rad
A =
1 9.15
180
+ tan
6.1

6.1 2 (12.2 ) + 9.15

tan 1 +
6.1 = 0.207 rad
B =
1 2 (12.2 ) 9.15
180
tan
6.10

2
r22 = (12.2 9.15 ) + 6.10 2 = 46.513 m 2

Figura 5. Variables necesarias para el clculo del

Finalmente dichos valores son sustituidos en la
incremento de esfuerzos debajo de un terrapln. expresin (30) para determinar el incremento de
esfuerzo vertical:
Con las expresiones (21) a la (32) es posible
determinar las tres componentes de esfuerzo (z, x,
9.15 (0.519 )
0.927 +
xz) producidas por el terrapln completo, 6 . 1

nicamente en funcin de su geometra, el material 6.1
que lo conforma y la ubicacin del punto bajo (9 . 15 12. 2 )
18 (5 ) 46.513
anlisis. z =

0 . 519 +
(2 (12.2 ) 9.15 )0.207
6 .1
3 EJEMPLO +
6.1 (12.2 9.15 )
A partir del terrapln mostrado en la Figura 6, se 46.513
determina la magnitud del incremento de esfuerzo z = 78.6 kPa
vertical (z) en el punto A.
2
Se calculan los valores de A, B, A, B y r2 con las
expresiones (21), (22), (23), (24) y (29)
respectivamente. 4 PROGRAMACIN DE LAS EXPRESIONES

1 12.2 9.15 Con el fin de mostrar la versatilidad de las

tan expresiones obtenidas, se gener una aplicacin en
6 .1 = 0.927 rad
A = el lenguaje de programacin FORTRAN, la cual
1 9.15 6.1 180
+ tan
6 .1

SOCIEDAD MEXICANA DE INGENIERA GEOTCNICA, A.C.

 MANICA M. 5

permite la determinacin del incremento de esfuerzo no puede tomar valores negativos, sin embargo, el
vertical en mltiples puntos simultneamente. valor de x s puede hacerlo para tomar en cuenta
Se tom el ejemplo presentado en la Figura 6 y se los puntos que estn a la izquierda del terrapln.
consider como origen del sistema coordenado el
extremo izquierdo del terrapln (tal como debe de
hacerse al utilizar las expresiones deducidas en la 6 CONLUSIONES
seccin 2). Se determin el incremento de esfuerzo
vertical desde x = -20.0 hasta x = 45.0 y desde z = A pesar de que a simple vista las expresiones aqu
0.0 hasta z = 40.0 metros con separaciones en presentadas puedan parecer de tediosa aplicacin,
ambas direcciones de un metro. As, el programa es importante hacer notar que una vez
determin el incremento de esfuerzo vertical en un implementadas en una hoja de clculo o en
poco ms de 2,700 puntos. Los resultados fueron cualquier lenguaje de programacin, se dispondr
procesados con un software para construir mapas de de una herramienta que permita la rpida
contornos. La Figura 7 muestra los resultados determinacin del incremento de esfuerzos en
obtenidos. En dicha figura tambin aparece indicada cualquier punto del medio, facilitando as el clculo
la ubicacin del punto que fue analizado en el de asentamientos.
ejemplo de la seccin 3, donde puede observarse
que en ambos clculos el resultado obtenido es el
REFERENCIAS
mismo.
Las ecuaciones obtenidas en el presente artculo, Gray, H. (1936). Stress distribution in elastic solids,
tambin fueron implementadas en una hoja de
Proceedings of the 1st International Conference
clculo de EXCEL, la cual puede descargarse en el
on Soil Mechanics and Foundation Engineering,
siguiente enlace:
Vol. 2, pp. 157-168.
http://miguelmanica.webs.com/descargas
Jones, C. A.; Stewart, D. I. and Danilewicz, C. J.
(2008). Bridge distress caused by approach
embankment settlement. Geotechnical
Engineering, 161 (2), pp. 63-74.
Jrgenson, L. (1934). The applications of elasticity
and plasticity to foundation problems, Journal of
Boston Society of Civil Engineers, 21, 206.
Newmark, N. M. (1940). Stress distribution in soils,
Proceedings of the Purdue Conference on Soil
Mechanics and its Applications, pp. 295-303.
Osterberg J. O. (1957) Influence values for vertical
stresses in a semi-infinite mass due to
embankment loading. Proceedings of the 4th
International Conference on Soil Mechanics and
Foundation Engineering, Vol. 1, pp. 393394.
Poulos, H. G. and Davis, E. H. (1974). Elastic
Solutions for Soil and Rock Mechanics, 1st
Edition, New York, John Wiley & Sons, Inc.
Tan, Y.C. & Gue, S.S. (2000) Embankment over
Figura 7. Mapa de contorno de la distribucin del Soft Clay Design and Construction Control,
incremento de esfuerzos verticales debajo de un terrapln. Seminar on Geotechnical Engineering, September
2000, Penang.
5 COMENTARIOS

Es importante notar que en las ecuaciones (21) a la

(24), el valor de z no puede ser igual a cero, ya que
conducira a una indeterminacin. Esto se resuelve
simplemente asignndole a z un valor muy
-5
pequeo cercano a cero tal como 10 m. De igual
forma a igual a cero conducira a una
indeterminacin en las ecuaciones (30), (31) y (32).
Aunque este problema puede ser resuelto de la
misma forma que para z igual a cero, a igual a
cero implicara una carga rectangular, por lo que no
tendra sentido el utilizar las soluciones aqu
presentadas. Tambin es importante resaltar que z

SOCIEDAD MEXICANA DE INGENIERA GEOTCNICA, A.C.