Robert Hooke (1635 1703)

Tomada de Gordon, J.E.

INTRODUCCIN A LA ELASTICIDAD PARA SUELOS

Luis Ortuo
Thomas Young (1773-1829)

ELASTICIDAD
El modelo ms simple para un suelo o una roca sera un material de CHILE (Continuo, Homogneo, Istropo, Lineal y Elstico). Slo requiere 2 parmetros: E, . La historia de su desarrollo desarrollo, tiene su gracia (ver Gordon Gordon, J J.E. E 1987):

1.- Robert Hooke (1635-1703) su sobrina, su amigo el relojero y. Newton 1676: Hooke publica La La verdadera teora de la elasticidad . Como subttulo, un aviso sin falsa modestia: Una dcima de una centsima parte de las invenciones que pretendo publicar publicar , y un jeroglfico: CEIIINOSSSTTUV ESTABLECE LA PROPORCIONALIDAD ENTRE FUERZA Y DESPLAZAMIENTO 1679: Posteriormente, en De De potentia restitutiva, or of a spring, resuelve el jeroglfico:

UT TENSIO SIC VIS (como la extensin, as la fuerza),

Ntese que no existe por el momento el concepto de tensin. Hay que esperar unos cuantos aos) Luis Ortuo

ELASTICIDAD 2.- Mdulo de elasticidad. Thomas Young (1773-1829) y los jeroglficos egipcios gp Young, que entre otras cosas dedic parte de su vida al estudio de los jeroglficos egipcios, defini su mdulo (E), (E) la constante constante de proporcionalidad entre tensiones y deformaciones, antes de que existieran los conceptos t de d tensin t i y d deformacin. f i P Por eso quizs no le entendieron muy bien en su momento: The modulus of elasticity of any substance is a column of the same substance, , capable p of producing a force on its base which is to the weight causing a certain degree of compression as the length of the substance is to the diminution of its length.
Tomados de Gordon, J.E.

E=
Luis Ortuo

ELASTICIDAD Material Goma Cartlago (humano) Tendn (humano) Polietileno, Nylon Madera (laminada) Madera (segn la fibra) Hueso (fresco) Hormign / Concreto Aleaciones de Mg E[1]
[2]

[ MPa ] 7

Algunos valores de E (internet) Material E[1] [2] [ MPa ] 50 000 Granito vidrio 70 000 Aleaciones de 70 000 Al Latn 105 000 110 000 Bronce Hierro colado < 175 000 190 000 Hierro forjado Acero 210 000 Zafiro 420 000

24 600 1400 7000 14 000 21000 27 000 42 000

Luis Ortuo

ELASTICIDAD

3.- Tensin, deformacin y coeficiente de Poisson

Si se aplica li a una probeta b t d de un d determinado t i d material elstico un incremento de tensiones en una direccin, por la ley de Hooke se tendr un incremento de deformacin en la misma direccin, de valor:

v =

v
E

v =

L L

Ahora bien, tambin se observa que, en general, los incrementos

de tensin en una direccin producen deformaciones en perpendicular a dicha direccin.

h = v =

Donde se denomina coeficiente de Poisson

Luis Ortuo

ELASTICIDAD
Tensin, deformacin y coeficiente de Poisson El efecto Poisson es de sobra conocido (y cumple una funcin fundamental en muchos materiales y estructuras). En palabras simples, es el efecto por el cul , en general un cuerpo adelgaza cuando es estirado o engorda al ser comprimido.

v =

v
E

h =

v
Luis Ortuo

ELASTICIDAD Algunos valores del coeficiente de Poisson. Se pueden observar aspectos curiosos: El corcho tiene =0 0 (bueno para estar abriendo y cerrando botellas) A la arcilla saturada saturada, Wikipedia le asigna = 0,5. Esto, para un material elstico, equivale a decir que no cambia bi d de volumen l ( (ver ms adelante) d l t ) Hay materiales con coeficiente de Poisson negativo. Se llaman auxticos y pueden ser muy tiles

material foam glass concrete sand cast iron steel stainless steel clay copper aluminium-alloy titanium magnesium saturated clay gold auxetics cork rubber
Fuente: Wikipedia

poisson's ratio 0.10 to 0.40 0.18-0.3 0.20 0.20-0.45 0.21-0.26 0.27-0.30 0.30-0.31 0.30-0.45 0.33 0.33 0.34 0.35 0.40-0.50 0.42 negative ~ 0.00 ~ 0.50
Luis Ortuo

ELASTICIDAD Materiales auxticos (figuras tomadas de internet) engordan al estirarse

Luis Ortuo

ELASTICIDAD 4.- Tensiones normales. Ley de Hooke generalizada.

Para un material istropo, suponiendo un incremento general de tensiones normales en tres planos perpendiculares (y aplicando el concepto de la ley de Hooke y el efecto Poisson, se tiene:

( ' y + ' z ) E' ' y = ( ' x + ' z ) E' E' 'z ' ( ' y + ' x ) z = E' E' E' 'y

x =

'x

'

OBSERVACIONES: Como tratamos de suelos, empleamos tensiones efectivas (tambin le podemos poner un apstrofe a E y a , para indicar que se trata de los parmetros de tensiones efectivas). En suelos y en rocas, las compresiones son positivas. Obsrvese que en un material elstico las deformaciones normales slo dependen Luis Ortuo de la tensiones normales.

ELASTICIDAD Tensiones normales. Ley de Hooke generalizada.

x=

' x
E' ' y

E' ' y = ( ' x + ' z ) E' E' ' z ' z = ( ' y + ' x ) E' E'

'

( ' y + 'z )

La deformacin volumtrica ser:

V = v = x + y + z V 1 2 ' [ ' x + ' y + ' z ] v = E'

Teniendo en cuenta q que el incremento de p presin efectiva media es Resulta:

p' =

' x + ' y + ' z

3

v =

3(1 2 ' ) p' E'

Luis Ortuo

ELASTICIDAD
OBSERVACIONES:

V 3(1 2 ' ) = v = p' V0 E'

V = 0 p' = 0 V0

Un suelo elstico cualquiera tendr en general un mdulo E y un coeficiente de Poisson . Ahora bien, si durante un proceso de carga no hay cambio de volumen, tampoco habr cambios en la tensin efectiva media. media Lo anterior es consecuencia de asumir que el suelo es elstico. Un definicin ms: El mdulo de compresibilidad volumtrica.

p' = Kv;

E' K= 3(1 2 ' )

Luis Ortuo

ELASTICIDAD Recordatorio de tensiones tangenciales

xy =

dv du + dx dy

El convenio de signos es que tensiones positivas originen deformaciones positivas, de manera que para suelos y rocas, al elegir compresiones positivas resulta:
l f l0 l0 l du = ( p. e.) l dx

Def normales : =

du dv Def . tan genciales : = dy + dx

Luis Ortuo

ELASTICIDAD Tensiones tangenciales

mx = 1 3
mx = 1 3=

1 =

( '2 + '3 ) E' E' '3 ' 3 = ( '1 + '2 ) E' E'

'1

'

1 + ' 1 + ' [ '1 '3 ] = 2 mx E' E'

mx =

2(1 + ' ) mx = mx E' G

Luis Ortuo

ELASTICIDAD Tensiones tangenciales en materiales CHILE

mx =

2(1 + ' ) mx = mx E' G

mx = G mx

G=
OBSERVACIONES:

E' 2(1 + ' )

G es el mdulo de rigidez transversal o de corte. Las deformaciones tangenciales dependen de los incrementos de tensiones tangenciales. Un incremento de tensin tangencial no produce cambio de volumen.
Luis Ortuo

ELASTICIDAD RESUMEN DE ELASTICIDAD ISTROPA

E' ' y = ( ' x + ' z ) E' E' ' z ' z = ( ' y + ' x ) E' E'
1 E' x ' y E ' ' z xy = E ' 2 0 yz 2 0 xz 2 0

x =

' x
E' ' y

'

( ' y + ' z )

xy = yz xz

2(1 + ' ) xy E' 2(1 + ' ) = yz E' 2(1 + ' ) = xz E'

0 x y 0 z 0 xy yz 0 xz 1 + ' E' 0

'

' '
E'

E' 1 E' ' E' 0 0 0

0 0 0 1 + ' E' 0 0

0 0 0 0 1 + ' E' 0

E' 1 E' 0 0 0

xy
2

= xy

Luis Ortuo

ELASTICIDAD Condiciones edomtricas o de deformacin lateral nula en un suelo elstico e istropo.

Tomando la vertical como direccin de 1 y de 1 : En condiciones edomtricas o unidimensionales resulta:

2= 3=0, 2 =3

3 =

' 3 ' ( '3 + '1 ) = 0 E' E' ' 3 ' (1 ' ) = '1 E' E'

' ' '1 = 3 1 '

Luis Ortuo

ELASTICIDAD Condiciones edomtricas o de deformacin lateral nula en un suelo elstico e istropo.

En definitiva, suponiendo que la presin vertical efectiva es la principal mayor mayor, y la horizontal la principal menor (suelo normalmente consolidado o ligeramente sobreconsolidado), el hacer la hiptesis de que el suelo es elstico obligatoriamente implica que el coeficiente de empuje al reposo es:

'H =

' ' 'V K0 = 1 ' 1 '

Luis Ortuo

ELASTICIDAD Mdulo de compresibilidad volumtrica (mv) y mdulo edomtrico Em en un suelo elstico e istropo.
Recordando definiciones y resultados previos:

Em =

'1 1

mv =

1 '1

' 2 = ' 3 =

' '1 1 '

'1 ' 1 = ( '3 + '3 ) E' E' '1 ' ' 1 = 2 '1 E E' 1 ' '1 2 '2 1 = 1 1 ' E'

1 1 2 '2 mv = = 1 '1 E ' 1 '

Em = '1 1 ' = E' 2 1 1 '2 '
Luis Ortuo

ELASTICIDAD Mdulo de deformacin E y coeficiente de Poisson deducidos de un triaxial con drenaje (CD) para un suelo elstico
En un triaxial CD, durante la aplicacin del desviador se mide la deformacin vertical (1) y la volumtrica (V/V0). Adems, en la fase de rotura se mantiene la presin efectiva de cmara constante:

'1 > 0 ' 2 = ' 3 = 0

1 =

'1 ' ( ' 2 + ' 3 ) '1 E' E' E ' = '1 1 1 = E'
Luis Ortuo

ELASTICIDAD Mdulo de deformacin E y coeficiente de Poisson deducidos de un triaxial con drenaje (CD) para un suelo elstico

'1 ' '1 1 = [ '2 + '3 ] = E' E' E' V '3 ' ' V0 2 = 3 = [ '1 + '3 ] = '1 = 1 2 ' E' E' E' 1 V '1 = (1 2 ' ) V0 E' V = 1 (1 2 ' ) V0 V = 1 + 2 + 3 V0

Luis Ortuo

ELASTICIDAD Mdulo de deformacin Eu y coeficiente de Poisson u deducidos de un triaxial sin drenaje para un suelo elstico
En un triaxial sin drenaje, durante la aplicacin del desviador se mide (1) y (u). De la curva (1 , 1 ) se puede d deducir d d i un mdulo d l de d d deformacin f i sin i d drenaje j Eu. (Ojo, que est definido en tensiones totales y no es E)

Eu =

1 1

Por otra parte, si la saturacin es completa, en un ensayo sin drenaje no hay cambio de volumen, y por lo tanto el coeficiente de Poisson en tensiones totales asociado a procesos sin drenaje, u , ser (Ojo, que no es ):

V 3(1 2 u ) = p = 0 V0 Eu

u =

1 2
Luis Ortuo

ELASTICIDAD Relacin entre mdulos E y , con y sin drenaje en suelo elstico e istropo
Como el agua no soporta tensiones tangenciales, el mdulo de corte del terreno ha de ser igual en condiciones con o sin drenaje:

G' =

1 E' Eu = Gu = ; con u = 2(1 + ' ) 2(1 + u ) 2

Eu = E '

2(1 + u ) 2(1 + ' )

Eu =

3 E' 2(1 + ' )

De nuevo, el resultado anterior es consecuencia directa de suponer el suelo elstico e istropo

Luis Ortuo

ELASTICIDAD Parmetros de presin intersticial en suelo elstico e istropo

u = B[ 3 + A( 1 3 )]
B=1 si saturacin completa y fluido incompresible:

u = 3 + A( 1 3 )

v = 0 p' = 0

( 1 u ) + ( 2 u ) + ( 3 u ) '1 + '2 + '3 =0 = p u = 0 3 3

p = u
Triaxial (2= 3):

u =

( 1 + 2 3 ) 1 = 3 + ( 1 3 ) 3 3

A=

1 3
Luis Ortuo

EL SOLIDO ELASTICO
ANISOTROPA (las propiedades y la respuesta dependen de la direccin)
Los tejidos, tejidos los vestidos ceidos y y. una ingeniossima introduccin del Prof Prof. Gordon (1987) a la anisotropa: El corte al bies, inventado por Mme Madeleine Vionnet) anisotropa Trama y urdimbre en vertical y horizontal. Hilos poco extensibles bajo el peso propio del tejido E elevado. Poco contraccin lateral pequeo Se necesitan lazos para ceir la ropa Trama y urdimbre a 45. Extensin mayor bajo el peso propio del tejido E pequeo Gran contraccin lateral grande La ropa se cie sola y no hacen falta lazos
Tomada de Gordon, J.E. , 1987

Luis Ortuo

EL SOLIDO ELASTICO
LA ANISOTROPA (las propiedades y la respuesta dependen de la direccin)
G&C II: La anisotropa tiene mucho inters para el estudio del suelo. En especial la anisotropa t transversal, l d dada d l la f forma de d d deposicin i i d de un suelo, l con l la di direccin i vertical ti l preferencial f i ld debido bid a la gravedad, la forma laminar de las partculas ms finas del suelo, etc).

Se necesitan 5 parmetros: Ev Eh vh (efecto de v sobre h) hh (efecto de h sobre h) (mod. Rigidez g transv. Planos verticales Gvh (

Habra un hv para el efecto de h sobre v, pero por simetra:

Ghh =

E 'h 2(1 + 'hh )

'vh h
E 'v

'hv h
E 'h
Luis Ortuo

EL SOLIDO ELASTICO
LA ANISOTROPA (las propiedades y la respuesta dependen de la direccin)
n' = E 'h E 'v Gvh E 'v

Ghh =

E 'h 2(1 + 'hh )

m' =

v = h 1 = h 2

' 'v 'hv 'h1 hv 'h 2 E 'v E 'h E 'h

'h1 'vh ' 'v hh 'h 2 E 'h E 'v E 'h

hh = vh = hv =

hh Ghh vh Gvh vh Gvh

'h 2 'vh ' = 'v hh 'h1 E 'h E 'v E 'h

Luis Ortuo

Ejemplo a resolver por los alumnos

Deducir qu se determina en un triaxial CD sobre un material anistropo transversal orientado como muestra la figura:
Solucin: En un triaxial CD, durante la aplicacin del desviador se mide la deformacin vertical (v) y la volumtrica (V/V0). Adems, en la fase de rotura se mantiene la presin efectiva de cmara constante: t t

' v > 0

' h 1 = ' h 2 = 0

Luis Ortuo

Ejemplo a resolver por los alumnos

Deducir qu se mide en un triaxial CD sobre un material anistropo transversal orientado como muestra la figura:
Solucin: En el ensayo, durante la aplicacin del desviador (h1) se mide la deformacin (h1) y la volumtrica (V/V0). Adems, en la fase de rotura se mantiene la presin efectiva de cmara constante: t t

'h1 > 0

'v = 'h 2 = 0

Luis Ortuo

EL SOLIDO ELASTICO Es aplicable el modelo elstico al suelo?

En un material elstico, en un ciclo cerrado de carga-descarga se recupera la forma y no se disipa energa. No se producen deformaciones irrecuperables (plsticas). ELASTICO

Lineal

No lineal

ELASTO-PLSTICO

Luis Ortuo

EL SOLIDO ELSTICO Es aplicable el modelo elstico al suelo?

La Mecnica del Suelo y de las Rocas hace uso de la teora de la Elasticidad para el clculo de los incrementos de tensin y la estimacin de deformaciones y asientos en una multitud de situaciones (cargas rgidas o flexibles con diferentes formas o reas, con terrenos estratificados de distinta g , excavaciones, , etc). ) rigidez, Desde el punto de vista de la estimacin de incrementos de tensin en el terreno (producidos por cimentaciones, por ejemplo), Burland, Broms & De Mello (1977) indican que, dado que la mayora de los suelos no se ajustan j p precisamente a un modelo elstico de comportamiento, p el empleo p de formulaciones elsticas puede producir una cierta sensacin inicial de malestar, de que algo no se estaba haciendo bien. que g gracias al desarrollo de p potentes herramientas de No obstante, los autores citados resaltan q clculo numrico se ha podido investigar con bastante precisin el grado de error que supone el empleo de modelos elsticos (asumiendo cualquier modelo o relacin tensin-deformacin complicado como si fuera la realidad, y comparndolo con lo que se deducira de un modelo elstico simple) simple). Las conclusiones fundamentales alcanzadas en este tipo de estudio tienen consecuencias prcticas de importancia y son las siguientes:

Luis Ortuo

EL SOLIDO ELSTICO Es aplicable el modelo elstico al suelo?

1.- El modelo ms simple de semiespacio elstico de Boussinesq (suelo lineal, elstico, istropo y homogneo) proporciona en la mayora de las ocasiones rdenes de magnitud razonables para las tensiones verticales originadas por una cimentacin. cimentacin

2.- Desde un p punto de vista p prctico, una excepcin p al p prrafo anterior se p produce cuando existe una capa rgida sobre un estrato ms blando. En estas circunstancias las tensiones derivadas de la teora de Boussinesq resultan conservadoras.

3.- Al contrario de lo que ocurre con las tensiones verticales, las tensiones horizontales y las tensiones cortantes que resultan del empleo del modelo elstico de Boussinesq estn sujetas a mayores errores.

Luis Ortuo

EL SOLIDO ELASTICO Es aplicable el modelo elstico al suelo?

En lo que respecta a las deformaciones, como se recordar, al menos para pequeos ciclos de descargarecarga, las curvas o ramas de descarga y recarga se aproximan bastante. Puede ser por tanto razonable suponer que coinciden y considerar elasticidad en un suelo OC.

Idealizacin de material elasto-plstico tipo strain hardening

Ensayo de compresin istropa en triaxial (tomado de Burland, J.B. (1987)

Tomada de Gonzlez Vallejo, L. et al (2000)

Luis Ortuo

EL SOLIDO ELASTICO Es aplicable el modelo elstico al suelo?

Por otra parte, Simons & Menzies (2005) indican que en bastantes ocasiones se puede asumir comportamiento elstico si los niveles de deformacin son suficientemente pequeos, lo que puede ocurrir en condiciones de carga con coeficientes de seguridad elevados, muy alejados de la rotura. (Estas consideraciones han dado lugar a un gran desarrollo de los mtodos geofsicos en los ltimos aos. Comentar). Como ejemplo sealan que la hiptesis anterior puede ser razonable para cimentaciones con coeficientes de seguridad frente al hundimiento habituales ,del orden de 3, pero no lo sera para terraplenes, en donde los coeficientes de seguridad empleados habitualmente rondan la cifra de 1,5.
Tomadas de FHWA-IF-02-034 (2002) ( )

En cualquier caso, cuando se vayan a emplear, los parmetros parmetros elsticos elsticos deben determinarse bajo condiciones de ensayo que simulen tanto el rango de tensiones como el nivel y tipo de deformaciones esperables en el problema real a analizar. Luis Ortuo

EL SOLIDO ELSTICO Modelos elsticos para el suelo.

El modelo bsico de semiespacio de Boussinesq se puede modificar para tener en cuenta varios efectos: Existencia de una capa rgida (sustrato de roca indeformable) Semiespacio de Boussinesq. Ez=E0

Capa elstica sobre suelo rgido Ez=E0

Aumento de la rigidez del suelo con la profundidad

Tomadas de G&C, II

Funcin creciente Ez=E(z)

Crecimiento lineal Ez=E0+mz

Luis Ortuo

EL SOLIDO ELSTICO Modelos elsticos para el suelo.

G&C, II: Modelo de Winkler: Se define como un medio en el que los desplazamientos verticales S de los puntos de la superficie bajo una presin p, son proporcionales a la presin, con un coeficiente de proporcionalidad K, llamado coeficiente de balasto o mdulo de reaccin vertical.

Tomado de G&C, II

El coeficiente de balasto tiene dimensiones de peso especfico. En la figura se ve que una cimentacin bajo una presin uniforme p se hunde en el semiespacio de Winkler como si estuviese flotando en un lquido de densidad K. Ntese que en este modelo slo asientan los puntos que reciben carga. Equivale a un suelo con E creciente linealmente (E=z), E=0 en superficie y coeficiente de Poisson = 0,5.too much! (demostracin en G&C, II). Es un modelo muy usado para clculo de estructuras (una concesin, casi capitulacin, de la Geotecnia a los calculistas de estructuras, para que se pueda Luis Ortuo sustituir el suelo por un muelle).

EL SOLIDO ELSTICO El semiespacio elstico de Boussinesq. Carga vertical puntual en superficie (Boussinesq, 1885).

Este fue el problema original estudiado por Boussinesq. Boussinesq Observaciones: La tensin vertical no depende de los parmetros elsticos del terreno E, En planos horizontales, la resultante de las tensiones pasa por el punto de aplicacin de la carga

rz = tan z
Tomadas de G&C, II

Luis Ortuo

EL SOLIDO ELSTICO El semiespacio elstico de Boussinesq. Carga vertical puntual en superficie (Boussinesq, 1885).
Otras figuras y detalles tomados de G&C II.

Tomadas de G&C, II

Luis Ortuo

EL SOLIDO ELSTICO El semiespacio elstico de Boussinesq. Carga vertical puntual en superficie (Boussinesq, 1885).
Asiento:

s=

Q (1 + ) 2(1 ) + cos2 2 E

En superficie (cos=0):

Q (1 2 ) s= E
Para 0 el asiento tiende a , lo que no corresponde a la realidad (G&C II recomienda no usar la expresin para puntos cercanos al de aplicacin de la carga)

Asiento en la vertical (cos=1):

s=

Q (1 + )(3 2 ) 2 E

Tomadas de G&C, II

Luis Ortuo

EL SOLIDO ELSTICO El semiespacio elstico de Boussinesq. Una curiosidad. Reciprocidad

En un d la carga g es: q qd

Q (1 2 ) s1 = E

(1 2 ) (1 2 ) s2 = qd = 2q E E 0

Una carga lineal circular produce un asiento en el centro que depende de su intensidad q pero no del radio de la circunferencia cargada. Si los asientos son iguales, la carga total es la misma:

Q (1 2 ) 2q(1 2 ) = Q = 2q E E
Luis Ortuo

EL SOLIDO ELSTICO El semiespacio elstico de Boussinesq. Carga vertical en el interior del semiespacio p de Boussinesq q (Mindlin, ( 1936)
G&C, , II: El problema p de Mindlin tiene su aplicacin ms extendida en el estudio de pilotes, que transmiten carga al terreno a todo lo largo del fuste. En el grfico (zona izquierda) se muestran las tensiones verticales cuando la carga se supone con distribucin uniforme en la vertical, desde la superficie a la profundidad c. Ntese no obstante Nt b t t que, en la l realidad, lid d l la distribucin de carga en un pilote no es uniforme con la profundidad.

Tomada de G&C, II

Luis Ortuo

EL SOLIDO ELSTICO El semiespacio elstico de Boussinesq. Carga vertical lineal en superficie

Ntese Nt que z no depende d d d de los parmetros elsticos (E, ). G&C, G&C II indica que este caso puede asimilarse a veces a una cimentacin corrida con una zapata estrecha estrecha.

Tomada de G&C G&C, II

Luis Ortuo

EL SOLIDO ELSTICO El semiespacio elstico de Boussinesq. CARGAS RGIDAS Y FLEXIBLES

G&., II: la distribucin de cargas depende de la interaccin tenso-deformacional entre la estructura total transmisora de la carga con su cimentacin y el terreno. Por ejemplo, la distribucin de presiones bajo la zapata de un pilar depende de la rigidez de la zapata frente a la rigidez del terreno, y no tiene por qu ser uniforme.

Carga flexible
Supngase que sobre una zapata infinitamente flexible, apoyada directamente sobre la superficie de un terreno horizontal, se coloca una sobrecarga uniforme. Por efecto de esta sobrecarga, el terreno y la zapata sufrirn un asiento. El asiento ser mayor en el centro que en los extremos y no se limitar al rea cargada, sino que se extender a ambos lados hasta una cierta distancia. Si la zapata es infinitamente flexible, no ser capaz de soportar momentos flectores, y en consecuencia la distribucin de presiones con que el terreno reaccione, ser idntica a la distribucin uniforme de presiones colocada sobre la zapata. Luis Ortuo

Figuras y texto procedentes de los borradores del C.T.E.

EL SOLIDO ELSTICO El semiespacio elstico de Boussinesq. Ejemplo de asientos de cargas flexibles.

Como ejemplo, j p , suponiendo p una carga g circular uniforme y flexible de intensidad p y radio R, situada en la superficie de un semiespacio elstico de parmetros E y ,se tiene: Asiento A i t bajo b j el l centro: t

Asiento A i t bajo b j el lb borde: d

1 2 s0 = 2 pR R E 2 sR = s0

Asiento medio:

16 pR R 1 2 sm = 0,85 s0 3 E

( Ntese que el asiento elstico es proporcional R)

Casos de carga flexible: depsitos (sobre todo los de base metlica), pilas de almacenamiento de mineral y otros productos a granel, etc.

Luis Ortuo

EL SOLIDO ELSTICO El semiespacio elstico de Boussinesq. Carga rgida

Supngase ahora que la zapata es infinitamente rgida. Al situar sobre ella la misma sobrecarga g uniforme, , se producir tambin un asiento, pero la rigidez de la zapata no permitir la forma del perfil de asientos que se obtena en el caso de la zapata flexible. En las regiones entre los puntos A y B, y C y D, el asiento ser mayor que el correspondiente a la zapata flexible; mientras que, entre B y C, el asiento de la zapata rgida ser menor.

Figuras y texto procedentes de los borradores del C.T.E y de G&C, II

Consecuentemente, las presiones entre A y B y entre C y D sern, sern en el caso de la zapata rgida, rgida superiores a las correspondientes a la zapata flexible y, por el contrario, entre B y sern menores. Resulta as una distribucin de presiones no uniforme, uniforme caracterizada por unos valores mximos en los extremos y un valor mnimo en el centro.
Luis Ortuo

EL SOLIDO ELSTICO El semiespacio elstico de Boussinesq. Carga rgida

Si el terreno se considerara elstico y de resistencia indefinida la presin bajo los bordes A y D de indefinida, la zapata rgida sera infinita. Dado que la resistencia del terreno es limitada, dichas presiones podrn ser elevadas, pero tendrn un valor finito. En el caso de arcillas (vase figura de arriba), la distribucin de presiones ser en general muy j a la terica del ejemplo j p anterior. Sin semejante embargo, la resistencia limitada del terreno producir en los extremos unas zonas de plastificacin que atenuarn las presiones de borde y las redistribuir hacia el centro de la zapata. En el caso de arenas, dado que la falta de confinamiento en el borde de la zapata, supuesta sta en superficie, no permitira el desarrollo de presiones elevadas, la distribucin tomar en general la forma parablica bli que se i indica di en l la fi figura d de abajo. b j
Figuras y texto procedentes de los borradores del C.T.E.

Luis Ortuo

EL SOLIDO ELSTICO El semiespacio elstico de Boussinesq. Carga en faja distribuida uniformemente

Tomada de G&C, II

Ntese que z no depende de los parmetros elsticos (E, ). Las direcciones principales coinciden con la bisectriz del ngulo 2 y su perpendicular.
Luis Ortuo

EL SOLIDO ELSTICO El semiespacio elstico de Boussinesq. Carga en faja distribuida uniformemente

G&C: Las tensiones principales dependen slo del ngulo 2, luego el lugar geomtrico de los puntos de igual tensin principal es el arco capaz de dicho ngulo respecto a los l b bordes d d de l la f faja. j

Tomada de G&C, II

G&C, G&C II: La integracin de la frmula de Boussineq conduce a un asiento infinito en todos los puntos del semiespacio, lo que se debe a la infinitud de la carga y del espesor del semiespacio. Es precisamente el hecho de que el estrato compresible tenga espesor limitado y que la carga en faja no sea infinita lo que hace que los asientos sean finitos.

Luis Ortuo

EL SOLIDO ELSTICO El semiespacio elstico de Boussinesq. Carga circular

Caso de carga lisa o extensible extensible. Tensiones bajo el centro

Tomada de G&C, II

Ntese que z no depende de los parmetros elsticos (E, ). Carga lisa: Movimiento horizontal libre en superficie: Tensin tangencial nula Carga rugosa: Movimiento horizontal nulo en superficie Luis Ortuo

EL SOLIDO ELSTICO El semiespacio elstico de Boussinesq. Carga circular

Caso de carga lisa o extensible. Tensiones verticales

Tomada de G&C, II

Luis Ortuo

EL SOLIDO ELSTICO El semiespacio elstico de Boussinesq. Carga circular

Caso de carga g lisa o extensible. Isobaras de tensin vertical La figura recoge de forma grfica los incrementos de tensin vertical (z) originados por una presin de intensidad q, distribuida uniformemente en superficie sobre un circulo de radio R. Se representa la mitad del semiespacio, con el eje de (profundidades relativas Z/R) ) ordenadas (p coincidiendo con el centro del crculo. Este grfico permite visualizar de forma sencilla ill e i intuitiva t iti el l concepto t ya apuntado t d en el curso sobre la limitacin del efecto de una carga no infinitamente extensa.

Luis Ortuo

EL SOLIDO ELSTICO El semiespacio elstico de Boussinesq. Carga circular

Caso de carga lisa o extensible. Isobaras de tensin vertical Frecuentemente el lmite de inters se circunscribe a puntos del terreno en los que el incremento de tensin vertical es del orden del 10% de la carga en superficie, para magnitudes habituales de dicha carga (ver CTE, por ejemplo). Ms all de esta zona, la sobretensin recibida por el terreno es, en la mayora de los casos, lo suficientemente pequea para que q sus efectos ( (asientos) ) sean como p comparativamente despreciables (aspecto que se debe comprobar en cualquier caso). La zona definida por el 10% de la carga en superficie forma un bulbo de tensiones que se extiende en profundidad hasta aproximadamente 2 dimetros (2 anchos para zapatas rectangulares) del rea de carga.
Luis Ortuo

De los apuntes sobre el edmetro (Recordatorio)

Cuando la carga es infinitamente extensa, el incremento de presin en superficie se transmite ntegramente en profundidad. Lo mismo se puede suponer (para una gran zona central de la capa compresible al menos) cuando el rea cargada es muy extensa en comparacin con el espesor de dicha capa compresible. Cuando la carga no es muy extensa en comparacin con el espesor de la capa compresible, el incremento de tensin que sufre el suelo necesariamente ir disminuyendo con la distancia distancia. (concepto introducido ya de forma cualitativa al hablar de la consolidacin) Evidentemente, el volumen de suelo que percibe la carga y asienta en consecuencia depender, entre otras cosas, de la extensin del rea cargada. A igualdad del resto de circunstancias, el asiento ser mayor cuanto mayor sea el rea cargada.

Luis Ortuo

De los apuntes sobre el edmetro (Recordatorio)

El bulbo de tensiones y su aplicacin a la profundidad de reconocimiento en la investigacin del terreno

Luis Ortuo

EL SOLIDO ELSTICO El semiespacio elstico de Boussinesq. Carga circular

Caso de carga flexible y lisa o extensible. Asientos.

Para una carga g circular, , uniforme y flexible, , de intensidad p y radio R, , situada en la superficie de un semiespacio elstico de parmetros E y ,se tiene: Asiento A i t bajo b j el l centro: t

Asiento A i t bajo b j el lb borde: d

1 2 s0 = 2 pR R E 2 sR = s0

Asiento medio:

16 pR R 1 2 sm = 0,85 s0 3 E

( Ntese que el asiento elstico es proporcional R)

Luis Ortuo

EL SOLIDO ELSTICO El semiespacio elstico de Boussinesq. Carga circular

Comparacin de cargas lisa y rugosa (tensiones bajo el centro).

Tomada de G&C, II

Para carga rugosa, z s depende de , pero su efecto no parece muy grande.

Luis Ortuo

EL SOLIDO ELSTICO El semiespacio elstico de Boussinesq. Carga rectangular vertical uniforme

El baco proporciona directamente ell incremento de tensin vertical (z) bajo la esquina de un rectngulo cargado con una tensin uniforme (p)

z = Ir p
De nuevo, ntese que z no depende de los parmetros elsticos (E, ).

NOTA: Corregir eje ordenadas en libro

Modificada de G&C, II

Luis Ortuo

EL SOLIDO ELSTICO El semiespacio elstico de Boussinesq. Carga rectangular vertical uniforme

En la mayora de los casos interesar conocer las tensiones producidas por cargas rectangulares no slo bajo sus esquinas, esquinas sino tambin en otros puntos del terreno (bajo el centro, un borde, o incluso bajo un punto alejado del rea cargada). cargada) Para ello se puede hacer uso del principio de superposicin de la teora de la elasticidad, elasticidad de forma que basta combinar de forma apropiada los incrementos de tensin producidos por una serie se e de rectngulos ect gu os debidamente deb da e te seleccionados.

Modificada de G&C, II

Luis Ortuo

EL SOLIDO ELSTICO El semiespacio elstico de Boussinesq. Carga rectangular vertical uniforme

Asiento A i t bajo b j la l esquina i d un rectngulo de t l de d lados l d a y b uniformemente if t cargado d con una carga vertical q (Scheicher, 1926)

sesquina

b(1 2 ) qb =K E
Modificada de G&C, II

Asiento en el centro (lados a/2, b/2mismo K):

scentro

b q (1 2 ) = 4K 2 = 2 sesquina q E
Luis Ortuo

NOTA: Corregir valor ordenada en libro

EL SOLIDO ELSTICO El semiespacio elstico de Boussinesq. Carga rectangular vertical uniforme

Tomada de G&C, II

Para P n , K tiende ti d a infinito i fi it (carga ( en faja). f j ) En la realidad, como ya se indic, el espesor de la capa compresible no es infinito. No suele superar unas decenas de metros, metros como mucho. mucho En esas circunstancias las cargas situadas a una distancia en horizontal de tres o cuatro veces la profundidad mxima de la capa compresible tienen ya una influencia despreciable en el asiento del punto estudiado. estudiado
Luis Ortuo

EL SOLIDO ELSTICO El semiespacio elstico de Boussinesq. Mtodo general para el clculo de tensiones (Newmark)
Un U crculo l cargado d uniformemente if t con una carga q produce a una profundidad z una tensin vertical bajo su centro:

z = q(1 cos3 )
Supongamos p g que p q para la p profundidad z de la figura, todo el crculo R1, cargado uniformemente con una carga unitaria q, da lugar a una tensin z=0,4q, y que todo el crculo R2, cargado uniformemente con una carga unitaria q, da lugar a una tensin z=0,5q . Entonces, la corona circular limitada por los dos crculos, cargada con la misma carga uniforme if q, dar d lugar l a una tensin t i z=0,1q
Luis Ortuo

z
q

= 0,4 tan =

R1 = 0,637 R1 = 0,637 z z

z
q

= 0,5 tan =

R1 = 0,766 R1 = 0,766 z z

EL SOLIDO ELSTICO El semiespacio elstico de Boussinesq. Mtodo general para el clculo de tensiones (Newmark)
Ap partir de las observaciones anteriores, , el mtodo de Newmark consiste en dibujar una serie de coronas con radios calculados de manera que la influencia de cada corona a la profundidad z a la que se quiere calcular el asiento, sea la misma. Para ello basta hallar los radios de crculos sucesivos cuya influencia vaya creciendo en progresin aritmtica. Por ejemplo: z/q R/z
0 0,1 0,2 0,3 04 0,4 0,5 0,6 0,7 0,8 0,9 1 0,000 0,540 0,801 1,036 1 274 1,274 1,533 1,835 2,219 2,774 3,817

Tomada de G&C, II

Luis Ortuo

EL SOLIDO ELSTICO El semiespacio elstico de Boussinesq. Mtodo general para el clculo de tensiones (Newmark)
Con C los l radios di anteriores, t i cada d corona entre t dos d crculos l proporciona una tensin de 0,1q. Si se subdividen ahora las coronas uniformemente por medio de radios, cada trapecio curvilneo transmitir la misma tensin. En el caso dibujado, con 10 sectores (10 radios) cada trapecio transmite 0,01q radios), 0 01q bajo el punto prefijado y a la profundidad prefijada. Dibujando el rea cargada, las coronas y los radios a una escala apropiada, basta por tanto con contar el nmero de trapecios curvilneos ocupados por el rea de carga para determinar fcilmente la tensin vertical buscada.
Tomada de G&C, II

Si la carga no es uniforme, habr que sumar las tensiones correspondientes de cada p trapecio. Cuando se necesita determinar z a varias profundidades, se puede emplear el mtodo de las influencias de Jimnez Salas (ver G&C, II).
Luis Ortuo

EL SOLIDO ELSTICO El semiespacio elstico de Boussinesq. Mtodo general para el clculo de tensiones (Newmark)
Se S dibuja dib j el l rea cargada d a escala l en papel l transparente de forma que la profundidad a la que se desea calcular la tensin coincida con el segmento AB de la figura. figura Se coloca el dibujo a escala sobre el baco, de forma que el centro de las coronas coincida con el punto en el que se desea calcular el asiento. Se cuenta el nmero N de sectores (completos ms fracciones estimadas) cubiertos por el rea cargada. Llamando q a la carga uniforme actuante, el incremento de presin vertical a la profundidad de clculo ser: v=0,001Nq
Tomada de Simons, N. & Menzies, B (2000)

Luis Ortuo

EL SOLIDO ELSTICO El semiespacio elstico de Boussinesq. Mtodo general para el clculo de asientos (Newmark)
Se S basa b en la l misma i id del idea d l clculo l l de d tensiones. t i El asiento i t producido d id en el l centro t de un crculo de radio R cargado con una presin uniforme q es:

1 2 s0 = 2qR E
Por tanto, el asiento producido en el centro de una corona cargada limitada por dos circunferencias de radios R y R, ser

1 2 s = 2qR E
Si la corona se divide en cierto nmero de partes iguales, N, el asiento debido a cada partes ser N veces menor. una de esas p

Luis Ortuo

EL SOLIDO ELSTICO El semiespacio elstico de Boussinesq. Mtodo general para el clculo de asientos (Nomograma 1 de Newmark)
Cada trapecio produce el mismo asiento. Se dibuja en papel transparente el rea cargada, a escala adecuada para que quepa superpuesta en el nomograma. g Se hace coincidir el centro del nomograma con el punto en donde se quiere calcular el asiento. Se cuenta el n de trapecios cubiertos por el rea cargada, estimando las fracciones. El asiento ser:

s = 0,02(1 2 )n q L / E
Tomada de G&C, II

L es la equivalencia en la escala empleada del segmento que figura como mdulo en el nomograma Luis Ortuo

EL SOLIDO ELSTICO Capa elstica homognea sobre base rgida.

Es un caso muy interesante por su analoga con la realidad, ya que frecuentemente existe a una cierta profundidad un estrato comparativamente mucho ms rgido (un sustrato de roca o una capa de suelo muy denso. La influencia de una capa rgida en profundidad es doble. Por un lado se elimina parte d l asiento del i t (que ( se supone despreciable d i bl en la l capa rgida) id ) y por otra, t se modifican difi l las tensiones originadas con respecto a la solucin del semiespacio de Boussinesq.

C Carga en f faja. j Interfaz I t f lisa. li (Egorov, (E 1939) 1939):

En la figura g se representan p las tensiones verticales para carga en faja de intensidad q y ancho 2a. Interfaz lisa (sin tensiones tangenciales en el contacto). Se observa que las presiones a una determinada profundidad crecen al disminuir el espesor de la capa compresible.

Tomada de G&C, II

Luis Ortuo

EL SOLIDO ELSTICO Capa elstica homognea sobre base rgida. Carga en faja. Tensiones verticales (interfaz rugosa. Poulos, 1967b):

Tomadas de G&C, II

Ntese q que las tensiones verticales no dependen p de E, , y apenas p de .

Luis Ortuo

EL SOLIDO ELSTICO Capa elstica homognea sobre base rgida. Carga en faja. Asientos:

Tomada de G&C, II

Ntese la dependencia de E y la gran influencia de , especialmente con capas elsticas de delgado espesor (h pequeo)
Luis Ortuo

EL SOLIDO ELSTICO Capa elstica homognea sobre base rgida. Carga circular (Milovic, 1970)

Tomadas de G&C, II

Primera Pi fi figura: I Boussinesq; I: B i II Adherencia II: Adh i perfecta:; f t III: III Contacto C t t liso li

Luis Ortuo

EL SOLIDO ELSTICO Capa elstica homognea sobre base rgida. Carga rectangular

Tomadas de G&C, II

Luis Ortuo

EL SOLIDO ELSTICO Capa elstica homognea sobre base rgida. Carga rectangular. Mtodo aproximado de Steinbrenner
Se supone p profundidad indefinida de la capa p p compresible p (semiespacio ( p de Boussinesq), con lo que, entre otras cosas, no se tiene en cuenta la modificacin de tensiones originada por la base rgida (se supone que no hay capa rgida) Se calcula el asiento de un punto situado a profundidad z bajo la esquina de un rectngulo cargado. Se calcula el asiento en superficie de la esquina del rectngulo cargado, s0: Se asume que el asiento bajo la esquina del rectngulo para el caso de capa rgida a profundidad z es:

s = s0 sz
Tomada de G&C, II

Luis Ortuo

EL SOLIDO ELSTICO Capa elstica homognea sobre base rgida. Carga rectangular. Mtodo aproximado de Steinbrenner

sz =

pb b [A1 (a, b, z ) B 2 (a, b, z )] 2E A = 1 2 ; B = 1 2 2

Se S calcula l l el l asiento i t en superficie fi i (z=0) y a la profundidad de la capa rgida (z) y se obtiene el acortamiento aproximado de la capa compresible (Otros bacos y tablas en el libro)

Tomada de G&C, II

s = s0 sz
Para varias capas de diversa compresibilidad, se calcula el asiento para distintas profundidades y capas con los parmetros elsticos correspondientes. Por diferencia se halla el asiento de cada una de ellas. Luis Ortuo

EL SOLIDO ELSTICO Ejemplo a resolver por los alumnos

Se carga un rectngulo de 10 x 40 m de lado con una carga uniforme q = 100 kPa. El terreno se puede modelizar en tres capas de E linealmente creciente con valores medios de 20, 20 30 y 40 MPa respectivamente, respectivamente que descansan sobre un sustrato rgido. rgido Suponiendo un coeficiente de Poisson =0,3 para las tres capas, calcular el asiento bajo la esquina del rectngulo.

Luis Ortuo

EL SOLIDO ELSTICO Sistema de varias capas deformables

Es un caso interesante por su analoga con algunas situaciones reales, como cuando existe una capa superior de arcilla endurecida por desecacin, o una capa granular densa de espesor pequeo (uno o dos metros) sobre otro material natural menos rgido. rgido Se aplica adems mucho en ingeniera de carreteras, ya que es habitual estudiarlas mediante sistemas elsticos multicapa.

Si t Sistema bicapa bi
En la figura se representa un caso particular de un rea circular cargada uniformemente sobre un sistema bicapa, siendo el espesor de la capa superior igual al radio del rea cargada. La figura proporciona las tensiones verticales bajo el centro del crculo.

Tomada de G&C, II

Luis Ortuo

EL SOLIDO ELSTICO Sistema de varias capas deformables Sistema bicapa

Se observa que: Cuanto ms rgida es la capa superior, menos tensin llega al estrato compresible subyacente (y en consecuencia asienta menos). La distribucin de Boussinesq es la ms desfavorable.

Tomada de G&C, II

G&C, II: La capa superior, firme, tiene como efecto el repartir las cargas sobre la inferior ms extensamente de lo que resultara en el caso de Boussinesq, a base de un momento flector en dicha capa superior que origina tracciones en la superficie. No es muy seguro que esta capa pueda resistir esas tracciones; si es de acarreos por su falta de cohesin; si es arcilla desecada, porque probablemente estar agrietada. Tan slo en el caso de una capa de verdadera roca sobre una capa de (mdulo de) elasticidad menor (por ejemplo, margas), podra tenerse cierta seguridad de este factor y, y an as, as es ms frecuente que la roca est fisurada o diaclasada y no se pueda confiar en su continuidad. En resumen, no es usual tener en cuenta este efecto en el clculo de las cimentaciones.
Luis Ortuo

EL SOLIDO ELSTICO Cargas rgidas sobre semiespacio elstico homogneo.

Como ya se ha visto, una zapata infinitamente rgida obliga a que los asientos sean iguales en el centro que en los bordes. Boussinesq calcul la distribucin de tensiones bajo una carga circular rgida, y obtuvo el resultado cualitativo de la figura adjunta. Para el caso terico de suelo puramente elstico, la tensin en los bordes resultara infinita. infinita Como el suelo tiene resistencia finita, el entorno de los bordes plastifica, limitndose su tensin y dando lugar a una redistribucin de la misma bajo el rea cargada (ver lnea de puntos en la figura).

Tomada de G&C, II

El inters de conocer esta distribucin de tensiones es en principio doble. Por una parte influir en los asientos y, por otra, ser necesaria para el clculo estructural de la zapata. zapata No obstante, la influencia en los asientos no es muy grande, y stos suelen estimarse p de forma aproximada. Con relacin a los esfuerzos, la distribucin tensional s resulta ciertamente Luis Ortuo relevante, como sugiere la propia figura,

EL SOLIDO ELSTICO Cargas rgidas sobre semiespacio elstico homogneo. Carga en faja. Carga vertical centrada

z =

a 1 ( x / a ) 2
P a

min =
Tomada de G&C, II

C Carga en f faja. j Momento. M t

z =

2 xM a 3

1 1 ( x / a )2
1 1 r2 a2

Placa circular (carga centrada): z = P 2 2a

1 2 P s= 2 aE
Luis Ortuo

EL SOLIDO ELSTICO Cargas rgidas sobre semiespacio elstico homogneo. Placa circular (carga centrada):

Rgida : 1 2 P 1 2 P s= = 0,5 coef = 0,5 E 2 aE a Fl ibl ( P = a 2 q) : Flexible 1 2 2 1 2 P Centro : s = 2qa = coef = 0,637 E E a 4 1 2 P Borde : s = 2 coef = 0,405 E a Coef medio c arg a flexible : 0,405 R lacin Re l i 1,042 ( 4% de d error )
Luis Ortuo

EL SOLIDO ELSTICO Cargas rgidas sobre semiespacio elstico homogneo.

Tomada de la GCOC

Luis Ortuo

EL SOLIDO ELSTICO Cimentaciones de rigidez no infinita.

G&C, II: En el caso de una placa de forma faja infinita, que soporta en su cara superior una carga uniforme (la distribucin de carga influye en la placa por no tener rigidez infinita) la distribucin de presiones de contacto depende de un coeficiente de rigidez infinita), de la placa: Donde:

N=

Ep

12(1 )
2 p

h3

o Ep: mdulo de elasticidad de la placa o p : coeficiente de Poisson de la placa o h : espesor de la placa

G&C, II: Depende tambin de un coeficiente que mide la rigidez del terreno de cimentacin: Donde: Ec o Ec: mdulo de elasticidad del tereno M= 2 o c : coeficiente de Poisson del terreno 1
c

G&C, II: Con ambos coeficientes se puede estimar en la mdulo de rigidez del conjunto, K, que gobierna la interaccin entre la placa y el terreno:
2 2N 1 E p 1 p K= = Ma 3 6 Ec 1 c2

h a

Donde: D d o a: semiancho de la placa

Luis Ortuo

EL SOLIDO ELSTICO Cimentaciones de rigidez no infinita.

2 2N 1 E p 1 p h = K= Ma 3 6 Ec 1 c2 a

Tomada de G&C,II.

Placa flexible en faja infinita sobre semiespacio de Bousinesq. Interfaz lisa. Carga repartida (izquierda) y carga concentrada (derecha). Borowicka, 1939. Tomada de G&C, G&C II. II
Luis Ortuo

EL SOLIDO ELSTICO Cimentaciones de rigidez no infinita.

Algunos comentarios y consideraciones adicionales obre la rigidez relativa cimiento-estructura

Tomadas de G&C, III.

Tomadas del CTE-DB-SE-C

Luis Ortuo

BIBLIOGRAFA
Burland, J.B., Broms, B. & De Mello, V.F.B.(1977): Behaviour of Foundations and Structures. State.of-the-art-report. Session 2. Proc. 9th ICSMFE. Tokyo. Vol. 2, pp 495546. Cdigo Tcnico de la Edificacin (CTE). Documento Bsico SE-C (Seguridad Estructural. Cimientos). Marzo, 2006 FHWA-IF-02-034 (2002): Evaluation of Soil and Rock Properties Gonzlez de Vallejo, L., Ferrer, M., Ortuo, L. & Oteo, C. (2002): Ingeniera Geolgica. Prentice Hall Geolgica Hall. Madrid Madrid. Gordon, J.E. (1987). Structures or why things dont fall down. Pelican Books. Penguin. Gua G de Cimentaciones C en Obras O de Carreteras C (2002). ( ) Direccin General G de Carreteras del Ministerio de Fomento. Series monogrficas. Jimnez Salas, J.A. et al. (1976): Geotecnia y Cimientos II y III. Editorial Rueda. M d id Madrid. Poulos, H.G.. & Davis, E.H. (1974): Elastic Solutions for Soil and Rock Mechanics. John Wiley & Sons, New York. Simons, N. & Menzies, . (2000): A short course in Foundation Engineering. Thomas Telford, London.
Luis Ortuo