Dinamica de Suelos (Alva)
Dinamica de Suelos (Alva)
Dinamica de Suelos (Alva)
DE
ET
LA
BO
R
DINMICA DE SUELOS
2
JORGE E. ALVA HUERTADO, PH. D.
PROFESOR PRINCIPAL UNIVERSIDAD NACIONAL DE INGENIERA
DINMICA DE SUELOS
3
Reservados todos los derechos. El contenido de esta obra esta protegido por la Ley, que
establece pena de prisin adems de las correspondientes indemnizaciones por daos y
perjuicios, para quienes reprodujeran, plagiaren, distribuyeren o comunicaren pblicamente, en
todo o en parte, una obra literaria, artstica o cientfica, o su transformacin, interpretacin o
ejecucin artstica fijada en cualquier tipo de soporte o comunicada a travs de cualquier otro
medio, sin la respectiva autorizacin del autor.
4
PRESENTACINN
Se documenta en estas notas parte del material utilizado en los cursos de Dinmica de
Suelos que se dictan a nivel de integrado y postgrado en la Facultad de Ingeniera civil de la
Universidad Nacional de Ingeniera en Lima, Per
Los cuatro primeros captulos son una adaptacin de las Notas de Clase del profesor R.
V. Whitman, dictadas en el Institutos tcnicos desarrollados por el autor y sus colaboradores a
lo largo de los ltimos aos.
El autor
5
Contenido
CAPITULO I
CAPITULO II
2.1 INTRODUCCIN...12
2.2 VIBRACIONES LIBRES......13
2.3 VIBRACIONES FORZADAS POR LA APLICACIN DE CARGAS
PERIDICAS.17
2.4 VIBRACIONES DEBIDAS A CARGAS TRANSITORIAS..22
2.5 VIBRACIONES FORZADAS PRODUCIDAS POR MOVIMIENTOS PERIDICOS DE
CIMENTACIN..26
2.6 VIBRACIONES DEBIDAS A CARGAS TRANSITORIAS CIMENTACIN.27
CAPITULO III
6
UNIVERSIDAD NACIONAL DE INGENIERA
CAPITULO IV
PROPAGACION DE ONDAS....................70
4.1 INTRODUCCION......................................................................................................................70
4.2 PROPAGACION DE ONDAS EN UN MEDIO INFINITO70
4.3 PROPAGACION DE ONDAS EN UN MEDIO SEMI-INFINITO.73
4.4 PROPAGACION DE ONDAS EN UN MEDIO ESTRATIFICADO.......76
4.5 PROPAGACION DE ONDAS EN BARRAS........78
REFERENCIAS......82
CAPITULO V
CAPITULO VI
22 DE ENERO DE 1582..104
14 DE FEBRERO DE 1619.104
31 DE MARZO DE 1650..105
12 DE MAYO DE 1664........105
10 DE FEBRERO DE 1716.105
28 DE OCTUBRE DE 1746.105
30 DE MARZO DE 1813.........105
20 DE AGOSTO DE 1857..........106
13 DE AGOSTO DE 1868..............106
24 DE JULIO DE 1912........106
24 DE DICIEMBRE DE 1937.................................................................................106
24 DE MAYO DE 1940..........................................................................................106
6 DE AGOSTO DE 1945.......................................................................................106
28 DE MAYO DE 1948..........................................................................................106
21 DE MAYO DE 1950..........................................................................................107
9 DE DICIEMBRE DE 1950...................................................................................107
12 DE DICIEMBRE DE 1953.................................................................................107
15 DE ENERO DE 1958........................................................................................107
17 DE OCTUBRE DE 1966...................................................................................107
19 DE JUNIO DE 1968..........................................................................................107
7
SECCIN DE POST-GRADO - FIC
31 DE MAYO DE 1970....................................................................................108
9 DE DICIEMBRE DE 1970.............................................................................108
20 DE MARZO DE 1972..................................................................................108
3 DE OCTUBRE DE 1974...............................................................................109
CONCLUSIONES.......109
AGRADECIMENTO.......109
REFERENCIAS......109
CAPITULO VII
8
UNIVERSIDAD NACIONAL DE INGENIERA
3 DE OCTUBRE DE 1974......................................................................................119
10 DE NOVIEMBRE DE 1980................................................................................119
7 DE ABRIL DE 1981.............................................................................................119
CONCLUSIONES.............119
AGRADECIMIENTO............120
REFERENCIAS............120
CAPITULO VIII
CAPITULO IX
CAPITULO X
9
2. METODOLOGAS PARA EL ANLISIS.......159
2.1 Compilacin de Whitman (1972) y Arya et al. (1981)......................159
2.2 Compilacin de Gazetas (1991................160
5. CONCLUSIONES Y RECOMENDACIONES.....165
REFERENCIAS...........165
CAPITULO XI
10
SECCIN DE POST-GRADO - FIC
Captulo I
PROBLEMAS DE INGENIERA QUE INVOLUCRAN A LA
DINMICA DE SUELOS
1
UNIVERSIDAD NACIONAL DE INGENIERA
La dinmica de suelos tuvo sus comienzos en Alemania, en 1930, cuando los ingenieros
de DEGEBO (Deutsche Forschungsgesellschaft fur Bodenmechanik) empezaron las primeras
investigaciones cientficas de problemas de cimentacin de maquinaria.
En los Estados Unidos se logr durante los aos 1955 a 1965 un conocimiento
fundamental de la relacin entre comportamiento, teora y propiedades esfuerzo-deformacin
del suelo, principalmente debido al apoyo e inters del gobierno federal en cimentaciones de
antenas grandes de radar.
Cuando recin se instala una turbina bien diseada, se originan fuerzas dinmicas
pequeas. Sin embargo, el desgaste conduce a desbalance en las partes rotativas, por lo que
se desarrollan fuerzas dinmicas.
2
SECCIN DE POST-GRADO - FIC
2. Qu magnitud de la fuerza dinmica debe ser aplicada al suelo durante la vida til
de la maquinaria?
Como ejemplo se tiene un compresor de gas de tipo pistn que est apoyado sobre un
bloque masivo de concreto de cimentacin, que descansa en arena media a densa. Tales
mquinas desarrollan fuerzas dinmicas desbalanceadas asociadas con la aceleracin de
subida y bajada del pistn; las magnitudes de dichas fuerzas pueden ser estimadas con
bastante precisin. En este caso, con el compresor en operacin, los trabajadores de la fbrica
notaron que los movimientos de la cimentacin eran excesivos. Se desarrollaron grietas en las
paredes del local del generador adyacente. El ingeniero de suelos que fue contratado para
recomendar las medidas correctivas debera responder a las siguientes interrogantes:
1. Cmo deberan medirse los movimientos dinmicos para establecer la razn del
movimiento excesivo?
3
UNIVERSIDAD NACIONAL DE INGENIERA
Existe una gran variedad de otros equipos industriales que causan cargas dinmicas en
las cimentaciones, tales como: prensas, vibradores, etc. El problema de disear o mejorar la
cimentacin de dichas maquinarias es similar a aqul para maquinaria reciprocante o rotativa,
con la diferencia que las cargas pueden no ser sinusoidales o peridicas.
Se ilustra una antena muy grande, con una antena tipo plato, apoyada sobre cuatro pies
de acero. Un par de ruedas en la base de cada pie corren en una pista circular de 113 pies de
dimetro, por lo tanto la superestructura rota alrededor de un eje vertical. La pista descansa en
un anillo de concreto. El suelo es una arena gravosa media a densa, compuesta por
fragmentos de coral. Durante el estudio de factibilidad para el sistema de radar, el cual condujo
al establecimiento del criterio final del diseo, se tuvieron que absolver las siguientes
preguntas:
4
SECCIN DE POST-GRADO - FIC
Una estructura muy especializada se utiliz en la plataforma del encendido de las diversas
etapas del cohete Saturno V en las Misiones Apolo. Estas plataformas estaban apoyadas en
pilotes de friccin en arenas y arcillas. Una pregunta importante fue: Causarn las vibraciones
producidas por el lanzamiento una penetracin adicional de los pilotes?. El emplazamiento y
diseo de los centros de control y otras estructuras cerca de las plataformas de lanzamiento
han requerido el estudio de la transmisin de ondas a travs del terreno y los mtodos de
aislamiento de edificaciones contra la vibracin.
Aplicaciones Civiles
Los estudios han indicado que el costo de excavacin utilizando explosivos nucleares
sera mucho menor que el costo de una excavacin convencional.
Construccin de Proteccin
4. Cmo se deberan disear las conexiones entre las estructuras para soportar los
desplazamientos relativos transitorios?
Cimentaciones de Edificios
6
SECCIN DE POST-GRADO - FIC
5. Permanecern estables los diques de los canales del agua de enfriamiento durante
un terremoto?
El efecto de amplificacin del suelo sobre roca fue particularmente importante durante el
sismo de Caracas de 1967. La profundidad de la roca base variaba grandemente en el Valle de
Caracas y a lo largo de la zona costera cercana; el dao severo y la prdida de vidas se
concentr en reas pequeas donde la frecuencia de resonancia de depsitos profundos de
suelo coincidi con la frecuencia de resonancia de edificios altos. Actualmente existe un gran
inters y esfuerzo para incorporar el efecto de las condiciones locales del suelo en los cdigos
de construccin.
Deslizamientos
7
UNIVERSIDAD NACIONAL DE INGENIERA
Presas de Tierra
Se tiene una presa localizada en una regin ssmica inmediatamente aguas arriba de un
rea muy poblada. El suelo en el lugar era una arcilla muy blanda con inclusiones de arena. El
diseo requera respuestas a las siguientes interrogantes en Dinmica de Suelos:
El comportamiento dinmico del suelo debe ser tomado en cuenta en la interpretacin del
hincado de pilotes por medio de martillos de impacto, ya que cada golpe del martillo causa una
onda de esfuerzo en el pilote y de all al suelo circundante. La penetracin durante cada golpe
es rpida y est relacionada a la resistencia dinmica del suelo. Muchas de las explicaciones
recientes sobre el comportamiento del hincado de pilotes toma en consideracin estos
fenmenos dinmicos.
Existen mquinas vibratorias para el hincado de pilotes. Estas operan estableciendo una
condicin resonante para las ondas de salen y llegan al pilote, o una condicin resonante para
el sistema suelo-pilote.
8
SECCIN DE POST-GRADO - FIC
Existen rodillos vibratorios que compactan el relleno hidrulico, en este caso por debajo de
los turbo-compresores. Fue necesario absolver las siguientes preguntas, con referencia a la
recompactacin del suelo:
1. Puede lograrse la densidad final por medio de rodillos vibratorios en la superficie del
relleno? o es necesario excavar y recompactar el suelo en capas?
2. Si el suelo debe excavarse y recompactarse, qu peso del vibrador debe ser
usado?, qu espesor de capa debe usarse? o cuntos pases se necesitan en cada
capa?
En este caso se decidi que el suelo debe excavarse y recompactarse en capas de cuatro
pies de espesor. Se encontr que esta forma de tratamiento era ms econmica que la
compactacin por vibroflotacin o el uso de pilotes.
Los rodillos vibratorios se utilizan para compactar todo tipo de suelo. Las vibraciones se
utilizan en ensayos de laboratorio para determinar las densidades mximas de suelos
granulares.
9
UNIVERSIDAD NACIONAL DE INGENIERA
Por otro lado, la necesidad de considerar la inercia cambia la manera de enfocar los
problemas. Esto es, los problemas en dinmicas de suelos son cualitativamente diferentes de
aquellos en mecnica de suelos, debido a la necesidad de considerar los efectos de las
fuerzas de inercia. Por lo tanto, cualquier curso de Dinmica de Suelos debe enfatizar
grandemente el rol de la inercia.
10
SECCIN DE POST-GRADO - FIC
REFERENCIAS
1) Barkan D.D. (1962), Dynamics of Bases and Foundations, Mc Graw-Hill, New York.
2) Major A. (1962), Vibration Analysis and Design of Foundations for Machines and Turbines,
Akademiai Kiado, Budapest.
3) Ohsaki Y. (196), The Effects of Local Soil Conditions upon Earthquake Damage, Proc.
Specialty Session Soil Dynamics, 7th International Conference on Soil Mechanics and
Foundation Engineering, Mexico.
4) Richart F.E., Woods R.D. y Hall J.R. (1970), Vibrations of Soil and Foundations, Prentice-
Hall, New Jersey.
5) Seed H.B. (1969), The Influence of Local Soil Conditions on Earthquake Damage, Proc.
Specialty Session Soil Dynamics, 7th International Conference on Soil Mechanics and
Foundation Engineering, Mexico.
11
SECCIN DE POST-GRADO - FIC
Captulo II
SISTEMAS LINEALES DE UN GRADO DE LIBERTAD
2.1 INTRODUCCIN
12
UNIVERSIDAD NACIONAL DE INGENIERA
1) Fuerza del resorte: Si el movimiento x es positivo hacia abajo, la fuerza hacia arriba
en la masa es:
kx
M&x&
M&x& + kx = 0 (2.1)
Se puede verificar, por sustitucin directa, que la siguiente expresin es una solucin
de la ecuacin diferencial 2.1
k k
x = A sen t + B cos t (2.2)
M M
t = 0 x = x o, x& = 0
Se encuentra que A debe ser cero y B=xo. Por lo tanto, el movimiento resultante de
esta perturbacin inicial particular es:
k
x = x 0 cos t (2.3)
M
13
SECCIN DE POST-GRADO - FIC
M
T = 2 ( 2.4)
k
1 k
fn = (2.5)
2 M
2 k
= = 2f = (2.6)
T M
t
x = x 0 cos t = x 0 cos 2 f n t = x 0 cos 2 (2.7)
T
1 2
E = k x0 (2.8)
2
pero no hay energa cintica. Por otro lado, cuando el movimiento es cero (x = 0 en t
= T/4, 3T/4, 5T/4, etc.), la velocidad es mxima.
(x& ) max = x 0
14
UNIVERSIDAD NACIONAL DE INGENIERA
NOTA: Este anlisis omite otros dos trminos de energa: el cambio en la energa
potencial de la masa conforme se mueve y el trabajo realizado por la fuerza esttica del resorte
actuando a travs del movimiento dinmico. Sin embargo, estos dos trminos se cancelan
mutuamente en todo momento.
cr = 2 kM (2.11)
D= (2.12)
cr
x = x 0 e -Dt (cos 1 t + D sen 1 t) (2.13)
1
1 = 1 - D2 (2.13a)
15
SECCIN DE POST-GRADO - FIC
2
D -Dt
x = x0 1 + e cos ( 1 t + ) (2.14)
1
donde :
= arc tg D
1
xi
= e Dt = e 2 D
x i +1
xi
= log = 2 D (2.15)
x i +1
x& + kx = 0
k
- t
x = x0 e
(2.16)
16
UNIVERSIDAD NACIONAL DE INGENIERA
Se estudia a continuacin el caso donde existe una fuerza peridica aplicada a la masa.
Esta carga puede describirse por una de las expresiones siguientes:
P = Po sen t
= Po sen 2 f t
Donde
Po = amplitud de la carga
= frecuencia circular de la carga aplicada
f = frecuencia de la carga aplicada
La solucin de esta ecuacin para el caso donde el sistema est inicialmente en reposo
.
(x = x = 0 para t =0) es:
2 2 2
1 - 2 sen t - 2D cos t + e - Dt 2D cos 1 t + 2D + 2 - 1 sen 1 t (2.18)
Po 1
x= 2
k 2 2
2
1 - + 4 D
- Dt
sen t - e sen t
P
x= o (2.19)
2
k
1-
17
SECCIN DE POST-GRADO - FIC
2
1 - 2 sen t - 2D cos t
P
x= o 2
(2.20)
k 2 2
2
1 + 4D
Po sen ( t - )
x= (2.21)
k 2 2
2
2
1 + 4D
donde:
2D
tg = (2.22)
2 - 2
Po
xo = DLF
k
1
DLF = (2.23)
2
2
2
2
1 + 4D
a. Conforme f/fn 0, DLF 1. Fsicamente esto significa que las cargas aplicadas
cambian muy lentamente, comparadas con la velocidad (representada por la
frecuencia natural) a la cual el sistema puede responder; de aqu que la carga sea
esencialmente esttica.
18
UNIVERSIDAD NACIONAL DE INGENIERA
c. puede responder. De aqu que la masa permanece inmvil y resiste la carga aplicada
enteramente por inercia.
1
DLFmax = (2.24)
2D 1 - D2
f. La frecuencia resonante es algo menor que la frecuencia natural y est dada por la
expresin:
r fr
= = 1 - 2D 2 (2.25)
fn
1
DLF = 2
(2.26)
f
1
fn
19
SECCIN DE POST-GRADO - FIC
Dividiendo por la energa almacenada dentro del sistema como energa de deformacin
y cintica, se obtiene la fraccin de la energa que se pierde en cada ciclo:
E 2
= 2 (2.28)
E k
Cuando el sistema est operando a la frecuencia natural, este valor llega a ser:
E
= 4 D = 2 = (2.29)
E
Esto es, la fraccin de la energa perdida por ciclo, es igual al doble del decremento
logartmico. La cantidad E/E a la frecuencia resonante se le conoce como capacidad de
amortiguamiento ().
La energa perdida debe ser proporcionada por la fuerza actuante. El retraso de la fase
entre la fuerza y el movimiento es tal que la integral pdx sobre cada ciclo produce el trabajo
neto sobre el sistema.
Masa excntrica: Los resultados anteriores dan la respuesta dinmica del sistema de
un grado de libertad, independientemente si Po vara o n con la frecuencia, por lo que
representa una solucin completa del sistema. Sin embargo, existe un tipo de carga que es
muy comn, por lo que es til tener un conjunto de grficos para dicho caso especial. Si la
maquinaria de rotacin tiene una masa desbalanceada, la fuerza dinmica ser:
P = M e I 2 sen t (2.30)
donde:
Me = masa excntrica
I = brazo del momento de la masa excntrica
La cantidad Me I2 tiene unidades de fuerza y corresponde a Po.
2
M e I 2 M I f
xo = DLF = e DLF
k M fn
20
UNIVERSIDAD NACIONAL DE INGENIERA
La cantidad (f/fn)2 DLF est graficada en la parte (b) de la Fig. 2.6. Algunas de las
caractersticas principales de estas curvas de respuesta son:
fr 1
= (2.31)
fn 1 - 2D 2
Estas propiedades tambin han sido resumidas en la Tabla 2.1. La Fig. 2.8 se aplica a
este caso especial, como al caso donde Po es independiente de la frecuencia.
21
SECCIN DE POST-GRADO - FIC
Tabla 2.1
PROPIEDADES DE LA RELACIN FACTOR DE CARGA DINMICA VS. FRECUENCIA
Sistema Masa
Fuerza Actuante
Excntrica
Respuesta adimensional
a f=0 1 0
Respuesta adimensional
a f 0 1
Relacin de frecuencia
1
resonante fr/fn 2
1 2D 1 2D 2
Respuesta adimensional SECCIN DE POST-GRADO - FIC
1 1
a f = fr 2
2D 1 D 2D 1 D 2
Respuesta adimensional 2
1 f
2 fn
para f << fr
1 f 2
fn 1 f
f >> fr fn
0 t<0
M&x& + x& + kx = (2.33)
Po t0
22
UNIVERSIDAD NACIONAL DE INGENIERA
x=
Po
k
[
1 - e -Dt cos t ] (2.34)
Carga rampa: Como se ilustra en la Fig. 2.10, una carga rampa involucra un incremento
lineal de carga, seguido de carga constante. La ecuacin gobernante para amortiguamiento
cero es:
0 t<0
t
M&x& + kx = Po 0t (2.35)
Po t
La solucin del movimiento debe obtenerse en dos pasos. Primero se obtiene una
solucin para 0 t , con las condiciones iniciales x = x& = 0 para t = 0:
Po t T
0t x= - sen t (2.36a)
k 2
P
x& = o (1 - cos ) = (A cos - B sen )
k
Po
t x= 1 +
T
[sen (t - ) - sen t ] (2.36c)
k 2
La respuesta para una carga rampa consiste de vibraciones libres superpuestas sobre la
solucin esttica, como se ilustra en la Fig. 2.10. La amplitud mxima de movimiento es funcin
de /T, como se muestra en la Fig. 2.11. Cuando el perodo natural del sistema es pequeo,
23
SECCIN DE POST-GRADO - FIC
0 t<0
M&x& + kx = Po 0t (2.37)
0 t>
Po
0 t x= (1 - cos t ) (2.38a)
k
2P
t x = o sen sen ( t + ) (2.38b)
k T
donde:
cos - 1
tg = (2.38c)
sen
Es importante mencionar que los movimientos causados por un pulso de corta duracin
(/T << 1) pueden obtenerse directamente con la aplicacin de la ecuacin impulso-momento.
El impulso es Po, ya que en este caso no se desarrolla fuerza en el resorte durante la carga. El
momento mximo durante la vibracin residual subsiguiente es M x0, de donde:
Po Po P 2
xo = = = o
M k k T
Este es el mismo resultado que el dado por la Ec. 2.38b para valores de /T pequeos.
24
UNIVERSIDAD NACIONAL DE INGENIERA
0 t<0
M&x& + x& + kx = Po sen t 0t (2.39)
0 t
Para 0 t , esta ecuacin es igual a la Ec. 2.17 y la solucin es la ecuacin 2.18. Para
amortiguamiento cero:
sen t - sen t
Po
0 t x= 2
(2.40a)
k
1-
La vibracin residual se determina igualando la vibracin libre y la forzada en t=. Para
amortiguamiento cero:
- sen
P
t = x= o 2
= A sen + B cos
k
1-
Po 1 - cos
x& = 2
= (A cos B sen )
k
1-
Po 2 (/) n
t x= sen sen ( t + ) (2.40b)
k 1 - (/ ) 2
2
donde:
- sen
tg = (2.40c)
1 - cos
Para el caso especial donde /=1, el numerador y denominador de las ecuaciones
(2.40) se anulan. Una expresin alternativa para este caso se obtienen aplicando la regla de
Hospital:
Po
0 t x= (sen t - t cos t ) (2.41)
2k
25
SECCIN DE POST-GRADO - FIC
- Po
t x= n cos t (2.42)
k
Las ecuaciones 2.41 y 2.42 se grafican en la Fig. 2.14a. La amplitud del movimiento
aumenta linealmente con el nmero de ciclos de fuerza aplicada, pero permanece finita en
tanto que la duracin total de la fuerza es finita. De este modo, una fuerza sinusoidal con =
causa movimientos muy grandes si se aplican muchos ciclos de la fuerza.
Las ecuaciones 2.41 y 2.42 se aplican para amortiguamiento nulo. Con amortiguamiento,
los movimientos se producen hacia los lmites indicados en la Fig. 2.14b. El nmero total de
ciclos de carga requeridos para lograr la respuesta total amortiguada es aproximadamente
1/2D. As, con 5% de amortiguamiento, el movimiento aumenta durante los primeros 3 ciclos
de carga y luego los ciclos sucesivos casi no aumentan el movimiento. Si 2n<1/2D, la
respuesta mxima es menor que para vibraciones forzadas, pero si 2n>1/2D, la mxima
respuesta transitoria es esencialmente la misma que para la respuesta de estado constante.
Esta seccin considera el caso cuando no se aplica fuerza a la masa, sino el apoyo del
resorte experimenta movimiento en la cimentacin (ver Fig. 2.16). En las siguientes
expresiones, x es el desplazamiento absoluto de la masa, s es el desplazamiento absoluto
del terreno e y es el desplazamiento relativo entre la masa y el terreno. Luego, y+s=x. Las
fuerzas actuantes en la masa son:
26
UNIVERSIDAD NACIONAL DE INGENIERA
S = So sen t
Por lo tanto, las soluciones de la seccin 2.3 son tambin aplicables, pero reemplazando
Po por MSo 2 y adems x por y. En particular, considerando solamente vibraciones forzadas
y omitiendo vibraciones libres, la amplitud yo de los desplazamientos relativos est dada por:
2
y o = S o DLF (2.44)
donde DLF est definida por la ecuacin 2.23. Por lo tanto, la Fig. 2.6b puede ser usada
como un grfico de yo/So versus f/fn. Se desarrolla resonancia a f/fn1, tanto para movimientos
peridicos de cimentacin como para fuerza peridica aplicada.
Movimiento seno-verso: Un pulso simple de movimiento del terreno seno verso (Fig.
2.17) est dado por:
0 t<0
S
2
S = o (1 - cos t) 0 t = (2.45)
2
0 t >
27
SECCIN DE POST-GRADO - FIC
2
3
t y = - S o 2 sen sen t + (2.46b)
1-
M &x& = k y max
&x& = 2 y max
28
UNIVERSIDAD NACIONAL DE INGENIERA
Sv = y max (2.48)
1
E max = MSv 2 (2.49)
2
Debido a las relaciones simples entre y max, Sv y &x& max , es conveniente graficar la
respuesta de un sistema dinmico en el papel grfico espectral mostrado en la Fig. 2.19. Con
este papel especial, es posible presentar simultneamente el espectro de respuesta de
desplazamiento relativo, seudo-velocidad y aceleracin absoluta. De este modo, si un sistema
con una frecuencia natural de 10 cps experimenta un desplazamiento mximo relativo de 1
pulgada, la aceleracin absoluta mxima (10 g) y la seudo velocidad (63 pulg/seg) pueden
leerse directamente del grfico.
El espectro de respuesta para un pulso simple de movimiento seno-verso del terreno ha
sido regraficado en la Fig. 2.20 para So=1 pulg. y =62.8 rad/seg (f=10 cps). Puede verse que
un sistema simple de un grado de libertad con una frecuencia natural de 5 cps, experimentara
un desplazamiento mximo relativo de 1.3 pulg (y &x& max = 3.2 g), mientras que con fn = 25 cps:
y max = 0.2 pulg y &x& max = 14 g. Varias relaciones muy importantes entre la respuesta y el
movimiento del terreno pueden entenderse con la ayuda de esta figura:
b) Para sistemas con / grande, el sistema muy rgido sigue al movimiento del terreno
muy de cerca, de modo que los desplazamientos relativos son pequeos. La
aceleracin absoluta mxima de la masa es similar a la mxima aceleracin del
terreno, pero puede ser algo amplificada debido a la dinmica del sistema. En el
ejemplo de la Fig. 2.20, la aceleracin mxima del terreno es 5.1 g, mientras que la
mxima aceleracin absoluta de estructuras con / grande vara de 10.2g a algo
ligeramente menor que 10 g.
Las relaciones entre movimiento del terreno y respuesta para / < <1 y para / > > 1
son vlidas para cualquier movimiento del terreno. Los ejemplos siguientes consideran los
factores que afectan la respuesta para / 1.
29
SECCIN DE POST-GRADO - FIC
2
y max = So 2 (2.50)
1-
Para / = 1/3, 3, 5,7,..., la respuesta mxima es realmente algo menor que la dada por
la ecuacin 2.50, ya que para estas condiciones los trminos de coseno en los parntesis de la
ecuacin 2.46a nunca son simultneamente la unidad. Para todos los otros valores de /,
incluyendo aquellos en la cercana de 1/3, 3, 5, 7,...., existir algn tiempo durante un tren
infinito de pulsos cuando los trminos de coseno sean la unidad simultneamente.
0 t<0
S 1 2
S o (1 - cos t + n (1 - cos n t) 0 t = (2.51)
2
0 t >
30
UNIVERSIDAD NACIONAL DE INGENIERA
31
SECCIN DE POST-GRADO - FIC
M M
x (b) x
(a)
k
k
M M
Resorte de extensin
Resorte
de corte
(a) Viga real
M Resortes de corte
32
UNIVERSIDAD NACIONAL DE INGENIERA
x
-Xoe-wDt
Xo
-Xo
-Xoe-wDt
Punto de tangencia casi en el punto de movimiento mximo
t
Movimiento libre
Movimiento forzado
Movimiento total
33
SECCIN DE POST-GRADO - FIC
5 5
D=0
D=0
D = 0.1
4 4
D = 0.1
3 3
Xo Xo M
Po / k Mel 0.2
D = 0.2
2 2
0.3 0.3
0.4 0.4
0.5 0.5
0.6 0.6
1 1
0 0
0 0.5 1.0 1.5 2.0 0 0.5 1.0 1.5 2.0
f / fn f / fn
34
UNIVERSIDAD NACIONAL DE INGENIERA
0.8
0.7
0.6
Relacin de amortiguamiento critico D
0.5
0.4
0.3
0.2
0.1
0 1 2 4 6 8 10 12
14 15
Xo k Xo M
Mel en
Po
resonancia
35
SECCIN DE POST-GRADO - FIC
D
0.1
0.2
180 0.3
Angulo de Fase,
0.4
0.5
0.6
1.0
90
0 1.0 2.0
0
x
P
Po
2
k
Po
Po
k
t t
P 2 T
10T =
Po = 4
Xk
3
Po 1
t 0
2 3 4
tiempo
36
UNIVERSIDAD NACIONAL DE INGENIERA
X max k
Po 1
0 1 2 3 4
/T
2
=
5
4
T
P
Po Xk
1
Po
=
T
10
t 0
0 2 3 4 5
X max k
Po 1
0
0 0.2 0.4 0.6 0.8
/T
37
SECCIN DE POST-GRADO - FIC
t
X max =
2
Xk Xk
Po Po
1
2D
t t
1
2D
t = 2
Po
= t
2.0
X max k
Po
1.5
Xmax
1.0
Xmax durante vibraciones residuales
0.5
0
0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0
/T
Figura 2.15 Mxima respuesta para pulso seno medio
38
UNIVERSIDAD NACIONAL DE INGENIERA
x
y
s
Figura 2.16 Sistema masa - resorte-amortiguador con movimiento de apoyo
2 2
= 1 1
T 4 2
1 1
y+s
So
0 0
1 1 3 1 3
2 t/ 2 2 2
-1
2 2
1 3
2
1 1
0 0
1 3 1 3
2 2 2 2
-1 -1
-2
So 1
0
0 t 1
Figura 2.17 Respuestas tpicas a un movimiento de cimentacin seno verso (un ciclo)
39
SECCIN DE POST-GRADO - FIC
YR
max 1
S
0
0 1 2 3 4 /T
Y + S max
1
So
0
0 1 2 3 4 /T
2
Y max
1
So
0
0 1 2 3 4 / T=
40
SECCIN DE POST-GRADO - FIC
500
50
0
50
00
0
25
g
10
25
00
00
250
g
0 s
10
g
50
es
n
0
pl
ci
g
az
50 a
er
am
25
el
0
ie
Ac
g
nt
25
o
100
,p
ul
10
g
0
g
10
50
g
0
Pseudo Velocidad, pulg/seg.
50 5.
25
g
5
2.
1.0
0
25
g
0
1.
5.
0
g
5
0.
2.
5
5
g
2
10 0.
1.
1
0.
0
g
50
0.
0
50
0.
g
5.0
0.
0.
0
2
25
25
5
50
0
g
02
0.
g
0.
0.
01
0.
0
2.5 10
10
0
01
0
g
0.
0
g
0.
0.
0.
00
50
0
5
5
00
00
g
0.
0
0.
1.0
0.10 0.25 0.50 1.0 2.5 5.0 10 25 50 100 250 500 1000
Frecuencia, cps.
41
Figura 2.19 Lamina espectral para graficar el espectro de respuesta
UNIVERSIDAD NACIONAL DE INGENIERA
500
n
n=4
De
10
00
250
sp
0
g
10
l az
am
i n
ie
ac
nt
er
o,
el
pu
Ac
lg
s.
100
10
0
10 g
Pseudo Velocidad, pulg/seg.
50
g
5
50
25
0
1.
g
2.
10
g
25
5
0.5
g
25
2.
n=1 0.
5
g
10
10
1
0.
g
5.0
0.
1
0.0
1
2.5
g
Seno verso
So = 1 pulgada
F = 10 cps
0.
01
1.0
g
0.10 0.25 0.50 1.0 2.5 5.0 10 25 50 100 250 500 1000
S 500
So
10
n,
250
00
i
ac
er
50
el
Ac
0
t
2
25
0
100
10
50
0
50
25
Respuesta al movimiento
combinado del terreno
10
25
5 ciclos a
50 cps.
5
D
5
2.
es
Pseudo Velocidad, pulg/seg.
2.
pl
5
10
az
am
ie
1 ciclo a 10 cps.
nt
0 01
1
1.
o,
0.
pu
lg
5 .
0.
.
0.
5
5.0
25
0.
0.
25
2.5 1
0. 00
1
0.
0.1
1.0
0.1 0.25 0.50 1.0 2.5 5.0 10 25 50 100 250 500 1000
&& max
&& max = S
X
Rango de frecuencias contenidos en el
Y max = S max terremoto.
Pseudo-Velocidad Sv (escala logartmica)
No amortiguado
Altamente amortiguado
10 10
8 8
Sv-pie/seg
Sv-pie/seg
6 6
4 4
2 2
0 0
0 0.5 1.0 1.5 2.0 2.5 3.0 0 0.5 1.0 1.5 2.0 2.5 3.0
PERIODO-SEGUNDOS PERIODO-SEGUNDOS
44
SECCIN DE POST-GRADO - FIC
Captulo III
SISTEMAS LINEALES DE VARIOS GGRADOS DE LBERTAD
3.1 INTRODUCCIN
Cuando un sistema dinmico tiene dos grados de libertad (2-GDL), la respuesta dinmica
puede evaluarse por solucin directa de las ecuaciones diferenciales; es decir, las mismas
tcnicas que se utilizan para sistemas de 1-GDL. Cuando un sistema tiene ms de dos grados de
libertad, es muy tedioso el obtener soluciones directas a las ecuaciones diferenciales, por lo que
se utiliza generalmente en ese caso el mtodo de los modos. En este captulo se ilustrarn los
dos mtodos para enfatizar algunos conceptos fundamentales y derivar frmulas simples
importantes. Mayores detalles se obtienen de las referencias presentadas.
1) Sistemas con dos masas, donde cada masa tiene solo un grado de libertad.
2) Sistemas con una masa, donde esta masa tiene dos grados de libertad para su
movimiento.
Estos diferentes tipos de sistemas se ilustran en la Fig. 3.1. El tipo 2a puede tratarse
simplemente como dos sistemas de un grado de libertad cada uno, utilizando los mtodos del
captulo anterior. En esta seccin y en la prxima se tratarn los tipos 1 y 2b.
45
UNIVERSIDAD NACIONAL DE INGENIERA
Las ecuaciones que describen la respuesta de este sistema pueden obtenerse del anlisis
que sigue. Las fuerzas actuantes en la masa inferior son:
M1 &x& 1 + (k 1 + k 2 ) x1 - k 2 x 2 = 0
(3.1)
M 2 &x& 2 + k 2 x 2 - k 2 x1 = 0
- M 1 A 1 2 + ( k 1 + k 2 ) A1 - k 2 A 2 = 0
(3.3)
2
- M 2 A 2 + k 2 A 2 - k 2 A1 = 0
2
A 2 = k 1 + k 2 - M1 = k2 (3.4)
A1 k2 k 2 - M 2 2
+
4 - k 1 k 2 + k 2 2 + k 1 k 2 = 0 (3.5)
M1 M2 M1 M 2
La solucin de la ecuacin 3.5 dar dos races reales y positivas de . La razn de amplitud
correspondiente a cada puede determinarse de la ecuacin 3.4. La magnitud real de cada
amplitud depender de la naturaleza de las condiciones de frontera que excitan las vibraciones
libres.
46
SECCIN DE POST-GRADO - FIC
Una solucin de la forma de las ecuaciones 3.2 es posible. Existen dos modos de vibracin posibles,
correspondientes a las dos races I y II de la ecuacin 3.5; el movimiento total puede obtenerse
por superposicin.
Ms an, una solucin de la forma de las ecuaciones 3.2 pero con B cos t en lugar de A
sen t, tambin dar las ecuaciones 3.3 a 3.4, y la solucin total de la vibracin libre sera de la
forma:
A
x 2 = 2 I (A 1I sen I t + B1I cos I t + ...
A1I
y + L
47
UNIVERSIDAD NACIONAL DE INGENIERA
Fuerzas Horizontales
Momentos Alrededor de Z
I = M L 2 + Io (3.7)
Los lados de la izquierda de las ecuaciones 3.6 son exactamente las ecuaciones de
vibracin libre en los modos horizontal y de cabeceo. Estas ecuaciones son acopladas, ya que
aparece en la ecuacin de y, y viceversa. Si el resorte restaurador y horizontal el amortiguador
estuviesen alineados con el CG, de modo que L = 0, entonces estas ecuaciones se desacoplaran.
M&y& + k y y = - ML &&
(3.8)
I&& + k = - ML &y&
* o
Un anlisis cuidadoso probar que y y deben estar exactamente en fase ( 180 fuera de fase) si no existe amortiguamiento.
48
SECCIN DE POST-GRADO - FIC
- 2 M B1 + k y B1 = ML B 2 2
(3.10)
2 2
- I B 2 + k B 2 = ML B1
Estas ecuaciones pueden rearreglarse del siguiente modo, dejando dos incgnitas: y la
relacin de amplitudes = B2 / B1
ML 2 - M 2 + k y
= B2 = = (3.11)
B1 - I 2 + k ML 2
(ML 2 )2 - (k - I 2) (k y - M 2) = 0 (3.12)
Io 4 - k + k y 2 + k k y = 0
(3.13)
I I M MI
2
y 2 2
2
2 I
y Io y
2
= 1+ 1 + - 4 (3.14)
2 Io I
donde:
k
= ; la frecuencia circular natural de solo cabeceo.
I
ky
y = ; la frecuencia circular natural de solo traslacin horizontal.
M
Puede usarse la Fig. 3.6 para determinar las dos races. La inspeccin de esta figura
mostrar:
I II
cuando y : y , y y
I II
cuando y : y y , y
49
UNIVERSIDAD NACIONAL DE INGENIERA
Por lo tanto, la frecuencia natural acoplada ms pequea es menor que las dos frecuencias
desacopladas, mientras que la frecuencia natural acoplada ms alta es mayor que las dos
frecuencias desacopladas.
Una vez que los valores de la frecuencia natural circular han sido hallados, los valores
correspondientes de la relacin de amplitudes = B2 / B1 , pueden encontrarse de la ecuacin 3.11
La relacin de amplitud ser diferente en cada una de las dos frecuencias naturales. Los
valores reales de B2 y B1 dependern de las condiciones iniciales que causen las vibraciones libres.
Soluciones de la forma:
y = A1 sen t , = A 2 sen t
tambin conducirn a la ecuacin 3.14 para las frecuencias naturales y a la ecuacin 3.11
para la razn de amplitud. La aplicacin de estas ecuaciones y grficos al problema de vibracin
libre acoplada se ilustra en el Ejemplo 3.2.
Esta seccin considera las vibraciones forzadas de un sistema de 2-GDL. Los principios se
ilustrarn por medio de un problema de vibracin acoplada en un sistema de una masa. Los mismos
principios se aplican a vibraciones desacopladas de un sistema de dos masas.
Asumamos que estas ecuaciones tienen soluciones de la forma (Nota: se asume que las
vibraciones libres son amortiguadas. Puede demostrarse que el ngulo de fase entre la fuerza
aplicada y la respuesta es cero 180, si es que se omite el amortiguamiento):
50
SECCIN DE POST-GRADO - FIC
(- M 2 + k y) y o - ML 2 o = P o
(3.17)
- ML y o + ( - I + k ) o = T o
2 2
(k - I 2) Po + ML 2 To
yo = 2 2 2 2
(k - I )(ky - M ) - (ML )
(3.18)
2
ML Po + (k y - M ) TO2
o =
(k - I 2) (ky - M 2) - (ML 2 )2
[ I 2 2 II 2
][
MIo (y ) - (y ) -
2
]
Por lo tanto, las ecuaciones 3.18 se pueden reescribir como:
I L
(2 - 2) P o + 2 T o
(MIo) Io
yo =
[( y ) - 2] [( yII ) 2 - 2]
I 2
(3.19)
L 2 2 2 To
P o + ( y - )
I Io
o = o I 2 2 II 2
[( y ) - ] [( y ) - 2]
51
UNIVERSIDAD NACIONAL DE INGENIERA
To LI
2
2
yo = 2 2
k Io [( y I ) - 2] [( y II) - 2]
(3.20)
2
T o I
2
(y - )
2
o = 2 2
k Io [( y I ) - 2] [( y II ) - 2]
Las cantidades en el lado derecho sirven como factores de carga dinmicos. Supongamos
que:
y I < < y < y II
Angulo de fase a la
Rango de Frecuencia fuerza motora
y
< y
I 0o 0o
y < < y
I 180o 180
y < < y
II 180o 0o
y <
II 0o 180
Las formas generales de los factores de carga dinmica se muestran en la Fig. 3.10
El punto alrededor del cual la cimentacin rota (el punto donde la traslacin horizontal es
cero) se encuentra a una altura h por encima del punto Z, tal que:
yo + h o = 0
2
h = L
2 2
- y
52
SECCIN DE POST-GRADO - FIC
> y, el centro de rotacin est por encima del centro de gravedad, aproximando gradualmente al
centro de gravedad conforme llega a ser muy grande.
El uso de los procedimientos presentados en esta seccin para calcular la amplitud de
movimiento de un sistema dinmico se ilustra en el Ejemplo 3.3.
Vibracin Forzada Acoplada de Sistema de una Masa Amortiguada: Con los trminos de
amortiguamiento, las ecuaciones de equilibrio dinmico son:
Po 0 , To = 0
y superponer los resultados
y = ( yo )T sen ( t - 1 - 2)
(3.22)
= ( o ) T sen ( t - 1)
Sustituyendo las ecuaciones 3.22 en las ecuaciones 3.21, expandiendo las funciones seno y
agrupando trminos en sen t y cos t, conduce a cuatro ecuaciones simultneas que pueden ser
resueltas para yo, o, 1 y 2. Esta solucin es:
2
(y o)T = L T o (3.23a)
Io 2
X +Z
2
2
T o
( y - 2 ) 2 + (2 D y y ) 2
( o)T = (3.23b)
Io X +Z
2 2
[ ][
X = ( y I ) 2 2 ( y II) 2 2 ] (3.23c)
2
I
Z =
2
[( 2 - 2 ) 2 (2 D y y ) 2 + (y 2 - 2 )2 (2 D ) 2
Io (3.23d)
ML2
+ (2 D y y ) 2 (2 D ) 2 + 2 (2 D y y ) (2 D ) 4 ]
I
2 D y y
tg 2 = 2
(3.23e)
y - 2
53
UNIVERSIDAD NACIONAL DE INGENIERA
[ ]
2
(2 D ) (y - 2 )2 + (2 D y y 2) + ML 4 (2 D y y )
2
tg 1 = I (3.23f)
[ ]
2
2 2 2 2 ML 4 ( 2 - 2 )
( - 2) (y - 2 ) + (2 D y y ) - y
I
y = (y o ) p sen ( t - 3)
(3.24)
= ( o) p sen ( t - 3 - 4)
donde:
I ( 2 - 2 )2 + (2 D )2
( yo )p = Po (3.25a)
M Io 2
X +Z
2
2
L
( o ) p = P o 2 2
(3.25b)
Io X +Z
2 D
tg 4 = 2 2
(3.25c)
-
(2 D ) [( - ) + (2 D ) ] + ML (2 D )
2
2 2 2 2 4
y y
tg 3 = I (3.25d)
( - ) [( - ) + (2 D ) ] - ML ( - )
2
2 2 2 2 2 2 4 2 2
y
I
donde X y Z estn dadas por las ecuaciones 3.23 c y 3.23 d. En el caso tpico, donde ambos
Po y To son diferentes de cero (cuando una fuerza horizontal acta en algn punto por encima de la
base de la cimentacin), las ecuaciones 3.22 y 3.24 pueden aadirse. Con el objeto de obtener la
amplitud del movimiento superimpuesto, las diferencias de fase deben tomarse en cuenta. El uso de
estas ecuaciones se ilustra en los Ejemplos 3.4 y 3.5. El trabajo numrico asociado es bastante
tedioso, por lo que se necesita un computador.
Ejemplo 3.1
k1 = k 2 = k M1 = M 2 = M
54
SECCIN DE POST-GRADO - FIC
k 3 5
2
=
M 2
k k
I = 0.618 II = 1.62
M M
A2 = k 2
=
A1 3 5 5 -1
k 1 -
2
A2 A2
= 1 = 1.62 = II = - 0.618
A1 I A1 II
Los dos patrones de movimiento se muestran en la Fig 3.3. Estos dos modos de vibracin se
denominan modos normales de vibracin.
Las relaciones de amplitud 1 y 2 son exactamente las mismas si se asume los movimientos
B cos t A sen t, es decir:
B2 = B2 =
I II
B1 I B1 II
La solucin completa para la vibracin libre del sistema de dos masas es:
k k k k
x1 = A1 I sen 0.618 t + B1 I cos 0.618 t + A1 II sen 1.62 t + B1 II cos 1.62 t
M M M M
k k k k
x 2 = 1.62 A1 I sen 0.618 t + B1 I cos 0.618 t - 0.618 A1 II sen 1.62 t + B1 II cos 1.62 t
M M
M M
A1 I = A1 II = 0
B1 I = 0.724 x o
B1 II = 0.276 x o
55
UNIVERSIDAD NACIONAL DE INGENIERA
k k
x1 = 0.724 x o cos 0.618 t + 0.276 x o cos 1.62 t
M M
k k
x 2 = 1.171 x o cos 0.618 t - 0.171 x o cos 1.62 t
M M
Ejemplo 3.2
Dados : Un sistema tal como se muestra en la Fig. 3.5, con los parmetros siguientes:
y = yo = 0.1 pie = o = o en t = 0
9.64 x 106
y = = 39.5 x 102 = 62.8 rad/seg
2440
(f n ) y = 10 cps
6.09 x 108
= = 19.3 x 102 = 43.9 rad/seg
3.15 x 105
(f n ) = 7 cps
2
y Io = 0.432
= 2.04
I
I I
y = 37.9 rad/seg (fn ) y = 6.04 cps
II II
y = 109.0 rad/seg (fn ) y = 17.5 cps
56
SECCIN DE POST-GRADO - FIC
Relaciones de Amplitud
ML 1 2440 x 8.56 1
I = x 2
= x
I 315,000 (1.158 )2 - 1
-1
y
1
= 0.0663 = 0.195 pies-1
0.34
1
II = 0.0663 x 2
= - 0.0801 pies-1
(0.403 ) - 1
Las vibraciones libres en cualquier modo involucran cabeceo y traslacin. Desde que I > II,
el cabeceo ser relativamente ms importante (en relacin a la traslacin horizontal) en el modo
primario que en el secundario.
La Fig. 3.8 presenta la naturaleza de los movimientos en dos puntos diferentes de la masa
vibrante.
57
UNIVERSIDAD NACIONAL DE INGENIERA
Ejemplo 3.3
Considere el sistema del Ejemplo 3.2 y asuma que este sistema est sujeto a la carga
mostrada en la Fig. 3.11. Esta fuerza aplicada fuera del CG puede ser resuelta en una fuerza
vertical a travs del CG ms un momento To = 15,830 pie-libra alrededor del eje horizontal a travs
de Z. La fuerza ser positiva hacia abajo; y el momento correspondiente tambin ser positivo.
La constante de resorte para el movimiento vertical se tomar como kx = 14.35 x 106 lb/pie.
La frecuencia natural para el movimiento vertical, que ser desacoplado si los resortes verticales
son simtricos con respecto al CG, es igual a:
1 14.35 x 106
(f n ) x = = 12.2 cps
2 2440
Movimientos verticales
Po 1
xo = 2
kx
1 -
x
donde x es la frecuencia circular natural no amortiguada para movimientos verticales: x =
2 (12.2) = 36.6 rad/seg. A la frecuencia normal de operacin = 2 (30) = 189 rad/seg:
El signo negativo significa que el movimiento est fuera de fase con la fuerza;
Las ecuaciones 3.20 se aplican (existe fuerza horizontal cero). En la frecuencia de operacin
normal:
yo = 0.00004405 pies
o = -0.00000458 radianes
Movimientos Totales
58
SECCIN DE POST-GRADO - FIC
Conceptos Bsicos
n=r
u i (t) = D n (t) n n (i) (1)
n =1
Dn(t) es el factor de carga dinmico para el modo n-simo. Dn(t) depende solamente de n,
Dn y el tiempo-historia de la carga aplicada, y es idntico a la respuesta de un sistema de 1-GDL que
tiene estos y D. As, cualquier solucin (para carga peridica o transitoria) ya obtenida para 1-
GDL puede ser utilizada. Este hecho proporciona al anlisis modal su mayor ventaja: la tarea de
evaluar el comportamiento complicado de un sistema complejo de n-GDL, se reduce a determinar y
superponer las respuestas de n sistemas de 1-GDL.
n=r
(u i )max (Dn )max n n (i )
n =1
o mejor an:
n =r
(u i )max [(D
n =1
)
n max n n (i )]
2
Usualmente solo unos cuantos modos deben retenerse, y a menudo solo el primer modo se
utiliza, especialmente para trabajo aproximado.
59
UNIVERSIDAD NACIONAL DE INGENIERA
Para este caso, la ecuacin (1) proporciona movimiento relativo y n es determinado por
..
cargas Mi S aplicadas a cada masa. El espectro de respuesta de Ymax puede ser usado para
determinar (Dn)max para cada modo.
Para trabajo aproximado (tal como diseo inicial) se asume algunas veces que existe
solamente un modo para un edificio de muchos pisos y que la forma caracterstica para este modo
es una lnea recta con desplazamiento mximo en el techo (esta ltima suposicin ahorra el tener
que calcular la forma caracterstica). Estas suposiciones corresponden a tener un coeficiente
ssmico que vara linealmente con la altura hasta un valor mximo de Sa en la parte superior del
edificio. As, el corte pico en la base Vmax del edificio es:
i= N
Vmax = Sa ( Zi / H ) Mr
i= i
donde Zi es la elevacin de la masa i-sima y H es la altura total del edificio. Si los picos son
igualmente espaciados y de igual masa
1
V max = Sa N M
2
1 F r
r nr
n =
n
2
M
r
r nr
2
M r
r nr
n =
M r
r nr
2
60
SECCIN DE POST-GRADO - FIC
Forma de Modo:
2
-
12 = 1 11 =
k2
2
= k 2 M 2 1
k1 + k 2 - M1 1 k2
2
22 = 1 21 = igual pero con 2
Factor de Participacin: F1 = 0 F2 = Po
Po
n = 2
(M1 n1 + M 2) n 2
2
n sen t
Dn (t) = 2 2
n -
2
2
1
2
2 1 1
+ - +
M +M M +M
2 2 M +M M +M
2 2
1 11 2 1 21 2 1 11 2 1 21 2
x 2 = Po sen t
(1 - ) (2 - )
2 2 2 2
1
(k1 + k 2 - 2 M1)
M1 M 2
61
UNIVERSIDAD NACIONAL DE INGENIERA
2a Una masa, con dos grados 2b Una masa, con dos grados
de libertad independientes de libertad acoplados
X2 M2
K2
X1
M1
K1
-1 0 +1
0 1 2
A2
A2
A1
A1
Segundo
Modo Modo
Fundamental
62
SECCIN DE POST-GRADO - FIC
1.0
K
1.171 Xo cos 0.618 t
M
-1.0
0.4
0
K
t
-0.4 -0.171 Xo cos 1.62 M
1 M
0 1 2 3 4 5 6 7 t 8 9
2 K
1.0
Movimiento
Resultante
-1.0
63
Figura 3.4 Movimiento de la masa superior en el ejemplo 3.1
UNIVERSIDAD NACIONAL DE INGENIERA
Io
CG
L
ky y
Z Y
K ,
0.9
2 0.8
15.0 y 1.00
0.4
0.3
0.6
Io
0.5 Io 0.5 I
10.0 0.80 0.4
0.6 I
0.3
0.8
5.0 0.9 0.60
0 0.40
0 1.0 2.0 3.0 4.0 5.0
2
y
0.20
0
0 1.0 2.0 3.0 4.0 5.0
2
y
Figura 3.6 Grafico para determinar las dos frecuencias naturales acopladas
Acoplado
Acoplado
Horizontal
Superior
Cabeceo
5 10 15 20
Frecuencia cps
64
SECCIN DE POST-GRADO - FIC
13
C En C Xc = (8.5) = 0.0483 (cos 33.9t - cos 1.09 t) en pies
Yc = Y+16 () = 0.1199 cos 33.9t - 0.0199 cos 1.09 t en pies.
16 Y
En D XD = Xc
X
YD = Y = 0.0291 cos 33.9 t + 0.0309 cos 1.09 t en pies.
D
Z
0.1
0.08
0.06
0.10
0.04 0.08
0.06
Yc - Pies
0.02 0.04
0.02
xc e YD - Pies
0 0
- 0.02
- 0.02 - 0.04
- 0.06
- 0.04 - 0.08
- 0.10
YC
- 0.06 SECCIN DE POST-GRADO - FIC
- 0.08
XC YD
- 0.1
65
Figura 3.8 Movimientos en el ejemplo 3.2
UNIVERSIDAD NACIONAL DE INGENIERA
CG
T
Z P
Yo o
I II I II
Y Y Y Y
Figura 3.10 Naturaleza general de las curvas del factor de carga dinmica de
movimientos acoplados
66
SECCIN DE POST-GRADO - FIC
Debido a Debido a
movimiento cabeceo
horizontal
CG 16 5.29 8.39
3.50
Z Debido a
movimiento
vertical
Movimiento
8.5 resultante
4.63
Debido a
cabeceo
67
UNIVERSIDAD NACIONAL DE INGENIERA
3.0 3.0
2.5 2.5
II
2.0 2.0
H H 1.5
1.5
L L
1.0 1.0
I
Y
0.5 0.5
0 0
0 0.5 1.0 1.5 2.0 2.5 0 5 10 15
o. H
/y Yo
2
I
II
1
0
0 0.5 1.0 1.5 2.0
H
1
L
0
0 0.5 1.0 1.5 2.0
. H
esttico Yo
69
SECCIN DE POST-GRADO - FIC
Captulo IV
PROPAGACION DE ONDAS
4.1 INTRODUCCION
Puesto que las vibraciones transmitidas por las cimentaciones (bien sean de las estructuras
hacia el suelo, como son las fuerzas de maquinaria, o del suelo hacia las estructuras, como es el
caso de sismos) se efectan siempre a travs de ondas, es muy importante conocer los distintos
tipos de ondas que se producen en el suelo y sus mecanismos de propagacin.
Los casos sealados se pueden analizar a partir del estudio de la propagacin de ondas, tanto
en semiespacios infinitos homogneos o estratificados, as como en barras de longitud finita.
El presente captulo no pretende cubrir el estado del arte en propagacin de ondas, sino
simplemente presentar los fundamentos que se requieren para el manejo de los conceptos que se
tratan en la dinmica de suelos. Al lector que le interese profundizar ms sobre el tema, podr
consultar las referencias sealadas al final del captulo.
En primer lugar se indicarn los tipos de ondas elsticas existentes en un medio infinito,
posteriormente se analizar la propagacin de ondas en un medio semi-infinito con caractersticas
tanto homogneas como las de un medio estratificado y finalmente se describir la propagacin de
ondas en barras.
En un medio infinito, homogneo e istropo, slo se pueden propagar los dos tipos de ondas
que corresponden a las dos nicas soluciones que se obtienen de las ecuaciones de movimiento,
que ms adelante se sealan; estas dos clases de ondas son las llamadas ondas de compresin,
primarias o dilatantes y las conocidas como ondas cortantes, secundarias o distorsionales.
70
UNIVERSIDAD NACIONAL DE INGENIERA
Partiendo del anlisis de equilibrio de un elemento pequeo como el mostrado por la Figura
4.1, se llega a las siguientes expresiones conocidas en la literatura como las ecuaciones de
movimiento (los pasos para llegar a las mismas se pueden ver en la Ref. 1)
2
u
= ( + G) + G 2 u (4.1)
t 2
x
2
v
= ( + G) + G 2 v (4.2)
t 2
y
2
w
= ( + G) +G2w (4.3)
t 2
z
donde:
2 2 2
2
= + + (operador laplaciano en coordenadas cartesianas)
x 2
y 2
z2
E
G= mdulo cortante
2(1 - )
relacin de Poisson
= x + y + z dilatacin cbica
2
2 = ( + 2G) 2
t
2
= vc 2 2 (ec. de onda de dilatacin cbica) (4.4)
t 2
71
SECCIN DE POST-GRADO - FIC
donde:
+ 2G
vc = (4.5)
(1 - ) E
D = + 2G =
(1 + ) ( 1 - 2 )
Derivando ahora la ecuacin 4.2 con respecto a z y a 4.3 con respecto a y, y eliminando
mediante la substraccin de las dos expresiones resultantes, se obtiene:
w v
2
w v
- = G 2 -
t y z
2
y z
o sea
2
2x = G 2 x (4.6)
t
donde:
w v
x = 2 - , o sea es la rotacin alrededor del eje x
y z
2
x
= vs 2 2 x (4.7)
t 2
donde:
es la velocidad de las llamadas ondas cortantes o equivolumtricas y
G
vs = representa la velocidad de propagacin de la rotacin x.
Las ecuaciones correspondientes a y y z se obtienen de manera similar a la ecuacin 4.7, y
se puede decir que la rotacin se propaga con la velocidad vs.
Adems de la velocidad con que se propagan cada una de estas ondas existentes en un
medio elstico infinito, llamadas ambas ondas de cuerpo, tienen la siguiente particularidad: en las
ondas compresionales el movimiento de las partculas tiene la misma direccin en que se propagan
(vase Figura 4.2), mientras que en las ondas cortantes los movimientos de las partculas son
perpendiculares a la direccin de su propagacin. La relacin entre las velocidades de estas dos
clases de ondas est dada por la expresin:
72
UNIVERSIDAD NACIONAL DE INGENIERA
vc 2(1 - )
= (4.8)
vs (1- 2 )
la cual implica que vc > vs para cualquier valor de , y que para = 0.5, vc adquiere un valor
terico de infinito.
En un medio semi-infinito existe una frontera que permite obtener una tercera solucin a las
ecuaciones de movimiento y as tener un tercer tipo de onda. Este tercer tipo corresponde a las
ondas superficiales llamadas de Rayleigh (en honor a quien las descubri), las cuales producen en
las partculas movimientos elpticos (Figura 4.2) y disminuyen rpidamente su amplitud con la
profundidad.
u= + y w= -
x z z x
donde y son funciones potenciales que resultan estar relacionadas respectivamente con la
dilatacin y la rotacin del medio, se obtiene, al sustituir u y w en las ecuaciones 4.1 y 4.3, las
siguientes expresiones:
2 2
2 + = ( + 2 G) (2 ) + G (2 ) (4.9)
x t z t 2 x z
y
2 2
2 - = ( + 2 G) (2 ) - G ( 2 ) (4.10)
z t x t 2 z x
+ 2G 2
2
= vc
2 2
= (4.11)
t 2
y
G 2
2
= = v s 2 2 (4.12)
t
2
Ahora bien, suponiendo una solucin del tipo de onda sinusoidal viajando en la direccin
positiva de las x, se puede escribir
= F (z) ei ( t- N x )
(4.13)
y
= G (z) ei ( t- N x )
(4.14)
73
SECCIN DE POST-GRADO - FIC
donde F(z) y G(z) son funciones que describen la variacin de la amplitud de la onda con la
profundidad, y N=2 /LR (conocido como nmero de onda); LR es la longitud de la onda generada.
Al sustituir los valores de y dados por las ecuaciones 4.13 y 4.14 dentro de las ecuaciones 4.11
y 4.12, y considerar la condicin de que la amplitud de la onda superficial tiende a cero con la
profundidad, los valores de F(z) y G(z) resultan iguales a:
2
N2 z
v 2
F(z) = A1 e
c
y
2
N2 z
v 2
G (z) = A 2 e
s
Los valores de A1 y A2 se obtienen de aplicar las condiciones de frontera relativas a que los
esfuerzos cortantes y normales en la superficie del semiespacio deben ser nulos. Aplicando dichas
condiciones se obtienen las siguientes expresiones:
2
( + 2 G ) N 2 - 2 - N 2
A1 vc
-1 = 0 (4.15)
2
A2
2 i GN N 2 - 2
vs
2
2 N2 - 2
iN
A1 vc
y 2
+1= 0 (4.16)
A2 2
2N - 2
vs
v vR 2
6 4 2 2
vR v vs
- 8 R + 24 - 16 s + 16 - 1 =0 (4.17)
vs vs vc vs vc
En la Figura 4.4 se muestra la relacin que guarda vR/vs y vc/vs para varios valores de la
relacin de Poisson ; obsrvese que vR es aproximadamente igual a vs, particularmente para
valores grandes de .
74
UNIVERSIDAD NACIONAL DE INGENIERA
2 2 2 2
N - 2N - 2
2
2 vc vs
N - 2 2 2
vc
u = A 1 N i - exp - (zN) + N x
N 2
2
N - 2
vs + 1
2
N
2
2
N - 2
vs
exp - (zN ) exp i ( t - N x) (4.18)
N
y
2
2
2 N - 2 2
2 2
2
vc N - N -
2 2
N vs - vc
w=A 1N exp - x
2
N2-
N
N
2
vs + 1
2
N
2
2
N - 2
vs
exp - (zN ) exp i ( t - N x) (4.19)
N
De la observacin de estas dos ecuaciones, se puede deducir que los trminos dentro de las
llaves representan la variacin respectiva de u y w con la profundidad. O sea:
u = U (z) A1 Ni e i ( t - Nx)
y
w = W (z) A1 N e i ( t - Nx)
La variacin de U(z) y W(z) con la profundidad para varios valores de , se indica en la Figura
4.5. La Figura 4.6 seala la interpretacin fsica del concepto de longitud de onda que interviene en
la figura anterior.
Debe sealarse que son las ondas Rayleigh las que trasmiten la mayor parte de la energa
generada por la vibracin de una zapata sobre la superficie de un semiespacio. (Cuando la zapata
es circular, el 67% de la energa es trasmitida por las ondas Rayleigh, mientras que las cortantes
trasmiten el 26% y las de compresin el 7% restante). Por otro lado, en comparacin con las ondas
de cuerpo, las amplitudes de las ondas Rayleigh disminuyen ms lentamente con la distancia r al
centro de la fuente de excitacin, mientras que la atenuacin de las ondas P y S en la superficie es
proporcional a 1/r, en las ondas Rayleigh es proporcional a 1 / r , la razn de esta diferencia se
75
SECCIN DE POST-GRADO - FIC
debe al concepto del llamado amortiguamiento radial que se estudia en el siguiente captulo. Lo
anterior hace, como se ilustrar posteriormente, que las ondas Rayleigh desempeen un papel muy
importante en la trasmisin de vibraciones en o cerca de la superficie.
Las ondas Rayleigh son generalmente fciles de reconocer ya que usualmente tienen una
amplitud grande con frecuencia relativamente baja, segn puede observarse en la Figura 4.7.
Sin embargo, con el objeto de considerar la divisin de la energa que se origina en el punto de
incidencia, es conveniente considerar primeramente el caso particular de la descomposicin de las
ondas P y S al llegar a una superficie libre. Para ello resulta a la vez conveniente tomar en cuenta
que las ondas cortantes S se pueden descomponer en una componente paralela a la superficie
(ondas SH), y en otra contenida en el plano vertical (ondas SV). La Figura 4.8 ilustra esta
descomposicin.
Cuando una onda dilatante P incide sobre la superficie libre del semiespacio, parte de la
energa se refleja a travs de una onda cortante SV y parte a travs de una onda P (Figura 4.9). El
ngulo de reflexin 1 de la onda SV est dado de acuerdo con la ley de Snell
vs
sen 1 = sen
vp
Al llegar una onda cortante SV a la superficie, toda la energa que se refleja se hace a travs
de: a) una onda SV con un ngulo de reflexin igual al de incidencia (Figura 4.10), y b) a travs de
una onda P cuyo ngulo de generacin est dado por:
vc
sen 1 = sen
vs
Existe un cierto ngulo de incidencia, llamado crtico, para el cual las ondas incidentes P y S se
reflejan horizontalmente (Figura 4.11); dicho ngulo depende nicamente de la relacin de Poisson.
76
UNIVERSIDAD NACIONAL DE INGENIERA
vc
Para onda dilatantes cr = sen -1 , y
vs
vs
Para ondas cortantes cr = sen -1
vp
La Figura 4.12 muestra la relacin entre cr y para el caso de ondas de incidencia SV.
Cuando los ngulos de incidencia son mayores, las componentes horizontal y vertical de los
movimientos del terreno se encuentran desfasadas, creando una vibracin del tipo elipsoidal; la
Figura 4.13 muestra que para s = 45 el movimiento es vertical y que para s =90 el movimiento se
reduce a cero.
En el caso de una onda SH que llega a la superficie, toda la energa que se refleja se hace a
travs de otra onda SH, la cual tiene un ngulo de reflexin igual al de incidencia (Figura 4.14). Esta
caracterstica hace que existan procedimientos especiales por medio de los cuales se generen este
tipo de ondas y se facilite la interpretacin de los datos obtenidos mediante los mtodos
geossmicos; el empleo de dichos mtodos se explicar en captulos posteriores.
Ahora bien, para el caso de llegar una onda a la superficie de contacto de dos estratos de
caractersticas diferentes, se tendr lo siguiente:
Al llegar una onda P sobre la superficie de contacto, se producen cuatro tipos de ondas segn
se ilustra en la Figura 4.15a; dos ondas SV (una reflejada [P-SV1] y otra refractada [P-SV2]) y dos
ondas P (una reflejada [P-P1] y otra refractada [P-P2]).
En cuanto a las ondas incidentes SH, parte de la energa es reflejada (ondas SH-SH1) y parte
refractada (SH-SH2), pero nuevamente slo a travs de ondas SH; la razn de no producir ondas P
se debe a que las ondas SH no tienen componente normal en el plano de contacto.
Los ngulos de reflexin y refraccin pueden calcularse a partir de la ley de Snell, de la cual se
obtiene la siguiente expresin:
donde:
vp1 y vp2 son respectivamente las velocidades de las ondas dilatantes en los medios
superior e inferior, y anlogamente.
vs1 y vs2 son las velocidades de las ondas cortantes de dichos medios.
77
SECCIN DE POST-GRADO - FIC
Cuando la velocidad de una onda reflejada o refractada es mayor que la velocidad de la onda
incidente, puede haber un ngulo de incidencia crtico para el cual la onda reflejada o refractada
ser horizontal; dicho ngulo se obtiene a partir de las expresiones 4.20. Por ejemplo, para una
onda dilatante P incidente,
vp 1
cr = sen -1
vp 2
Existen en la literatura frmulas y grficos que proporcionan la cantidad de energa que se
trasmite a travs de cada una de las ondas reflejadas o refractadas; vanse por ejemplo las
Referencias 1 y 2.
Cuando existen varios estratos se tendrn mltiples refracciones y reflexiones, segn puede
observarse en la Figura 4.16, y el problema de propagacin de ondas se vuelve ms complejo.
Cuando el estrato superior es menos rgido que el que lo subyace, se puede generar otro tipo de
onda superficial conocida como onda Love; estas ondas son originadas por las reflexiones totales
mltiples de la capa superior, y son ondas que se desplazan horizontalmente y producen
movimientos transversales horizontales. Ewing (Ref. 3) define a esta clase de onda como "la onda
cortante polarizada horizontalmente, atrapada en la capa superficial y originada por las reflexiones
totales mltiples". Jones (Ref. 4) demuestra que para altas frecuencias de excitacin, la velocidad de
propagacin de las ondas Love se aproxima asintticamente a la velocidad de propagacin de las
ondas cortantes en el estrato superior, mientras que para bajas frecuencias dicha aproximacin se
refiere a la velocidad de las ondas cortantes en el estrato inferior.
Cuando las ondas dilatantes o compresionales se propagan en medios que no son infinitos, las
condiciones de frontera modifican las ondas generadas haciendo que stas sean un poco diferentes
a las sealadas hasta ahora. Por ejemplo, las ondas compresionales que se propagan a travs de
una barra donde pueden haber expansiones
libres en el sentido transversal, tienen una velocidad de propagacin que resulta, segn se
demuestra ms adelante, aproximadamente igual a:
E
vL = (4.21)
Esta velocidad es menor que la velocidad vc dada por la ecuacin 4.5; la razn de ello es que
en un medio infinito o semi-infinito no existen desplazamientos normales a la direccin en que se
propagan estas ondas, mientras que en una barra dichos desplazamientos son factibles. A esta
clase de ondas compresionales en barras se les conoce en la literatura con el nombre de ondas
longitudinales.
La obtencin de la ecuacin 4.21 se puede hacer a partir del anlisis de fuerzas actuando en un
elemento de barra de longitud x (Figura 4.17), que tiene una seccin transversal de rea A.
78
UNIVERSIDAD NACIONAL DE INGENIERA
2 u
x A = x A -
2
x t
Simplificando la expresin anterior se obtiene
2
u
+ =0 (4.22)
x t 2
Esta misma ecuacin se puede expresar en otros trminos, de la siguiente manera. De la teora
de elasticidad se tiene:
= E (4.23)
u
donde: =-
x
2 2
u u
E = (4.24)
x 2
t 2
que es la llamada "ecuacin de ondas en una dimensin". La solucin a esta ecuacin es del
tipo
E
u = f x t (4.25)
E
u = sen x t
E
u = cos x t
E
u = x t
El significado fsico de las implicancias de dicha solucin se muestra en la Figura 4.18. Para
un tiempo cualquiera t1 (que puede ser t1=0), se tiene un cierto tipo de desplazamiento caracterizado
por una funcin que satisfaga la ecuacin 4.25; posteriormente, en el tiempo t2, se observar
exactamente el mismo tipo de desplazamiento pero en un lugar diferente. Es decir, el tipo de
movimiento que se observa es precisamente como el de una onda que se desplaza a una velocidad
E /
vL= . Analticamente lo anterior se puede comprobar de la siguiente manera; supngase el
signo negativo de la ecuacin 4.25, y que t2 = t1 + t; se tiene entonces:
79
SECCIN DE POST-GRADO - FIC
u |t = t1 = f (x - vL t )
u |t 2 = t1+ t = f [ (x + x) - vL (t + t ) ]
u |t 2 = t 1+ t = f [ (x + vL t - vL t - v L t ) ] = f (x - vL t )
u
x=E
x
de donde se obtiene que
u=
x x =x t
vL
E E
. u x vL
u= = (4.26)
t E
Obsrvese en esta ltima expresin que la velocidad de la partcula depende del valor del
esfuerzo aplicado, mientras que la velocidad de propagacin de ondas depende slo de las
propiedades del material.
Ahora bien, al analizar las ondas cortantes en barras, siguiendo un procedimiento similar al
descrito para las ondas compresionales, se llega a que la ecuacin de onda est dada por la
siguiente expresin:
2 2
= vs 2 (4.27)
t 2
x 2
G
vs = (4.28)
es la velocidad con la que se propagan las ondas cortantes en barras.
Este valor, como puede notarse, resulta ser igual al obtenido en el anlisis de propagacin de
ondas en un medio infinito o semi-infinito.
80
UNIVERSIDAD NACIONAL DE INGENIERA
Obsrvese que conociendo las velocidades vL vs, los mdulos E y G se pueden obtener
respectivamente mediante las ecuaciones 4.21 y 4.28.
donde:
2 U n2
+ U=0 (4.30)
x 2 vL 2
U = A3 cos
n x + n x
A 4 sen
vL vL
cos
n l = o
vL
de donde se deduce que
vL
n = (2 n- 1) , n = 1,2,3... (4.31)
2l
81
SECCIN DE POST-GRADO - FIC
(2 n- 1) x
U = A 4 sen (4.32)
2 l
En la Figura 4.19 se muestran los tres primeros modos de vibracin de una probeta circular y
el significado fsico de la constante A4. Al sustituir la ecuacin 4.32 en 4.29, se obtiene la forma
general de los desplazamientos:
(2 n - 1) x (2 n- 1) vL t (2 n- 1) vL t
u = sen (A 1 ) n cos + (A 2 ) n sen (4.33)
2 l 2 l 2 l
Para otras condiciones de frontera o para el caso de vibraciones torsionales se podr seguir
el procedimiento descrito y obtener expresiones anlogas a la ecuacin 2.33. La expresin
correspondiente a la frecuencia circular natural bajo excitaciones torsionales, considerando las
mismas condiciones de frontera (un extremo fijo y el otro libre), resulta exactamente la misma que la
dada por la ecuacin 4.31, slo que en vez de vL interviene vs.
REFERENCIAS
1) Richart, F.E., Hall, J.R., y Woods, R.D. (1970), "Vibrations of Soils and Foundations",
Prentice-Hall.
3) Ewing, W.M., Jardetsky, W.S., y Press, F. (1975), "Elastic Waves in Layered Media",
McGraw-Hill Book Co, New York.
82
UNIVERSIDAD NACIONAL DE INGENIERA
xz y
( xz + z)
z
x
xy z
xy
XZ
x ( xy + y)
y
x x
( x + x )
x
83
SECCIN DE POST-GRADO - FIC
(a)
(b)
(c)
(d)
84
UNIVERSIDAD NACIONAL DE INGENIERA
Superficie
Frente de ondas
5
G
=
VS
V
4
Valores de
Ondas P
Ondas S
Ondas R
0
0 0.1 0.2 0.3 0.4 0.5
Relacin de Poisson,
Figura 4.4 Relacin entre VS, VP, y VR, contra la relacin de Poisson
85
SECCIN DE POST-GRADO - FIC
0.0
0.2
0.6
LR
z
= 0.25 = 0.25
= 0.33
Longitud de onda
= 0.33
Profundidad
= 0.40
0.8
= 0.50 = 0.40
= 0.50
1.0
1.2
Componente vertical W (z) [ ]
1.4
86
UNIVERSIDAD NACIONAL DE INGENIERA
Amplitud
Tiempo
u
Onda P Onda S Onda R
(+ en la direccin de propagacin)
t
(a) 1
(+ hacia abajo)
1
t
(b)
87
SECCIN DE POST-GRADO - FIC
SV Plano vertical
de incidencia
S
E
SH
Plano perpendicular
Rayo incidente
al rayo incidente
P
1
SV 1
P P
SV SV
Figura 4.9 Reflexin en la superficie de Figura 4.10 Reflexin de una onda incidente
una onda incidente P SV en una superficie libre
88
UNIVERSIDAD NACIONAL DE INGENIERA
cr cr
SV SV
40
cr
30
20
Para ondas SV
10
0
0.25 0.30 0.35 0.40 0.45 0.50
89
SECCIN DE POST-GRADO - FIC
d d
d d
Superficie
a a a
a
= 30 = 34
= 0 = 20
d Superficie
a
= 30 16
d d d d
Superficie
a a a a
= 40 = 45
= 37
= 37 1/2
d d d d d=0
a Sup. a a
a
a
= 85 = 90
= 50 = 63 = 75
90
UNIVERSIDAD NACIONAL DE INGENIERA
SH SH
Medio 2
P2, vP2 , vS2
e
e
f f
P-P2 SV-P2
P-SV2 SH-SH2
f SV-SV2
Punto de excitacin
(P
1)
(P - S
)
- P2
)
)
+ x
x
Area A
x
x
VL (t2 t1)
Figura 4.18 Desplazamientos observados en los tiempos t1 y t2, para una funcin
del tipo sealado por la ec. 4 - 5
(a)
A4 x
u1 = A 4 sen (n = 1)
2l
3 x
A4 u2 = A 4 sen (n = 3)
(b) 2l
5 x
u3 = A 4 sen (n = 5)
2l
A4
l
Figura 4.19 Primeros tres modos naturales de vibracin de una
barra con un extremo fijo y el otro libre 92
SECCIN DE POST-GRADO - FIC
Captulo V
DDISTRIBUCIN DE MXIMAS INTESIDADES
SSMICAS OBSERVADAS EN EL PER Jorge E. Alva Hurtado (1)
RESUMEN
INTRODUCCIN
La primera etapa del Proyecto SISRA (Mitigacin de Daos Causados por un Terremoto en
la Regin Andina), patrocinado por el Centro Regional de Sismologa para Amrica del Sur
(CERESIS), consisti en recopilar y evaluar los datos histricos y geolgicos referentes a la
sismicidad andina y la informacin instrumental de eventos ssmicos, para estimar el peligro
ssmico de la regin. El compromiso mnimo del Grupo de Peligro Ssmico del Proyecto SISRA era
confeccionar a nivel regional un mapa de Intensidades Mximas en la escala Mercalli Modificada.
Tambin se prepararan separadamente mapas de Distribucin de Areas de Licuacin de Suelos y
de Deslizamientos Producidos por Sismos en la Regin Andina. En este trabajo se presenta el
mapa de Distribucin de Mximas Intensidades Ssmicas Observadas en el Per.
(1) Profesor Asociado, Facultad de Ingeniera Civil, Universidad Nacional de Ingeniera, Lima, Per.
(2) Ingeniero Civil, Universidad Nacional de Ingeniera, Lima, Per.
(3) Ingeniero Civil, ELECTROPERU S.A. Lima, Per.
Ponencia presentada en el V Congreso Nacional de Ingeniera Civil, Tacna, 12-18 Noviembre de 1984.
93
UNIVERSIDAD NACIONAL DE INGENIERA
En este mapa se expresan las intensidades mximas por medio de curvas isosistas en la
escala Mercalli Modificada, incluyendo eventos histricos de importancia ocurridos en el Per
hasta el 31 de Diciembre de 1981. La principal diferencia de la versin reciente del mapa de
Intensidades Mximas con la versin anterior (Alva Hurtado y Torres Cabrejos, 1983) es la
incorporacin de tres mapas de isosistas de sismos histricos preparados por Silgado (1983) y la
revisin actualizada de la informacin disponible. El mapa presentado para el Per ha sido
compatibilizado con los mapas de intensidades mximas de los pases vecinos por el Coordinador
Regional del Grupo de Peligro Ssmico del Proyecto SISRA, Ing. Juan Carlos Castano, del
Instituto Nacional de Prevencin Ssmica (INPRES) de la Repblica Argentina.
DOCUMENTACIN BSICA
METODOLOGA UTILIZADA
94
SECCIN DE POST-GRADO -
Con el objeto de documentar toda modificacin realizada por los autores a los valores de
intensidades publicadas, se presenta una lista que resume dichas modificaciones:
Sismo del 24 de Julio de 1912. Se disminuyeron los valores de intensidades, hasta valores
mximos de VIII MM y el epicentro se situ al oeste del mapa original (Silgado, 1977).
95
UNIVERSIDAD NACIONAL DE INGENIERA
El mapa presentado a escala 1:5'000,000 en este artculo fue dibujado teniendo como base el
mapa original a escala 1:2'000,000, por lo que pueden haberse incorporado algunas imprecisiones
en la localizacin de las curvas.
CONCLUSIONES Y RECOMENDACIONES
96
SECCIN DE POST-GRADO - FIC
Existe suficiente informacin de intensidades ssmicas del pasado en la Costa, Sierra y Selva
Alta del Per para poder elaborar un mapa de distribucin de intensidades. Sin embargo, dicha
informacin es escasa en la zona de Selva Baja.
Adems de la alta sismicidad existente en la Costa del Per, se aprecia una gran actividad
ssmica en la denominada Zona Subandina, localizada en la Selva Alta.
Existe alguna similitud entre el mapa presentado y el mapa de Zonificacin Ssmica del Per
de las Normas de Diseo Sismo-Resistente en vigencia (RNC, 1977). La zona 3 (sismicidad baja)
coincide con la zona de intensidades menores que IV MM en el mapa. Como se indic
anteriormente, existe una gran escasez de informacin en esta rea. La zona 2 (sismicidad media)
se asemeja a la zona de intensidad V MM, an cuando existen diferencias en los departamentos
de Huancavelica, Ayacucho, Cuzco y Puno. El mapa presentado indica en dichas zonas mayor
sismicidad. La zona 1 (sismicidad alta) est representada en el mapa presentado por reas de
intensidad mayor que VI MM. Esta zona presenta un gran rango de valores de intensidad, por lo
que podra subdividirse.
AGRADECIMIENTOS
Se agradece a los Drs. Enrique Silgado y Leonidas Ocola por haber puesto a disposicin
de los autores las publicaciones relevantes. Igualmente se agradece el apoyo brindado a este
Proyecto por el Ing. Alberto Giesecke, Director del Centro Regional de Sismologa para Amrica
del Sur. Se agradece a los Ings. Edgar Arvalo, Abel Ordoez y Luis Chang, ex-alumnos del
primer autor por su colaboracin en la bsqueda de la informacin pertinente; de igual manera se
agradece a la Srta. Elizabeth Zegarra por el paciente mecanografiado.
REFERENCIAS
1. Alva Hurtado, J.E. (1981) "Bibliografa sobre los Terremotos Peruanos: Aspectos Ingenieriles",
Revista el Ingeniero Civil, Ao 3, No. 15, pp. 50-55.
97
UNIVERSIDAD NACIONAL DE INGENIERA
4. Deza, E. (1971), "The Pariahuanca Earthquakes Huancayo, Peru: July-October 1969", Recent
Crustal Movements, Royal Society of New Zealand, Bulletin 9, pp. 77-83.
6. Espinosa, A.F., Husid, R., Algermissen, S.T. y De las Casas, J. (1977), "The Lima Earthquake
of October 3, 1974: Intensity Distribution", Bulletin of the Seismological Society of Amrica, Vol.
67, No. 5, pp. 1429-1439.
7. Giesecke. A., Ocola, L. y Silgado, E. (1980), "El Terremoto de Lima del 3 de Octubre de 1974",
Informe preparado por el Centro Regional de Sismologa para Amrica del Sur a UNESCO,
Contrato SC/RP-601.013.
8. Guzmn, V. (1983), "El Riesgo Ssmico en Proyectos de Ingeniera Civil", Tesis de Grado,
Programa Acadmico de Ingeniera Civil, Universidad Nacional de Ingeniera, Lima, Per.
9. Hansen, W. (1968), "Earthquake Activity Near the Tarapoto Rio Nieva Highway", Woodward-
Clyde and Associates, Report to ENKAY S.A.
10. Herrera Cano, J. y Giuliani, H. (1975), "El Terremoto de Lima (Per) del 3 de Octubre de 1974 -
Espectro de la Componente N 82W del Acelerograma", Informe preparado por el Centro
Regional de Sismologa para Amrica del Sur a UNESCO, Contrato SC/RP - 601.013.
11. Huaco, P., Lescano, W., Bueno, A., Vsquez, M. y Chvez, J. (1974), "Evaluacin de
Intensidades, Sismo del 3 de Octubre de 1974 en la zona de Lima-Ro Pisco", Instituto
Geofsico del Per, Lima, Per.
12. Huaco, D., Rodrguez, A. y Rodrguez, L., "Sismicidad de Arequipa", Informe Indito, Instituto
Geofsico del Per, Lima, Per.
13. IGP (1974), "Boletn del Instituto Geofsico del Per del Sismo del 3 de Octubre de 1974",
Instituto Geofsico del Per, 18 de Octubre de 1974.
14. Kuroiwa, J. y Deza, E. (1968), "Daos Causados en Moyobamba por el Sismo del 19 de Junio
de 1968", Facultad de Ingeniera Civil, Universidad Nacional de Ingeniera, Lima, Per.
15. Ministerio de Agricultura y Alimentacin (1978), "Desarrollo Integral de los Recursos Hdricos
de Tacna y Moquegua", Vol. IV, Estudios Bsicos Generales, Apndice 6, Sismologa, Oficina
de Estudios, Estudios Integral de las Cuencas de Tacna y Moquegua, Diciembre.
17. Ocola, L. (1982), "Mapa Preliminar Intensidades Mximas (MSK) 1471-1974", Proyecto SISAN,
Instituto Geofsico del Per, Lima, Per.
98
SECCIN DE POST-GRADO - FIC
18. Orihuela, P. (1981), "Influencia de las Condiciones Locales en el Sismo de Arequipa 16.02.79",
Tesis de Grado, Programa Acadmico de Ingeniera Civil, Universidad Nacional de Ingeniera,
Lima, Per.
21. Shikiya, J. y Alemn, H. (1974), "Evaluacin de Intensidades y Daos producidos por el Sismo
del 3 de Octubre de 1974 en la Zona de San Juan", Institu to Geofsico del Per, Lima, Per.
100
22. Sieberg, A. (1930), "Los Terremotos en el Per", Captulo VI, Geologa del Per de G.
Steinmann, Heidelberg Carl Winster Universitat-Shuchhan-Lunc, p. 406.
23. Silgado, E. (1946), "Datos Sismolgicos del Per", Boletn No. 3, Instituto Geolgico del Per,
Lima, Per.
24. Silgado, E. (1947), "Datos Sismolgicos del Per-1946", Boletn No. 7, Instituto Geolgico del
Per, Lima, Per.
25. Silgado, E. (1948a), "Datos Sismolgicos del Per-1947", Boletn No. 11, Instituto Geolgico del
Per, Lima, Per.
26. Silgado, E. (1948b), "Datos Sismolgicos del Per-1948", Boletn No. 13, Instituto Geolgico del
Per, Lima, Per.
27. Silgado, E. (1951), "The Ancash Earthquake of November 10, 1946", Bulletin of the
Seismological Society of America, Vol. 41, No. 2, pp. 83-99.
28. Silgado, E. (1952a), "Estadstica Ssmica del Per, 1949-1950", Datos Sismolgicos del Per
1949-1950, Boletn No. 4, Instituto Nacional de Investigacin y Fomento Mineros, Ministerio de
Fomento y Obras Pblicas, Lima, Per, pp. 7-26.
29. Silgado, E. (1952b), "El Sismo del 9 de Diciembre de 1950", Datos Sismolgicos del Per 1949-
1950, Boletn No. 4, Instituto Nacional de Investigacin y Fomento Mineros, Ministerio de
Fomento y Obras Pblicas, Lima, Per, pp. 47-51.
30. Silgado, E. (1953), "Datos Sismolgicos del Per-1951", Boletn No. 8, Instituto Nacional de
Investigacin y Fomento Mineros, Ministerio de Fomento y Obras Pblicas, Lima, Per.
31. Silgado, E. (1957a), "El Movimiento Ssmico del 12 de Diciembre de 1953", Anales del Primer
Congreso Nacional de Geologa, Sociedad Geolgica del Per, Lima, Per, Parte II, Tomo 32,
pp. 225-238.
99
UNIVERSIDAD NACIONAL DE INGENIERA
32. Silgado, E. (1957b), "Datos Sismolgicos del Per 1952-1955", Boletn de la Sociedad
Geolgica del Per, Lima, Per, Tomo XXIX.
33. Silgado, E. (1973), "Historia de los Sismos ms Notables Ocurridos en el Per (1513-1970)",
Geofsica Panamericana, Vol. 2, No.1, pp. 179-243.
34. Silgado, E. (1977), "Datos Macrossmicos de Terremotos Ocurridos en Territorio Peruano 1912-
1974. Mapas de Isosistas y Apreciaciones", Centro Regional de Sismologa para Amrica del
Sur, Lima, Per.
35. Silgado, E. (1978), "Historia de los Sismos ms Notables Ocurridos en el Per (1513-1974)",
Instituto de Geologa y Minera, Boletn No. 3, Serie C, Geodinmica e Ingeniera Geolgica,
Lima, Per.
36. Silgado, E. (1983), "Isosistas de los Sismos del 14 de Febrero de 1619, 28 de Octubre de 1746
y el 13 de Agosto de 1868", Comunicacin Personal.
37. Silgado, E., Fernndez-Concha, J. y Ericksen, G.E. (1952), "El Terremoto del Cuzco del 21 de
Mayo de 1950", Datos Sismolgicos del Per 1949-1950, Boletn No. 4, Instituto Nacional de
Investigacin y Fomento Mineros, Ministerio de Fomento y Obras Pblicas, Lima, Per, pp. 27-
46.
38. Silgado, E. y De las Casas, J. (1971), "Seismological Report on Chira Piura Project",
Submitted to Energoprojekt, Lima, Per.
39. Umlauff, A.F. (1915), "La Regin Ssmica de Caravel", Boletn de la Sociedad Geogrfica de
Lima, Tomo XXXI, Segundo Trimestre, Lima, Per.
100
SECCIN DE POST-GRADO - FIC
101
UNIVERSIDAD NACIONAL DE INGENIERA
102
SECCIN DE POST-GRADO - FIC
Captulo VI
BREVE HISTORIA DEL FENMENO DE
LICUACIN DE SUELOS EN EL PER
Dr. Ing. Jorge E. Alva Hurtado (1)
RESUMEN
INTRODUCCIN
El Grupo de Evaluacin del Peligro Ssmico del Proyecto SISRA (Sismicidad de la Regin
Andina), patrocinado por el Centro Regional de Sismologa para Amrica del Sur (CERESIS), se
fij como objetivo obtener un producto final que expresara de alguna manera niveles relativos del
peligro ssmico y que fuera homogneo para la regin andina. Se tom en cuenta el desarrollo
alcanzado en cada pas en esta clase de estudios y el tiempo disponible para ejecutar el trabajo.
El resultado final en comn con los pases andinos consisti en confeccionar mapas regionales a
partir de los antecedentes histricos, graficando intensidades mximas, reas de licuacin de
suelos y zonas de deslizamiento asociadas a terremotos.
103
UNIVERSIDAD NACIONAL DE INGENIERA
DOCUMENTACIN BSICA
22 DE ENERO DE 1582
Silgado (1978), basado en la obra de Fray Vctor Barriga, "Los Terremotos de Arequipa",
indica un terremoto que dej en ruinas a la ciudad de Arequipa. Despus del mismo toda la ciudad
qued anegada por la mucha agua que corri por sus calles, probablemente debido al
asentamiento y compresin de la napa fretica. En la costa se sumergi un arroyo de agua que
entraba al Puerto de Islay. Tambin se referencia este fenmeno en la Historia del Per del Padre
Rubn Vargas Ugarte, S.J. La mxima intensidad de este sismo fue de X MM.
14 DE FEBRERO DE 1619
Terremoto en el norte del Per, que arruin los edificios de Trujillo. Silgado (1978) refiere la
crnica del Padre Calancha que indica que se agriet la tierra en varias partes, de las cuales
surgi un lodo negruzco. Adems Feijoo deca que el material viscoso y pestilente expelido
envuelto en agua gredosa de que se formaron ros que corrieron por las campias de la Villa del
Santa, el de Barranca y otros, tieron sus corrientes con tal maligna inundacin. De los relatos se
desprende que hubo otros efectos en las vertientes, en algunos lechos secos de ros apareci
104
SECCIN DE POST-GRADO - FIC
agua, mientras que en otros riachuelos el agua dej de correr. La mxima intensidad de este
sismo fue de IX MM.
31 DE MARZO DE 1650
Terremoto en el Cuzco que derrib todos los templos y la mayor parte de las edificaciones.
Segn Silgado (1978), la tierra se agriet en varios lugares, observndose disturbios en el nivel
fretico de las aguas de escorrenta cerca del pueblo de Oropesa. La mxima intensidad de este
sismo fue de IX MM.
12 DE MAYO DE 1664
Terremoto en Ica. Silgado (1978) refiere el relato del cura de la Iglesia Parroquial de San
Jernimo, Licenciado Cristbal Rodrguez Alvarez: "otra vez abrise la tierra por muchas partes,
en los montes y campos se desuna la tierra formando abras y horribles profundidades que
parecan bocas para tragarnos. Corri el ro en ms de seis riegos de agua, rebosaron algunos
pozos de la ciudad, arrancndose de raz muchos y grandsimos rboles, sauces y espinos". La
mxima intensidad de este sismo fue de X MM.
10 DE FEBRERO DE 1716
Terremoto en Pisco que La Barbinais relata: "al cuarto de hora tembl de nuevo la tierra, que
abrindose en algunos lugares, expeli chorros de polvo y agua con ruido pavoroso" (Silgado,
1978). La mxima intensidad de este sismo fue de IX MM.
28 DE OCTUBRE DE 1746
Terremoto en Lima y tsunami en el Callao. Bachmann (1935) indica que en la quebrada
mineral del ro Viseca, de la Provincia de Lucanas, se abri la tierra y salieron sabandijas;
reventando tambin un volcn de agua caliente que lo inund todo. En las Misiones de
Cajamarquilla de los Franciscanos, en Pataz, reventaron dos volcanes de cieno y lodo. La
mxima intensidad de este sismo fue de X MM.
1747
Toribio Polo en su "Sinopsis de los Temblores y Volcanes del Per" seala que en dicho
ao ocurri un fuerte sismo, que por haberse producido en una regin apartada de Puno no se
pudo obtener fecha. Se sabe que ocasion grandes destrozos en Ayapata, provincia de Carabaya;
agua cenagosa brot de la tierra y pereci mucha gente (Silgado, 1978). La mxima intensidad de
este sismo fue de VIII MM.
30 DE MARZO DE 1813
Silgado (1978) refiere un terremoto en Ica que destruy casas y templos, muriendo 32
personas. Se formaron grandes grietas en el cauce del ro, del cual surgi gran cantidad de lodo.
La mxima intensidad de este sismo fue de VII MM.
20 DE AGOSTO DE 1857
Silgado (1978) indica un fuerte sismo en Piura que destruy muchos edificios. Se abri la
tierra, de la cual emanaron aguas negras. Daos menores en el puerto de Paita. La mxima
intensidad de este sismo fue de VIII MM.
105
UNIVERSIDAD NACIONAL DE INGENIERA
13 DE AGOSTO DE 1868
Terremoto acompaado de tsunami en Arica. Silgado (1978) refiere el relato de Toribio
Polo: "Agrietamientos del suelo se observaron en varios lugares, especialmente en Arica, de los
que brot agua cenagosa". Bachmann (1935) reporta que en Sama y Locumba se perdi gran
parte de las cosechas y la tierra se abri a trechos en hondas grietas que vomitaban agua
cenagosa. La mxima intensidad de este sismo fue de XI MM.
24 DE JULIO DE 1912
Terremoto en Piura y Huancabamba. En el cauce seco del ro Piura se formaron grietas
con surgencia de agua, otros daos afectaron el terrapln del ferrocarril. En el puerto de Paita se
produjeron agrietamientos del suelo (Silgado, 1978). La mxima intensidad de este sismo fue de
VIII MM.
24 DE DICIEMBRE DE 1937
Terremoto en las vertientes orientales de la Cordillera Central. Afect los pueblos de
Huancabamba y Oxapampa. Silgado (1978) indica que en el Fundo Victoria se abri una grieta de
la que eman abundante cantidad de agua que arras corpulentos rboles, aumentando el caudal
del ro Chorobamba. La mxima intensidad de este sismo fue de IX MM y la magnitud fue de Ms =
6.3.
24 DE MAYO DE 1940
Terremoto en la ciudad de Lima y poblaciones cercanas. Valencia (1940) reporta que en el
Callao quedaron efectos del sismo, sobretodo en terrenos formados por relleno hidrulico. En
estas zonas el terreno se agriet y brot a la superficie masas de lodo semilquido. Las grietas del
terreno atravezaron algunas construcciones. La mxima intensidad de este sismo fue de IX MM y
su magnitud fue de Ms = 8.0.
6 DE AGOSTO DE 1945
Fuerte temblor en la ciudad de Moyobamba y alrededores. De acuerdo a Silgado (1946),
se formaron algunas grietas en la quebrada de Shango. Posteriormente, el temblor del da 8
produjo nuevas grietas vecinas a las primeras, una de ellas semicircular de 15 m. de dimetro y 4
cm. de separacin, de las cuales emanaron aguas cargadas de limo durante dos das. Las grietas
se presentaron tambin en los bordes de los barrancos en Tahuisco, cerca del ro Mayo y en la
quebrada Azungue. A unos cinco km. de los baos sulfurosos y a diez km. de la ciudad se haba
producido la aparicin de nuevos manantiales. La mxima intensidad de este sismo fue de VII
MM.
28 DE MAYO DE 1948
Fuerte sismo destructor en Caete. En las inmediaciones del lugar denominado Calavera
se produjeron varios deslizamientos en terrenos pantanosos. En las faldas del Cerro Candela se
formaron grietas, observndose en el lugar pequeos derrumbes debido a la saturacin del
terreno (Silgado, 1978). La mxima intensidad de este sismo fue de VII MM y su magnitud fue de
Ms = 7.0.
106
SECCIN DE POST-GRADO - FIC
21 DE MAYO DE 1950
Terremoto en la ciudad del Cuzco. Silgado, Fernndez-Concha y Ericksen (1952) notaron
en el lado sur del Valle, al sureste del pueblo de San Sebastin, una zona de extensa fisuracin.
Tambin observaron dos pequeas fracturas en una zona pantanosa situada a 300 m. al sur de
San Sebastin, de las cuales surgi agua y arena durante el terremoto. Los hoyos producidos por
la eyeccin tenan cerca de 2 m. de dimetro y la arena alrededor de la fractura un espesor de 1 a
2 cm. Durante el movimiento ssmico estas fracturas y otras producidas a lo largo del cerro,
vertieron chorros de agua que alcanzaron 1 a 2 m. de altura. El nivel de la napa fretica se levant
en el lado sur del Valle. Areas que haban estado casi secas antes del terremoto, aparecieron
cubiertas con 10 a 40 cm. de agua, semana y media despus del sismo. El agua en un pozo de la
Hacienda San Antonio subi a 1.80 m. por encima de su nivel normal, despus del terremoto. La
mxima intensidad de este sismo fue de VII MM y su magnitud fue de Ms = 6.0.
9 DE DICIEMBRE DE 1950
Fuerte temblor en Ica. En el Fundo La Vela se produjeron algunas pequeas grietas en el
terreno de sembro, de los cuales se dice, sali agua hasta unas horas despus del sismo
(Silgado, 1951). La mxima intensidad del sismo fue de VII MM y su magnitud de Ms = 7.0.
12 DE DICIEMBRE DE 1953
Un fuerte y prolongado movimiento ssmico afect la parte noroeste del Per y parte del
territorio ecuatoriano. Silgado (1957) indic que se produjeron grietas largas en los terrenos
hmedos. Se apreciaron eyecciones de lodo en la quebrada de Bocapn, en los esteros de Puerto
Pizarro y en otros lugares. En Bocapn, que haba estado seco antes del movimiento, corri
momentneamente agua a causa de los surtidores. En Puerto Pizarro se originaron chorros de
agua de 60 cm. de altura y grietas. La mxima intensidad de este sismo fue de VIII MM y su
magnitud fue de Ms = 7.8.
15 DE ENERO DE 1958
Terremoto en Arequipa. Silgado (1978) indic agrietamiento del terreno cerca de la zona
de Caman, con eyeccin de aguas negras. La mxima intensidad del sismo fue de VIII MM y su
magnitud fue de Ms = 7.0.
17 DE OCTUBRE DE 1966
La ciudad de Lima fue estremecida por un sismo. En la Hacienda San Nicols, a 156 km. al
norte de Lima, aparecieron numerosas grietas y de varias de ellas surgi agua de color amarillo
(Silgado, 1978). La mxima intensidad de este sismo fue de VIII MM y su magnitud fue de mb =
6.3.
19 DE JUNIO DE 1968
Terremoto en Moyobamba. Kuroiwa y Deza (1968) describieron agrietamientos del suelo,
surgimiento de arena y agua por las grietas y grandes deslizamientos de tierra en la regin
epicentral. Los fenmenos de agrietamientos y surgimiento de agua fueron los ms numerosos,
especialmente a lo largo de las mrgenes del ro Mayo. Martnez Vargas (1969) present vistas
del afloramiento de arenas en forma de conitos de 10 a 20 cm. de dimetro producidos por el
fenmeno de licuacin en la terraza de Moyobamba. La mxima intensidad de este sismo fue de
VIII MM y su magnitud fue de Ms = 6.9.
107
UNIVERSIDAD NACIONAL DE INGENIERA
31 DE MAYO DE 1970
Terremoto que afect todo el departamento de Ancash y sur de La Libertad. Ericksen et al
(1970) y Plafker et al (1971) indicaron que en Casma, Puerto Casma y en zonas cercanas al litoral
en Chimbote, se produjo desplazamiento lateral del terreno causado por licuacin de depsitos
deltaicos y de playa, ocasionando grietas en el terreno que derrumbaron las estructuras que las
cruzaron. Las reas ms extensas de volcanes de arenas se formaron a lo largo del ro Casma,
entre Casma y Puerto Casma. Los volcanes tenan un crter central de unos cuantos centmetros
a 1 m. de dimetro, cercados por un montculo de arena y limo de hasta 15 m. de dimetro. Se
produjeron eyecciones de agua de un metro de altura. La zona central de Chimbote fue
evidentemente un rea de licuacin de suelos, as como de compactacin diferencial de la
cimentacin. El puente de Casma fue daado por licuacin de la cimentacin de los estribos. En
Chimbote y Casma y a lo largo de la Carretera Panamericana se notaron subsidencias
superficiales producto de la licuacin. La seccin residencial de Puerto Casma mostr evidencias
de asentamientos y eyeccin de agua. Cluff (1971) report fallas del terreno en Chimbote debido a
depsitos de playa saturados y sueltos. En Casma se produjo compactacin diferencial y
desplazamiento lateral del terreno debido a licuacin. Se produjeron inundaciones del terreno por
agua fretica, debido a la compactacin diferencial. En muchas reas se produjeron volcanes de
arenas y eyeccin de agua por existir nivel fretico alto. Berg y Husid (1973) indicaron evidencia
de licuacin de suelos en la cimentacin del Colegio Mundo Mejor, en Chimbote. Carrillo (1970)
indic descensos en los terraplenes de acceso de casi todos los puentes de la Carretera
Panamericana y asentamientos en las plataformas del Terminal Martimo de Chimbote. Tambin
se present evidencias del fenmeno de licuacin en los depsitos de arenas saturadas en la calle
Elas Aguirre en Chimbote y en el km. 380 de la Carretera Panamericana, cerca de Samanco.
Corporacin Hidrotcnica y C. Lotti (1979) reportaron licuacin generalizada en Puerto Casma,
produciendo agrietamientos de suelo y eyecciones de agua con arena. En Chimbote se produjeron
numerosos casos de licuacin y Puerto Casma se inund totalmente. Morimoto et al (1971)
describieron el fenmeno de licuacin de suelos en la ciudad de Chimbote. En la zona pantanosa
se produjo licuacin generalizada, con grietas debido a compactacin diferencial; y en la zona
aluvial licuacin subsuperficial con grietas y volcanes de arena. La mxima intensidad del sismo
fue de IX MM y su magnitud fue de Ms = 7.8.
9 DE DICIEMBRE DE 1970
Terremoto en el noroeste del Per. En el rea de Querecotillo en terraza fluvial y aluvial se
form un sistema de grietas en echeln, de longitud de 500 m. con aberturas de 0.30 m. y saltos
de 0.25 m. Se not efusin de arena formando sumideros de 0.60-1.00 m. de dimetro. Cerca al
casero La Huaca se agriet el suelo, brotando arena y lodo. En Tumbes cerca al Puerto Cura, en
las terrazas fluviales, se observ efusin de aguas negras acompaadas de arena que salieron a
la superficie a travs de grietas (Taype, 1971). La mxima intensidad de este sismo fue de IX MM
y su magnitud fue de Ms = 7.1.
20 DE MARZO DE 1972
Sismo en el nororiente. Segn Perales y Agramonte (1972), en el rea urbana de Juanju
se produjo el fenmeno de licuacin de suelos con sumideros alineados de hasta 1 m. de
dimetro. En la Carretera Marginal se produjeron asentamientos. Las aguas subterrneas variaron
su nivel esttico en ms de un metro. Se inspeccionaron dos pozos de agua que al momento de la
108
SECCIN DE POST-GRADO - FIC
visita se encontraban secos y taponeados con arena. La mxima intensidad de este sismo fue de
VIII MM y su magnitud fue de Ms = 6.9.
3 DE OCTUBRE DE 1974
Terremoto en Lima. Segn Huaco et al (1975) y Giesecke et al (1980), ocurrieron
fenmenos locales de licuacin en el valle de Caete, donde el nivel fretico es muy superficial. El
fenmeno local ms importante se encontraba en la Cooperativa La Quebrada, cubriendo un rea
de 30,000 m2. Maggiolo (1975) indic licuacin generalizada en Tambo de Mora, asociada a una
subsidencia o hundimiento, con densificacin posterior a lo largo de 4 km. paralelos a la lnea de
playa. En la zona norte se desarrollaron eyecciones de agua con arena a travs de volcanes de
arena. Espinosa et al (1977) indicaron posibles asentamientos diferenciales en El Callao debido a
licuacin de suelos y Moran et al (1975) presentaron vistas de posible licuacin en Ancn. La
mxima intensidad del sismo fue de VIII MM y su magnitud fue de Ms = 7.5.
CONCLUSIONES
Se considera que el mapa de reas de licuacin de suelos presentado puede estar incompleto y
no ser totalmente representativo de dicho fenmeno en el Per. Se espera que en el futuro se
realicen estudios adicionales que aporten nuevas evidencias sobre la ocurrencia de este
fenmeno, que servirn para modificar o completar el mapa propuesto.
AGRADECIMIENTO
Se agradece el apoyo brindado por el Centro Regional de Sismologa para Amrica del Sur
para la ejecucin de este trabajo, en especial a su director el Ing. Alberto Giesecke.
REFERENCIAS
1. Alva Hurtado, J. (1981), "Bibliografa sobre los Terremotos Peruanos: Aspectos Ingenieriles",
El Ingeniero Civil, No. 15, pp 50-51, Lima, Per.
2. Alva Hurtado, J. (1982), "Estudio del Potencial de Licuacin del Suelo de Cimentacin de las
Obras del Proyecto Central Hidroelctrica de Curumuy", Informe presentado a la Gerencia de
Estudios de Electroper.
3. Alva Hurtado, J. y Orihuela, P. (1982), "Estudio de Licuacin y Amplificacin Ssmica en
Chimbote, Per", IV Congreso Nacional de Ingeniera Civil, Chiclayo, Per.
4. Alva Hurtado, J. y Torres Cabrejos, R. (1983), "Mapa de Distribucin de Mximas
Intensidades Ssmicas en el Per", Revista TECNIA de la Universidad Nacional de Ingeniera,
Lima, Per.
109
UNIVERSIDAD NACIONAL DE INGENIERA
110
SECCIN DE POST-GRADO - FIC
20. Plafker, G., Ericksen , G.E. y Fernndez-Concha, J. (1971), "Geological Aspects of the May
31, 1970, Peru Earthquake", Bulletin of the Seismological Society of America, Vol. 61, No. 3,
pp. 543-578.
21. Silgado, E. (1946), "El Temblor Destructivo de Moyobamba del 6 de Agosto de 1945", Datos
Sismolgicos del Per, Boletn 3, Instituto Geolgico del Per, pp. 29-38.
22. Silgado, E. (1952), "El Sismo del 9 de Diciembre de 1950", Datos Sismolgicos del Per
1949-1950, Boletn No. 4, Instituto Nacional de Investigacin y Fomento Mineros, Ministerio
de Fomento y Obras Pblicas, pp. 47-51.
23. Silgado, E. (1957), "El Movimiento Ssmico del 12 de Diciembre de 1953", Anales del Primer
Congreso Nacional de Geologa, Parte II, Boletn de la Sociedad Geolgica del Per, Tomo
32, pp. 225-238.
24. Silgado, E. (1978), "Historia de los Sismos ms Notables Ocurridos en el Per (1513-1974)",
Instituto de Geologa y Minera, Boletn No. 3 Serie C, Geodinmica e Ingeniera Geolgica,
Lima, Per.
25. Silgado, E., Fernndez-Concha, J. y Ericksen, G.E. (1952), "El Terremoto del Cusco del 21 de
Mayo de 1950", Datos Sismolgicos del Per 1949-1950, Boletn No. 4, Instituto Nacional de
Investigacin y Fomento Mineros, Ministerio de Fomento y Obras Pblicas, pp. 27-46.
26. Taype, R.V., (1971), "Observaciones Geolgicas de los Daos Producidos por el Sismo del 9
de Diciembre de 1970 en los Departamentos de Piura y Tumbes", Servicio Nacional de
Geologa y Minera, Lima, Per.
27. Valencia, R. (1940), "El Terremoto del 24 de Mayo de 1940, sus Efectos y sus Enseanzas",
Revista de la Universidad Catlica del Per, Tomo VIII, pp. 294-329, Lima, Per.
111
UNIVERSIDAD NACIONAL DE INGENIERA
81 79 77 75 73 71 69
E C U A D O R C O L O M B I A
2
1953
PTO. PIZARRO
1970
BOCAPAN
1970
4
1970
1970
QUERECOTILLO
LA HUACA
1857
PIURA
1912
M OYOBAM BA
1945
6
1968
O
1972
J UANJUI
C
1746 B R A S I L
8 PATAZ
E
1970
CHIM BOTE 1970
CASMA
BARRANCA 1970
O
PTO.
CASMA
1970
10
1937
CHURUBAMBA
HUACHO 1966
1974
ANCON
12 LIMA 1974
1948
1974
CAETE
TAM BO CUSCO
P
1974 1950
DE MORA
1747
OROPESA
A
1716
PISCO
14 1690
AYAPATA
1664
C
1813
ICA 1746
R. VISECA
I
F
C
I
16
J.E. ALVA HURTADO (1983)
O AREQUIPA
1582
CAM ANA
1958
ISLAY
1582
AREA DE LICUACION 1868
LEYENDA :
Locu mba 1868
AREA DE PROBABLE LICUACION
18 SAMA
ARICA
ESCA L A : 1 : 5'000,000
100 80 60 40 20 0 100 100 km. CHILE
1868
112
SECCIN DE POST-GRADO - FIC
Captulo VII
MAPA DE AREAS DE DESLIZAMIENTOS
POR SISMOS EN EL PERU
(2)
Jorge E. Alva Hurtado
(1)
Luis A. Chang Chang
RESUMEN
INTRODUCCIN
El Grupo de Evaluacin del Peligro Ssmico del Proyecto SISRA (Sismicidad de la Regin
Andina), patrocinado por el Centro Regional de Sismologa para Amrica del Sur (CERESIS), se
fij como objetivo obtener del Proyecto un producto final que expresara de algn modo los niveles
relativos del peligro ssmico en el rea andina, y que fuera adems homogneo.
Se tom en consideracin el desarrollo alcanzado en cada pas del rea andina en esta
clase de estudios y el tiempo disponible para realizar el proyecto. El resultado final en comn con
los pases andinos del Grupo de Peligro Ssmico consisti en preparar mapas regionales a partir
de los antecedentes histricos existentes, graficando intensidades mximas, reas de licuacin de
suelos y zonas de deslizamientos asociados a terremotos.
Ponencia presentada en el V Congreso Nacional de Mecnica de Suelos e Ingeniera de Cimentaciones, Lima, 1315
Abril de 1987.
113
UNIVERSIDAD NACIONAL DE INGENIERA
DOCUMENTACIN BSICA
1513 A 1515
Silgado (1978), en base a la obra de Esquivel y Navia, reporta que en esos aos ocurrieron
grandes sismos, acompaados de formidables deslizamientos de tierras en el Cuzco.
22 DE ENERO DE 1582
Terremoto en Arequipa. Fray Victor Barriga (1915) comenta que "los cerros conmovidos
ayudaron al temor con el ruido, desgalgo de piedras, la tierra y el mucho polvo pareca que unos
con otros se apresuraban a sepultar a los vivientes".
28 DE FEBRERO DE 1600
31 DE MARZO DE 1650
Terremoto en el Cuzco. Silgado (1978) se refiere a la obra de Esquivel y Navia para indicar
que se produjeron grandes deslizamientos en las partes altas de Pisac y Paucartambo; uno de
ellos repres el curso del ro Apurimac. La tierra se agriet en varios lugares.
6 DE ENERO DE 1725
Movimiento ssmico en Trujillo. Silgado (1978) indica que en los nevados de la Cordillera
Blanca el sismo origin la rotura de una laguna glaciar, la cual desbordndose, arras un pueblo
cercano a Yungay, muriendo 1,500 personas.
114
SECCIN DE POST-GRADO - FIC
28 DE OCTUBRE DE 1746
Terremoto en Lima y tsunami en El Callao. Silgado (1978) indic que en algunos parajes
de Lucanas (Ayacucho) ocurrieron agrietamientos del terreno y deslizamientos debido a este
sismo.
4 DE MARZO DE 1904
Intenso movimiento ssmico en Lima. Silgado (1978) indica que en Matucana hubo
desprendimientos del material meteorizado de la parte alta de los cerros, y en Pasamayo fueron
profusos los deslizamientos en los acantilados de arena.
6 DE AGOSTO DE 1913
11 DE SETIEMBRE DE 1914
Se repiti en Caravel otro violento sismo que la volvi a dejar en escombros y ocasion
derrumbes en los cerros que rodeaban a la ciudad (Umlauff, 1915).
28 DE DICIEMBRE DE 1915
Segn Umlauff (1915), fortsimo temblor en Caravel donde cuarte viviendas, agriet parte de las
torres de la iglesia y caus deslizamientos en las partes sueltas de los cerros.
8 DE FEBRERO DE 1916
Sismo de origen coordillerano. Silgado (1978) indica que afect Ayacucho y Huancavelica.
En Ocros se cayeron viviendas y hubo deslizamientos en la parte alta de los cerros.
9 DE ABRIL DE 1928
Terremoto en Ayapata, Puno. Silgado (1978) refiere que en el pueblo de Ollachea se oan
detonaciones subterrneas y luego venan sacudimientos de tierra. Del nevado Allinocapae se
desprendieron con estrpido enormes masas de hielo. Los derrumbes de los taludes del ro
Esquilaya formaron una represa natural que ms tarde se rompi.
14 DE MAYO DE 1928
115
UNIVERSIDAD NACIONAL DE INGENIERA
24 DE DICIEMBRE DE 1937
Terremoto en las vertientes orientales de la Cordillera Central. Silgado (1978) anot que
hubo grandes deslizamientos en los cerros boscosos de Chontabamba y Chorobamba, en una
extensin aproximada de 50 km. En los caminos de herradura a Huancabamba se produjeron
enormes grietas y derrumbes de cerros.
2 DE JULIO DE 1938
24 DE AGOSTO DE 1942
Terremoto en Nazca. Silgado (1978) indica derrumbes en los cerros de los tramos de la
carretera de acceso al interior, especialmente a la altura de Puquio y en otros lugares de la
vertiente occidental, as como al noroeste de Cangallo, en la vertiente atlntica.
10 DE NOVIEMBRE DE 1946
Terremoto en Ancash. Heim (1949) destaca que grandes derrumbes se produjeron en las
quebradas de Pelagatos, Shuitococha, Huaychihuaco, Llama, San Miguel, Mayas, Mirador y
Quiches. Desprendidas de los cerros Novillo, Bandera y Acobamba. Silgado (1951) comenta que
en la quebrada Pelagatos se estim en unos 25 millones de metros cbicos la masa del material
grantico descompuesto que haba cado al valle. En el derrumbe de Huaychayaco, en la parte alta
de Quiches, se desprendieron de las faldas del cerro unos 5 millones de metros cbicos de
material calizo-margoso. Un deslizamiento de grandes proporciones sepult al casero de
Acobamba, donde murieron 217 personas. Otros derrumbes notables ocurrieron en Portachuelo,
cerca del nevado Pelagatos y sobre el ro Maran.
1 DE NOVIEMBRE DE 1947
Terremoto en Satipo. Silgado (1948) comenta que desde los vuelos de reconocimiento
entre San Ramn y Satipo se observ infinidad de derrumbes y desprendimientos de las partes
altas de las zonas boscosas hacia el fondo de las quebradas, los cuales causaron represamientos.
Deslizamientos en Carrizales, Calabazas y Cerro Pelado.
14 DE FEBRERO DE 1948
Violento temblor en Quiches. Silgado (1978) indica que los derrumbes de los cerros
ocasionaron la muerte de 7 personas, numerosos heridos y la interrupcin de varios tramos de los
caminos de herradura.
28 DE MAYO DE 1948
116
SECCIN DE POST-GRADO - FIC
21 DE MAYO DE 1950
Terremoto en el Cuzco. Silgado et al (1952) indican que en el lado sur del valle, al sur-este
del pueblo de San Sebastin, se observ en una longitud de 5 km. una zona de extensa
fisuracin. Las grietas eran de forma irregular, de abertura variable y de 10 a 50 m. de longitud.
Los deslizamientos fueron de magnitud reducida y ocurrieron principalmente en algunas laderas
escarpadas y en los taludes de gravas aluviales.
9 DE DICIEMBRE DE 1950
12 DE DICIEMBRE DE 1953
Un fuerte y prolongado movimiento ssmico afect la parte noroeste del Per y parte del
territorio ecuatoriano. Silgado (1957) indic deslizamientos de material suelto en El Alto, en los
alrededores de Zorritos y de las partes altas del Can del ro Tumbes.
21 DE ABRIL DE 1954
18 DE FEBRERO DE 1956
18 DE FEBRERO DE 1957
Movimiento ssmico sentido a lo largo de la costa, desde Huarmey hasta Chincha. Silgado
(1978) anota que en las cercanas del pueblo de Sayn, en el ro Huaura, los deslizamientos de
grandes bloques de piedra rompieron el muro de contencin de un canal de irrigacin. Derrumbes
de arena en los acantilados de Pasamayo.
15 DE ENERO DE 1958
Terremoto en Arequipa. Silgado (1978) indica que por efecto, tanto del Misti como de los
cerros vecinos, se desprendieron enormes piedras que al rodar cuesta abajo, levantaron gran
polvareda. Hacia la cordillera, los derrumbes daaron en varios tramos la lnea del ferrocarril a
Puno. La Carretera Panamericana entre Chala y Arequipa fue cubierta en varios trechos por
deslizamientos de magnitud variable, siendo la zona ms afectada entre Atico y Caman.
117
UNIVERSIDAD NACIONAL DE INGENIERA
3 DE MAYO DE 1962
17 DE OCTUBRE DE 1966
La ciudad de Lima fue estremecida por un sismo. Las vas de acceso al interior,
especialmente el km. 51 de la Carretera a Canta y el km. 22 de la Carretera Central quedaron
bloqueadas a consecuencia de los derrumbes. En la costa hubo deslizamientos de material suelto
de los acantilados de Chorrillos, Miraflores y Magdalena, que causaron gran polvareda. En
Pasamayo se produjeron desprendimientos de arena elica (Lomnitz et al, 1967).
19 DE JUNIO DE 1968
Terremoto en Moyobamba. Desde el avin, Kuroiwa y Deza (1968) pudieron divisar varios
deslizamientos en zonas de fuertes y moderadas pendientes en los cerros de Angaisha. Martnez
Vargas (1969) indic derrumbes en el rea afectada.
1 DE OCTUBRE DE 1969
31 DE MAYO DE 1970
5 DE MAYO DE 1971
118
SECCIN DE POST-GRADO - FIC
14 DE OCTUBRE DE 1971
20 DE MARZO DE 1972
3 DE OCTUBRE DE 1974
10 DE NOVIEMBRE DE 1980
17 DE ABRIL DE 1981
Sismo en Ayacucho. Ocola (1981) anota que en Opancca, Paccha y Ticllas se han podido
observar deslizamientos de rocas y derrumbes. Tambin en los cerros Moyo Orcco y Suyto Orcco
y en las quebradas de Molinos y Cruzhuayro.
CONCLUSIONES
119
UNIVERSIDAD NACIONAL DE INGENIERA
AGRADECIMIENTO
Se agradece el apoyo brindado por el Centro Regional de Sismologa para Amrica del Sur
(CERESIS) para la ejecucin del trabajo presentado, en especial a su Director el Ing. Alberto
Giesecke Matto.
REFERENCIAS
1. Alva Hurtado J.E. (1981), "Bibliografa sobre los Terremotos Peruanos: Aspectos Ingenieriles",
El Ingeniero Civil, Ao 3, No. 15, pp. 50-51.
2. Alva Hurtado J.E. (1983), "Breve Historia del Fenmeno de Licuacin de Suelos en el Per",
IV Congreso Nacional de Mecnica de Suelos e Ingeniera de Fundaciones, Lima-Per.
3. Alva Hurtado J.E., Meneses Loja J. y Guzmn Len V. (1984), "Distribucin de Mximas
Intensidades Ssmicas Observadas en el Per", V Congreso Nacional de Ingeniera Civil,
Tacna-Per.
4. Barriga V.M. (1951), "Los Terremotos en Arequipa 1582-1868", Biblioteca Arequipa, Tomo
Sptimo.
5. Descourt B. (1600), "El Terremoto de 1600 en Arequipa", Carta escrita en 1600, Coleccin
Zugger, Londres.
6. Deza E. (1971), "The Pariahuanca Earthquake, Huancayo, Peru. July-October 1969", Recent
Crustal Movements, Royal Society of New Zealand, Bulletin 9, pp. 77-83.
7. Ericksen G.E., Plafker G. y Fernndez-Concha J. (1970), "Preliminary Report on the
Geological Events Associated with the May 31, 1970 Peru Earthquake",U.S. Geological
Survey Circular 639.
8. Esquivel y Navia D. (1740), "Apuntes Histricos del Per y Noticias Cronolgicas del Cuzco
(Anales). Perodo 1543-1610", Cuzco.
9. Giesecke A., Ocola L. y Silgado E. (1980), "El Terremoto de Lima del 3 de Octubre de 1974",
Informe Preparado por el Centro Regional de Sismologa para Amrica del Sur a UNESCO
Contrato SC/RP 601-013.
10. Heim A. (1949), "Observaciones Geolgicas en la Regin del Terremoto de Ancash en
Noviembre de 1946", Vol. Jubilar de la Sociedad Geolgica del Per, Fascculo 6, Lima.
11. INGEMMET (1983), "Memoria Descriptiva del Mapa de Deslizamientos del Per", Proyecto
SISRA-CERESIS, Instituto Geolgico, Minero y Metalrgico, Lima-Per.
12. Kuroiwa J. y Deza E. (1968), "Daos Causados en Moyobamba por el Sismo del 19 de Junio
de 1968", Facultad de Ingeniera Civil, Universidad Nacional de Ingeniera, Lima-Per.
13. Lomnitz C., Esteva L., Flores R. y Cabr R. (1967), "Informe de la Misin de Reconocimiento
Sismolgico del Terremoto del 17 de Octubre de 1966", UNESCO, Pars.
14. Martnez Vargas A. (1969), "Apreciaciones Geolgicas de Campo sobre los Efectos del
Terremoto del 19 de Junio de 1968 en Moyobamba y Alrededores del Nor-Oriente Peruano",
Primer Congreso Nacional de Sismologa e Ingeniera Antissmica, Lima-Per.
120
SECCIN DE POST-GRADO - FIC
15. Martnez Vargas A. (1971), "Anlisis sobre los Huaycos, Aluviones y Alud-Aluviones ante el
Efecto Ssmico en el Per", III Congreso Nacional de Ingeniera Ssmica, Acapulco, Mexico.
16. Ocola L. (1981), "Actividad Ssmica en la Zona de Ayacucho y Programa de Estudios",
Direccin de Investigacin Cientfica de Geofsica Aplicada, Instituto Geofsico del Per, Lima-
Per.
17. Perales F., Perez G. y Taype V. (1972), "Apreciaciones Geolgicas en el Area Afectada por el
Sismo del 14 de Octubre de 1971", Servicio de Geologa y Minera, Ministerio de Energa y
Minas, Lima-Per.
18. Perales F. y Agramonte J. (1972), "Reconocimiento Geotcnico entre Juanju y Moyobamba
con Motivo del Sismo del 20 de Marzo de 1972", Servicio de Geologa y Minera, Ministerio de
Energa y Minas, Lima-Per.
19. Plafker G., Ericksen G.E. y Fernndez-Concha J. (1971), "Geological Aspects of the May 31,
1970, Peru Earthquake", Bulletin of the Seismological Society of America, Vol 61, No. 3, pp.
543-578.
20. Silgado E. (1948), "Datos Sismolgicos del Per-1947", Boletn No. 11, Instituto Geolgico del
Per, Lima-Per.
21. Silgado E. (1951), "The Ancash Earthquake of November 10, 1946", Bulletin of the
Seismological Society of America, Vol 41, No. 2, pp. 83-99.
22. Silgado E. (1952), "El Sismo del 9 de Diciembre de 1950", Datos Sismolgicos del Per 1949-
1950, Boletn No. 4, Instituto Nacional de Investigacin y Fomento Mineros, Ministerio de
Fomento y Obras Pblicas, Lima-Per, pp. 47-51.
23. Silgado E. (1957), "El Movimiento Ssmico del 12 de Diciembre de 1953", Anales del Primer
Congreso Nacional de Geologa, Parte II, Boletn de la Sociedad Geolgica del Per, Tomo
32, pp. 225-238.
24. Silgado E. (1957), "Datos Sismolgicos del Per 1952-1953", Boletn de la Sociedad
Geolgica del Per, Tomo XXIX, Lima-Per.
25. Silgado E, (1978), "Historia de los Sismos ms Notables Ocurridos en el Per (1513-1974)",
Instituto de Geologa y Minera, Boletn No. 3, Serie C, Geodinmica e Ingeniera Geolgica,
Lima-Per.
26. Silgado E. y Castro L. (1956), "Informe sobre los Sismos de Carhuaz", Instituto Nacional de
Investigacin y Fomento Mineros, Ministerio de Fomento y Obras Pblicas, Lima-Per.
27. Silgado E., Fernndez-Concha J. y Ericksen G.E. (1952), "El Terremoto del Cuzco del 21 de
Mayo de 1950", Datos Sismolgicos del Per 1949-1950, Boletn No. 4, Instituto Nacional de
Investigacin y Fomento Mineros, Ministerio de Fomento y Obras Pblicas, Lima-Per, pp. 27-
46.
123
28. Umlauff A.F. (1915), "La Regin Ssmica de Caravel", Boletn de la Sociedad Geogrfica de
Lima, Tomo XXXI, Segundo Trimestre, Lima.
121
UNIVERSIDAD NACIONAL DE INGENIERA
81 79 77 75 73 71 69
E C U A D O R C O L O M B I A
2
ANGAISHA 1968
CHAM AYA 1928
6
O
SAPOSOA 1972
C
B R A S I L
8
E
M AYAS
1946 PELAGATOS
1946
SIHUAS QUICHES
A
YUNGAY
1725,1970
N
CARHUAZ
1956
O
RECUAY 1970
10
HUANCABAMBA
1937
CHONTABAMBA
1937
PAUCARTAMBO
1937
1962
YUNGUL
SAYAN 1957 SATIPO 1947
TARMA 1938
PASAMAYO 1904,1957
PARIAHUANCA 1969
12 M IRAFLORES 1966,1974
CUSCO
HUAYTARA 1950 1942
CANGALLO
1513-1515,1950
A
OLLACHEA 1928
14 TORAYA 1971
C
1746
LUCANAS
PUQUIO 1942
I
F
CARAVELI 1 913-1914-1915
16
J.E. ALVA HURTADO O OCOA 1958
AREQUIPA 1582,1958
AREA DE DESLIZAMIENTO
LEYENDA :
AREA DE PROBABLE DESLIZAMIENTO
18
ESCA LA : 1 : 5'000,000
100 80 60 40 20 0 100 100 km. CHILE
122
SECCIN DE POST-GRADO - FIC
Captulo VIII
PELIGRO SSMICO EN EL PERU
INTRODUCCIN
El Per est comprendido entre una de las regiones de ms alta actividad ssmica que
existe en la tierra, por lo tanto est expuesto a este peligro, que trae consigo la prdida de vidas
humanas y prdidas materiales. Es necesario efectuar estudios que permitan conocer el
comportamiento ms probable de este fenmeno para poder planificar y mitigar los grandes
efectos que trae consigo. Una forma de conocer el probable comportamiento ssmico de un lugar
es mediante la evaluacin del peligro ssmico en trminos probabilsticos, es decir predecir las
posibles aceleraciones que podran ocurrir en un lugar determinado.
En las normas de diseo se especifican las cargas ssmicas, por lo que no es necesario
realizar investigaciones detalladas de la actividad ssmica del rea donde se construirn
estructuras comunes. El coeficiente de diseo ssmico a ser usado en el diseo ssmico pseudo-
esttico se determina en base a la zona, condicin del suelo e importancia de la estructura. Si la
estructura es flexible, la carga ssmica se modifica tomando en cuenta su periodo fundamental.
Sin embargo, cuando se planifican estructuras importantes, deben evaluarse sus capacidades de
resistir terremotos en base a estudios detallados de peligro ssmico. Tales estructuras incluyen:
grandes presas, puentes con luces grandes, tneles y centrales nucleares. Tambin se necesitan
estudios detallados para la evaluacin del peligro ssmico en una zona grande por urbanizar.
Ponencia presentada en el VII Congreso Nacional de Mecnica de Suelos e Ingeniera de Cimentaciones, Lima, 6-10
Diciembre 1993.
123
UNIVERSIDAD NACIONAL DE INGENIERA
SISMOTECTNICA
SISMICIDAD
Sismicidad Histrica
124
SECCIN DE POST-GRADO - FIC
Sismicidad Instrumental
Dentro de la metodologa para el clculo del peligro ssmico se considera que los eventos
ssmicos presentan una distribucin de Poisson, que se caracteriza por suponer independencia
entre los tiempos de ocurrencia, ya que cada uno de los sismos se considera como un evento
aislado e independiente. Por ello es necesario depurar del catlogo todas las rplicas y
premonitores, quedando los sismos como eventos principales.
mb = 3.30 + 0.40 Ms
En la Figura 1 se presentan todos los hipocentros del Catlogo Ssmico SISRA (1963-
1990) y los rasgos neotectnicos indicados por Machar et al (1991). Se observa que la actividad
ssmica en la zona Norte y Centro del pas est distribuida en dos fajas ssmicas longitudinales a
los Andes; una occidental a los Andes y exclusivamente producto de la subduccin con
hipocentros mayormente superficiales y algunos intermedios; y la otra, oriental a los Andes que
involucra tanto a procesos de subduccin (para hipocentros de profundidades intermedias, hasta
300 Km), como tambin a procesos secundarios, tal como la accin compresiva del escudo
brasilero contra el cinturn andino. Estas dos fajas ssmicas se unen en la zona de transicin
sismotectnica (13o-14o Sur), para constituir una sola amplia faja ssmica en la regin
sismotectnica del Per (Deza, 1990).
125
UNIVERSIDAD NACIONAL DE INGENIERA
Existe una actividad ssmica superficial causada por el proceso de reajuste tectnico del
Aparato Andino. En la Figura 1 se observa agrupamientos importantes de eventos en algunas
estructuras neotectnicas, tales como las fallas de Huaytapallana, fallas ubicadas en la sierra
central y en Moyobamba, en donde la actividad ssmica se encuentra en los primeros 40 Km de
profundidad. Los sismos recientes e histricos de Ayacucho, Cusco, Urcos y norte del lago
Titicaca, son manifestaciones de esta zona ssmica, muy superficial y destructiva (Ocola, 1989).
Cada una de las fuentes sismognicas tiene caractersticas propias definidas por sus
parmetros sismolgicos: magnitud mnima de homogeneidad (Mmin), pendiente de la distribucin
Gutenberg-Richter (b), tasa media anual de actividad ssmica (m) y magnitud mxima (Mmax). Las
escalas de magnitud ms utilizadas son mb y Ms. Dependiendo de la escala utilizada, los sismos
muestran valores asintticos a partir de una cierta magnitud (Idriss, 1985). Para evitar este
problema de saturacin de la magnitud se utilizar la magnitud M definida como max{mb, Ms}.
126
SECCIN DE POST-GRADO - FIC
Log N = a - b M
donde:
N = 10 a e- M
donde:
= b ln 10
LEYES DE ATENUACIN
127
UNIVERSIDAD NACIONAL DE INGENIERA
atenuacin de aceleraciones utilizada es la propuesta por McGuire (1974) para la Costa Oeste de
los Estados Unidos y ha sido empleada para las fuentes asociadas a sismos continentales.
n
P(A) = P (A/ E i) * P (E i)
i
Esta es la expresin que resume la teora desarrollada por Cornell en 1968, para analizar
el peligro ssmico. La evaluacin de esta integral es efectuada por el programa de cmputo RISK
desarrollado por McGuire (1976) en el clculo del peligro ssmico.
NIVEL DE CONFIDENCIA
donde:
t : tiempo de vida til
128
SECCIN DE POST-GRADO - FIC
Los movimientos de diseo que el ingeniero debe seleccionar estn asociados a un nivel
de excedencia suficientemente pequeo durante la vida til de la edificacin. En la Tabla 2 se
muestran valores representativos de criterios empleados en la seleccin de movimientos ssmicos
de diseo (Grases, 1989). La seleccin de los movimientos ssmicos dependen del tipo de obra.
Calculados los parmetros sismolgicos de las fuentes (Mmin, Mmax, , ), las profundidades
representativas de los hipocentros de las fuentes y seleccionadas las leyes de atenuacin, se
calcularon las aceleraciones horizontales mediante el programa RISK en una malla de puntos
(malla de 50x50 Km aproximadamente) en todo el territorio peruano y reas vecinas. En las
Figuras 4 y 5 se muestran los mapas de isoaceleraciones con un 10% de excedencia para 50 y
100 aos de vida til.
129
UNIVERSIDAD NACIONAL DE INGENIERA
Los resultados que muestran las Figuras 4 y 5 tienen una buena correlacin con el mapa
de Mximas Intensidades Ssmicas Observadas (Alva et al, 1984), en el cual se observa que las
zonas de Tumbes, Piura, Lima, Arequipa, Tacna y el Norte de Chile tienen intensidades entre VIII
y IX Mercalli Modificada y las intensidades ms bajas se presentan en la zona oriental con valores
por debajo de V MM.
CONCLUSIONES Y RECOMENDACIONES
Considerando que el territorio peruano se halla ubicado en una de las regiones de ms alto
ndice de actividad ssmica de la tierra, ha sido necesario evaluar apropiadamente el peligro
ssmico existente, prediciendo probabilsticamente las aceleraciones mximas que podran ocurrir
en cualquier punto del pas, utilizando leyes de atenuacin de aceleraciones y correlacionando la
sismicidad y la tectnica para determinar las fuentes sismognicas y sus respectivos parmetros
sismolgicos.
130
SECCIN DE POST-GRADO - FIC
Los resultados obtenidos en el presente estudio pueden ser utilizados para fines de
regionalizacin ssmica y otros estudios tales como, anlisis de vulnerabilidad, riesgo ssmico,
efectos de amplificacin y obtencin del espectro corregido de diseo, etc. Este estudio no es un
trabajo final, pues existen parmetros que cambiarn a medida que avancen las investigaciones,
produciendo mejores estimaciones del peligro ssmico en el Per.
REFERENCIAS
1. Aiquel A. (1990), "Hacia una Nueva Regionalizacin y Clculo del Peligro Ssmico en Chile",
Tesis de Grado, Universidad de Chile, Santiago, Chile.
2. Alva J., Meneses J. y Guzmn V. (1984), "Distribucin de Mximas Intensidades Ssmicas
Observadas en el Per", V Congreso Nacional de Ingeniera Civil, Tacna, Per.
3. Bender B. y Campbell K. (1989), "A Note on the Selection of Minimun Magnitude for use in
Seismic Hazard Analysis", Bulletin of the Seismological Society of America, Vol 79, No. 1,
pags 199-204.
4. Berrocal J. (1974), "South American Seismotectonics from SAAS Data", Thesis submitted for
the Degree of Doctor of Philosophy in the University of Edinburg.
5. Berrocal J., Deza E. y Shikiya J. (1975), "Estudio de Sismicidad para el Proyecto de
Derivacin del Ro Mantaro a Lima", Informe del Instituto Geofsico del Per a
ELECTROPERU S.A.
6. Bolt B. (981), "Terremotos", University of California, Berkeley.
7. Bonilla F. y Ruiz M. (1992), "Evaluacin del Peligro Ssmico en el Ecuador", Tesis de Grado,
Escuela Politcnica Nacional, Quito, Ecuador.
8. Casaverde L. y Vargas J. (1980), "Zonificacin Ssmica del Per", II Seminario
Latinoamericano de Ingeniera Sismo-Resistente, Organizacin de Estados Americanos y
Pontificia Universidad Catlica del Per, Lima, Per.
9. Cornell A. (1968), "Engineering Seismic Risk Analysis", Bulletin of the Seismological Society
of America", Vol 58, No. 5 pags. 1538-1606.
131
UNIVERSIDAD NACIONAL DE INGENIERA
10. Deza E. (1969), "Estudio Preliminar sobre las Zonas de Transicin que separan posibles
Regiones Sismotectnicas del Margen Occidental de Sudamrica: Zona de Transicin en el
Per", I Congreso Nacional de Sismologa e Ingeniera Antissmica, Lima, Per.
25. SISRA (1985), Catlogo de Terremotos para Amrica del Sur, Vol 7a, 7b y 7c, Proyecto
SISRA, CERESIS, Lima, Per.
132
SECCIN DE POST-GRADO - FIC
VIDA
PROBABILIDAD TIEMPO DE
TIPO DE OBRA UTIL
DE EXCEDENCIA RETORNO (aos)
(t aos)
133
UNIVERSIDAD NACIONAL DE INGENIERA
134
SECCIN DE POST-GRADO - FIC
135
UNIVERSIDAD NACIONAL DE INGENIERA
136
SECCIN DE POST-GRADO - FIC
137
UNIVERSIDAD NACIONAL DE INGENIERA
138
SECCIN DE POST-GRADO - FIC
Captulo IX
DISEO SSMICO DE PRESAS DE
TIERRA Y ENROCADO
Jorge E. Alva Hurtado (1)
Miguel Infantes Quijano (1)
RESUMEN
Se presenta una revisin de los mtodos existentes para realizar el diseo ssmico de presas
de tierra y enrocado. En primer lugar se presenta el mtodo pseudo-esttico que utiliza un
coeficiente lateral ssmico. Despus se ilustra el mtodo simplificado de deformacin permanente
y finalmente se describe el anlisis de estabilidad dinmico en base a la respuesta ssmica de la
presa. Los casos de ocurrencia de licuacin de suelos en materiales susceptibles no sern
tratados en este artculo. Con el objeto de ilustrar la metodologa se presenta su aplicacin a la
Presa Palo Redondo del Proyecto Chavimochic, diseo a nivel de factibilidad.
1. INTRODUCCION
Cualquier presa bien construida puede soportar sismos moderados, con aceleraciones
mximas de 0.2 g, sin daos. Las presas de material arcilloso con cimentacin arcillosa o rocosa
pueden soportar sismos fuertes con magnitudes Richter de 8.25 y aceleraciones mximas de 0.35
a 0.8 g, sin daos aparentes. Las presas de enrocado con pantalla de concreto se mantienen
secas y son capaces de soportar sismos extremadamente fuertes, con pequeas deformaciones.
El esfuerzo de diseo debe concentrarse en las presas que pueden presentar dao por sismo
severo o estn constituidas por materiales granulares saturados, que pueden perder su resistencia
o sufrir aumento de presin de poros. En este caso se requiere un anlisis dinmico, que
proporcione una base de diseo ms confiable.
139
UNIVERSIDAD NACIONAL DE INGENIERA
2. MTODO PSEUDO-ESTATICO
En este mtodo se emplea el clculo del factor de seguridad de la estabilidad de taludes por
equilibrio lmite, adicionando un coeficiente lateral ssmico. Uno de los principales problemas con
este mtodo es determinar el coeficiente lateral ssmico, que depende de la sismicidad del pas y
es de naturaleza semi-emprica (Seed y Martn,1966). Ruesta et al (1988) han propuesto valores
del coeficiente lateral ssmico para presas de tierra y enrocado en el Per, que se presentan en la
Figura 1. Los valores propuestos son consistentes con los valores utilizados en presas de tierra y
enrocado diseadas y construidas en el Per y en el mundo. Adems, el anlisis de estabilidad
pseudo-esttico sirve tambin para determinar el valor de la aceleracin de fluencia que se
emplea en el mtodo de deformacin permanente.
En este mtodo la deformacin permanente producida por un sismo se emplea como criterio
de diseo. Es un mtodo racional simple, que se aplica a presas formadas por suelos arcillosos
compactos, arenas secas y suelos granulares densos. En estos casos existe poco potencial de
desarrollo de presin de poros, se desarrollan deformaciones pequeas y el material retiene su
resistencia esttica. En este mtodo se evala la respuesta dinmica para calcular las
deformaciones permanentes. La falla ocurre en una superficie de deslizamiento bien definida con
comportamiento elstico hasta la falla y luego ocurre comportamiento perfectamente plstico.
Este mtodo fue propuesto originalmente por Newmark (1965). Sarma (1975) propuso usar
el modelo para analizar los efectos de las fuerzas de inercia y presin de poros en el factor de
seguridad, la aceleracin crtica y el desplazamiento. Makdisi y Seed (1977) desarrollaron un
mtodo simplificado basado en el concepto original de Newmark. La descripcin del mtodo es la
siguiente:
1. Se determina la aceleracin de fluencia, es decir el coeficiente lateral ssmico con el cual una
superficie potencial de falla desarrollar un factor de seguridad igual a la unidad. Los valores
de la aceleracin de fluencia estn en funcin de la geometra del terrapln, la resistencia del
material y la localizacin de la superficie potencial de deslizamiento.
2. Se determinan las aceleraciones producidas por el sismo en la presa por medio de un anlisis
de respuesta dinmica. Se emplean tcnicas de elementos finitos con propiedades del suelo
dependientes del nivel de deformacin o tcnicas unidimensionales ms simples. De estos
anlisis se determinan los tiempo-historia de aceleraciones promedio para las superficies
potenciales de falla.
140
SECCIN DE POST-GRADO - FIC
En este mtodo de anlisis se evala en primer lugar los esfuerzos estticos existentes en la
presa antes de la ocurrencia del sismo, tomando en consideracin la influencia de la secuencia
constructiva, el llenado de la presa y las condiciones de drenaje. Esta evaluacin se realiza
mediante el mtodo de elementos finitos con programas como el FEADAM 84 (Duncan et al,
1984) o el GEOSOFT (Parra, 1996) desarrollado en el CISMID de la UNI. Generalmente se
utilizan parmetros hiperblicos estticos para los materiales constituyentes de la presa (Duncan
et al, 1980).
Despus de calculados los esfuerzos estticos, se evala la respuesta ssmica por medio de
procedimientos analticos o mtodos numricos. El programa SHAKE (Schnabel et al, 1972) y las
versiones ms modernas de SHAKE91 (Idriss y Sun, 1992) y SHAKE21 (Kagawa, 1995) evalan
la respuesta mediante el anlisis unidimensional de propagacin de ondas con el uso para el
suelo de modelos lineales equivalentes (Seed e Idriss, 1970; Seed et al, 1984).
Para calcular las deformaciones permanentes producidas en una presa por la accin
ssmica, tanto por los mtodos simplificados como por el anlisis de estabilidad dinmico se
desarrollaron algoritmos y hojas de clculo, Infantes (1999).
141
UNIVERSIDAD NACIONAL DE INGENIERA
El Proyecto Chavimochic requiere de un futuro embalse para regular los caudales derivados
del Ro Santa. Este embalse se ubicar en la Quebrada Palo Redondo, antes que el canal
ingrese al tnel Intercuencas (Figura 5). El embalse tendr un volumen total de 370 millones de
metros cbicos, de los cuales 70 corresponden al volumen muerto de sedimentos. La Presa Palo
Redondo tendra una longitud de coronacin de 770 m y cerca de 480 m en la base, con una
altura mxima de 115 m sobre la cimentacin de roca. En el estudio de factibilidad de la Presa
Palo Redondo se han considerado tres alternativas de presa para el embalse:
El cuerpo principal de la presa estar formado por material de la Quebrada Palo Redondo,
del cual deber descartarse las partculas mayores de 12". En el talud aguas abajo, a modo de
relleno estabilizador, se colocar el mismo material, pero con bloques de hasta 24" de dimetro.
El material del cuerpo de presa se apoyar directamente sobre el material de quebrada previa
limpieza de los 0.5 m superiores y se colocar en capas de 1 m compactado con rodillo vibratorio.
El material aluvial de cimentacin tiene un espesor variable de 5 a 20 m. La pantalla de concreto
deber cimentarse en roca, la que requiere una cortina de inyecciones de impermeabilizacin y
consolidacin. Existen adems materiales de transicin, filtros y materiales de proteccin.
Resisten
Material Densidad cia
(KN/m3) C
(KPa) ()
Roca de Cimentacin 24 100 45
Material de Cimentacin 21 0 36
Cuerpo Principal de Presa 21 0 38
Bloque Estabilizador 22 0 42
Material de Transicin 20 0 36
142
SECCIN DE POST-GRADO - FIC
Se realiz el anlisis de estabilidad para la condicin al final de la construccin, por ser sta
la condicin ms desfavorable y estar la presa permanentemente seca. En el estudio original se
recomend el uso de un coeficiente lateral ssmico de 0.20 (Castillo y Alva Hurtado, 1993).
Para realizar los anlisis de estabilidad de taludes para las condiciones establecidas se utiliz
el programa de cmputo SLOPE/W. Los factores de seguridad calculados se detallan en la Tabla
2.
Si se comparan estos resultados con los factores mnimos establecidos por el US Corps of
Engineers para presas de tierra, se podr observar que los taludes de esta presa son adecuados.
Los parmetros utilizados en este anlisis, as como las curvas de reduccin del mdulo
cortante y factor de amortiguamiento se han estimado considerando que el material de la
cimentacin corresponde a un suelo gravoso de mal a bien graduado (Seed et al, 1984). El anlisis
de respuesta ssmica se realiz con la componente horizontal N 82 O del sismo de Lima-Per de
Octubre de 1974. Segn el Estudio de Peligro Ssmico de la zona (Castillo y Alva Hurtado, 1993),
para el sismo de diseo se ha determinado un valor de aceleracin mxima de 0.38 g con una
magnitud Ms igual a 7.5; el acelerograma del sismo anteriormente nombrado fue escalado a este
valor. Los parmetros dinmicos empleados se presentan en Infantes (1999).
Este procedimiento fue desarrollado por Makdisi y Seed (1977). Para realizar este anlisis
se tom en cuenta que el material del cuerpo de la presa es un material gravoso, para lo cual se
utilizaron los valores de los factores de reduccin del mdulo cortante y la razn de
amortiguamiento en funcin de la deformacin cortante efectiva correspondientes a este material.
Luego se determina una velocidad de corte inicial (Vs), que fue de 330 m/s. Con esto se
determina una relacin G/Gmax y una deformacin cortante inicial. Los resultados finales despus
de 5 iteraciones se presentan en la Tabla 3.
143
UNIVERSIDAD NACIONAL DE INGENIERA
Se observa que las aceleraciones mximas en la cresta son altas. Esto se debe a que el
material es gravoso. Este valor de amplificacin utilizando este procedimiento se ha obtenido en
este tipo de material en otras oportunidades.
A. ANALISIS UNIDIMENSIONAL
La seccin que se utiliz corresponde al eje de la presa. Su estratificacin parte desde una
cota de 230.00 m.s.n.m. donde se encuentra el basamento rocoso, hasta llegar a los 250.00
m.s.n.n.m, que corresponde al material gravoso de cimentacin. El material del cuerpo de la
presa alcanza hasta una cota de 345.00 m.s.n.m., que corresponde a la corona de la presa. Los
valores de los factores de reduccin del mdulo cortante y la razn de amortiguamiento en funcin
de la deformacin cortante efectiva correspondientes al material de cimentacin y al material del
cuerpo de la presa son los mismos que en el caso anterior. Se obtuvo una aceleracin en la
cresta de 0.47 g con el programa SHAKE21 y de 0.41 g con el programa SHAKE91.
148
B. ANALISIS BIDIMENSIONAL
En la Figura 7 se ilustra la malla de elementos finitos que se utiliz para los anlisis
bidimensionales, tanto estticos como dinmicos. Esta malla est compuesta por 769 nudos y
727 elementos. Para la evaluacin de los esfuerzos estticos se utiliz el modelo hiperblico. La
estimacin de los parmetros hiperblicos utilizados en el modelo constitutivo empleado se ha
realizado a partir de la informacin existente en la literatura para materiales similares analizados
(Duncan et al., 1980). Se utiliz el programa FEADAM84 (Duncan et al., 1984). En la Figura 8 se
muestra la distribucin de los esfuerzos cortantes obtenidos.
144
SECCIN DE POST-GRADO - FIC
similares. Por ejemplo, se obtuvo una aceleracin en la cresta de 0.66 g con QUAD4M y de 0.60
g con GEOSOFT.
Este anlisis se realiz con los mtodos de Newmark, Sarma y Makdisi y Seed, cuyos
resultados se presentan a continuacin.
A. MTODO DE NEWMARK
Se consider ocho superficies potenciales de falla, cuatro de ellas ubicadas aguas arriba y
las dems ubicadas aguas abajo. Son superficies que pasan de la corona a 1/4, 1/2 , 3/4 y una
altura total de la presa. Se calcul la aceleracin de fluencia, Ky, para cada una de estas
superficies potenciales de deslizamiento, es decir la aceleracin que desarrollara un factor de
seguridad igual a la unidad. Para esto se emple el programa SLOPE/W.
Tabla 4. Aceleracin de Fluencia (Ky) de las Superficies de Falla - Presa Palo Redondo
Sismo
Profundidad/Altura Total Lima 10/74
(y/H) N 82 0
Aguas Arriba Aguas Abajo
1/4 0.739 0.465
1/2 0.482 0.237
3/4 0.327 0.171
1 0.258 0.165
145
UNIVERSIDAD NACIONAL DE INGENIERA
Sismo
Lima 10/74
Profundidad/Altura Total
N 82 0
(y/H)
Aguas Arriba Aguas Abajo
(cm) (cm)
1/4 41.196 1.532
1/2 13.253 0.000
3/4 1.498 0.000
1 0.094 0.000
Las deformaciones obtenidas son muy pequeas, esto indica que la presa tiene una
capacidad muy grande de resistir sismos extremadamente fuertes.
B. MTODO DE SARMA
Para este clculo se consider un ngulo de friccin del material del cuerpo de la presa de
38 y un periodo predominante de vibracin de 1.20 seg. En la Tabla 7 se presenta el resumen de
resultados para cada uno de los taludes de la presa.
Deformacin
y/H Kc K Km Kc / Km Permanente
(cm)
0.294 1.35 0.513 0.573 80.002
0.266 1.00 0.380 0.700 36.376
0.249 0.70 0.266 0.936 4.057
1 0.215 0.50 0.190 1.132 0.084
150
B. Talud Aguas Abajo
Deformacin
y/H Kc K Km Kc / Km Permanente
(cm)
0.314 1.35 0.513 0.612 75.515
0.285 1.00 0.380 0.750 28.920
0.255 0.70 0.266 0.954 2.326
1 0.232 0.50 0.190 1.221 0.083
Para aplicar este mtodo se tom en cuenta los resultados obtenidos con el mtodo
simplificado de Makdisi y Seed para el clculo de la aceleracin mxima en la cresta y el periodo
natural, por ser sta la ms conservadora en trminos de la aceleracin mxima calculada. Las
aceleraciones de fluencia (Ky) se pueden obtener con algn procedimiento o programa de
146
SECCIN DE POST-GRADO - FIC
equilibrio lmite. Para este caso se utilizaron las mismas calculadas para el Mtodo de Newmark.
En la Tabla 8 se presenta las deformaciones permanentes para cada una de las profundidades de
la superficie potencial de deslizamiento para una magnitud Ms igual a 7.5.
Deformacin
Kmax Kv
y/H Ky Kmax Kmax.g.To Permanente
max Kmax
(cm)
1/4 0.29 0.85 0.959 0.302 0.140 97.089
1/2 0.27 0.60 0.677 0.399 0.070 34.267
3/4 0.25 0.44 0.496 0.504 0.038 13.641
1 0.23 0.35 0.395 0.583 0.022 6.282
U Deformacin
Kmax Kv
y/H Ky Kmax Kmax.g.To Permanente
max Kmax
(cm)
1/4 0.32 0.85 0.959 0.334 0.118 81.832
1/2 0.29 0.60 0.677 0.428 0.058 28.392
3/4 0.26 0.44 0.496 0.524 0.032 11.488
1 0.23 0.35 0.395 0.583 0.021 5.997
151
Segn Seed (1969), la deformacin permanente inducida por un sismo en una presa no
debe exceder de 3 pies, que es aproximadamente 1 m. Se observa que la presa cumple con este
requisito para las condiciones de diseo.
6. CONCLUSIONES Y RECOMENDACIONES
147
UNIVERSIDAD NACIONAL DE INGENIERA
* La zonificacin del coeficiente ssmico propuesto por Ruesta et al. (1988), no exime al
diseador de la realizacin de estudios de peligro ssmico para determinar dicho coeficiente en
el diseo ssmico definitivo de presas importantes.
* Para presas de enrocado los coeficientes ssmicos propuestos son menores que para presas
de tierra debido a que stas presentan un mayor ngulo de friccin interna y generalmente no
presentan problemas de presin de poros.
* Para las presas construidas de o sobre suelos granulares saturados, sueltos a medianamente
densos (presas de relave), sometidas a movimientos fuertes, una causa principal de dao o
falla es el incremento en la presin de poros en el suelo granular y la posible prdida de la
mayor parte de su resistencia, producto de este aumento de presin de poros (licuacin). No
es posible predecir este tipo de falla mediante el anlisis pseudo-esttico. En este caso se
requiere un procedimiento de anlisis dinmico, que proporcione una base ms confiable para
evaluar el comportamiento ssmico de este tipo de presa. Esto significa realizar un anlisis de
estabilidad post-sismo. La descripcin de este mtodo est fuera del alcance de este artculo.
* El anlisis riguroso requiere determinar los esfuerzos estticos iniciales antes que se produzca
el movimiento ssmico. Para este caso se concluye que los mtodos de elementos finitos son
adecuados para simular la secuencia constructiva y la respuesta dinmica posterior. El
148
SECCIN DE POST-GRADO - FIC
programa GEOSOFT que integra el anlisis esttico y dinmico es conveniente y fcil de usar,
por lo que se recomienda su empleo.
* Del anlisis del caso presentado se concluye que el procedimiento de evaluacin propuesto fue
el ms indicado pues se verific la validez de ste con mtodos ms sofisticados.
* Cabe indicar que la experiencia del proyectista juega un rol muy importante en la seleccin de
los procedimientos, programas, mtodos y parmetros ms adecuados para el diseo ssmico
de presas de tierra y enrocado.
REFERENCIAS
1) Alva Hurtado J.E. (1990), " Anlisis Dinmico de la Presa Palo Redondo-Nivel Factibilidad.
Alternativa de Materiales Gruesos con Pantalla de Concreto". Informe presentado a GMI
S.A., Proyecto Chavimochic.
153
2) Ambraseys N.N. y Sarma S.K. (1967), "The Response of Earth Dams to Strong
Earthquakes", Geotechnique, London, Vol. 17, pp. 180-213.
3) Bishop A.W. (1955), "The Use of the Slip Circle in the Stability Analysis of Earth Slopes",
Geotechnique, Vol. 5, pp. 7-17.
4) Castillo J. y Alva Hurtado J.E. (1993), "Peligro Ssmico en el Per", VII Congreso Nacional
de Mecnica de Suelos e Ingeniera de Cimentaciones", Lima, pp. 409-431.
5) Duncan J.M., Byrne P., Wong K.S. y Mabry P. (1980), "Strength Stress-Strain and Bulk
Modulus Parameters for Finite Element Analysis of Stresses and Movements in Soil
Masses", Geotechnical Engineering Research Report N UCB/GT 80-01, Department of Civil
Engineering, University of California, Berkeley.
6) Duncan J.M., Seed R.B., Wong K.S. y Ozawa Y. (1984), "FEADAM 84: A Computer
Program for Finite Element Analysis of Dams", Department of Civil Engineering, Virginia
Polytechnic Institute and State University, Blacksburg.
7) Geoslope International Ltd (1998), "Users Guide for Slope/W for Slope Stabillty Analysis".
8) Hudson M., Idriss I.M. y Beckae M. (1994), "QUAD4M: A Computer Program to Evaluate the
Seismic Response of Soil Structures using Finite Element Procedures and Incorporating a
Compliant Base", Department of Civil and Environmental Engineering, University of
California.
149
UNIVERSIDAD NACIONAL DE INGENIERA
9) Idriss I. M. y Sun J.I (1992), "SHAKE91 A Computer Program for Conducting Equivalent
Linear Seismic Response Analysis of Horizontally Layered Soil Deposits", National Institute
of Standards and Technology of Maryland and University of California, Davis.
10) Infantes Quijano M.A., (1999), "Anlisis Ssmico de Estructuras de Tierra", Tesis de Grado,
Facultad de Ingeniera Civil, Universidad Nacional de Ingeniera.
11) Kagawa T. (1995), "SHAKE 21: New Approach to Earthquake Response Analysis of
Horizontally Layered Soil Stress using the Equivalent Lineal Method", Wayne State
University, Detroit.
12) Makdisi F.I. y Seed H.B. (1977), "A Simplified Procedure for Estimating Earthquake-Induced
Deformations in Dams and Embankments", Report EERC. 77-19, University of California,
Berkeley.
16) Sarma S.K. (1975), "Seismic Stability of Earth Dams and Embankments", Geotechnique,
Vol. 25, No. 4, pp. 743-761.
17) Seed H.B. (1979), "Considerations in the Earthquake-Resistant Design of Earth and Rockfill
Dams", Geotechnique, Vol. 29, No 3 pp. 215-263.
18) Seed H.B. y Martin G.R. (1966), "The Seismic Coefficient in Earth Dam Design", Journal of
the Soil Mechanics and Foundation Division, ASCE, Vol. 92 SM3, pp. 25-58.
19) Seed H.B., Makdisi F.I. y De Alba P. (1978), "Performance of Earth Dams During
Earthquakes", Journal of the Geotechnical Engineering Division, ASCE, Vol. 104, GT7, pp.
967-994.
20) Seed H.B., Wong R.T., Idriss I.M. y Tokimatsu K. (1984), "Moduli and Damping Factors for
Dynamic Analysis of Cohesionless Soils", report EERC 84-14, University of California,
Berkeley.
21) Seed H.B. e Idriss I.M. (1970), "Soil Moduli and Damping Factors for Dynamic Response
Analysis", Report N EERC 70-10, University of California, Berkeley.
22) Schnabel P.B., Lysmer J. y Seed H.B. (1972), "SHAKE: A Computer Program for
Earthquake Response Analysis of Horizontally Layered Soils", Report N EERC 72-12,
University of California, Berkeley.
150
SECCIN DE POST-GRADO - FIC
81 79 77 75 73 71 69
0
ECUADOR COLOMBIA
2
III 4
BRASIL
II
10
OCANO
PACIFICO I III
12
14
B
COEFICIENTE SSMICO PROPUESTO
PARA PRESAS PEQUEAS Y O
(Ruesta, P.,Diaz, J. Y Alva, J., 1988)
MEDIANAS L
PRESAS PRESAS
16
ZONA DE TIERRA DE ENROCADO I
I 0.15 0.25 0.10 0.20
V
I
II 0.10 0.15 0.05 0.10
A 18
III 0.05 0.10 0.05
CHILE
151
UNIVERSIDAD NACIONAL DE INGENIERA
Viga de Corte
(rango de todos los datos)
0.4
y/h
0.6
Promedio de todos
0.8 los datos
1.0
0 0.2 0.4 0.6 0.8 1.0
Kmax / max
152
SECCIN DE POST-GRADO - FIC
10
M~8 1/4
M~7 1/2
U/Kmax gT0 - Segundos
0.1
0.01
M~6 1/2
0.001
0.0001
0 0.2 0.4 0.6 0.8 1.0
Ky / Kmax
153
UNIVERSIDAD NACIONAL DE INGENIERA
10
M~8 1/4
1
M~7 1/2
U/Kmax gT0 - Segundos
M~6 1/2
0.1
0.01
0.001
0.0001
0 0.2 0.4 0.6 0.8 1.0
Ky / Kmax
154
SECCIN DE POST-GRADO - FIC
COLOMBIA
ECUADOR
BRASIL
Localizacin del
La Libertad Proyecto
Ancash
OCEANO
PACIFICO
BOLIVIA
CHILE
155
UNIVERSIDAD NACIONAL DE INGENIERA
Material de Transicin
1.75 1 1.50
1 1 Bloque Estabilizador
Pantalla de Concreto 1
250.00 m.s.n.m.
(Terreno Natural)
156
SECCIN DE POST-GRADO - FIC
xy (KPa)
+3.10E+002
+2.50E+002
+2.00E+002
+1.00E+002
+0.00E+000
-5.00E+001
-1.50E+002
-2.00E+002
-3.00E+002
-3.30E+002
+2.20E+00
+2.00E+00
+1.80E+00
+1.60E+00
+1.40E+00
+1.20E+00
+8.00E+00
+4.00E+00
+0.00E+00
+6.00E-001
+5.00E-001
+4.50E-001
+4.00E-001
+3.50E-001
+3.25E-001
+3.00E-001
+2.75E-001
+2.50E-001
157
SECCIN DE POST-GRADO - FIC
Captulo X
ESTUDIO COMPARATIVO DE METODOLOGAS PARA EL
ANLISIS DE CIMENTACIN DE MAQUINARIAS
Rafael Salinas Basualdo (1)
Jorge E. Alva Hurtado (1)
RESUMEN
1. INTRODUCCION
158
UNIVERSIDAD NACIONAL DE INGENIERA
K = f (G, , B, L, R) . . (1)
Donde:
159
SECCIN DE POST-GRADO - FIC
Gazetas (1983, 1991), teniendo en cuenta los efectos que produce el empotramiento
(efectos de foso o excavacin y de contacto lateral), ha presentado una serie de expresiones con
las cuales son calculados las rigideces dinmicas y el amortiguamiento del sistema. Tales
expresiones son el resultado del anlisis de cimentaciones con diferentes condiciones y formas de
cimentacin, empleando modelos analticos relativamente ms refinados. Para considerar todos
los efectos, las rigideces tienen la forma:
Donde:
160
UNIVERSIDAD NACIONAL DE INGENIERA
El programa de cmputo tiene la ventaja de hacer posible los anlisis veloces a nivel de
diseo, con pequeas correcciones en los datos de entrada, ya sean cambios en las propiedades
de la cimentacin, por ejemplo, o los anlisis paramtricos en el caso que se quiera conocer
la influencia de la variacin de algunas caractersticas del suelo, por ser este factor el ms
importante y, frecuentemente, el de mayor dispersin. Asimismo, permite conocer la variacin de
los parmetros dinmicos y las amplitudes segn se realice el anlisis con rigideces estticas y
amortiguamientos independientes de la frecuencia, mediante la compilacin de Whitman, o con
rigideces dinmicas y amortiguamientos dependientes de la frecuencia de operacin, mediante la
compilacin de Gazetas. Tales evaluaciones ampliarn la capacidad de decisin del ingeniero al
momento de adoptar las dimensiones definitivas de diseo.
A partir del desarrollo de este trabajo, se presenta una secuencia del desarrollo de estudios
para la cimentacin de maquinarias vibratorias, que se resume en las siguientes etapas:
a. Conocimiento previo del conjunto de la infraestructura general del entorno de la cimentacin.
b. Conocimiento de las caractersticas de operacin, dimensiones y distribucin en planta de la
maquinaria. Determinacin de las fuerzas dinmicas, en caso de falta de informacin al
respecto.
c. Definicin de los criterios de diseo.
d. Definicin del perfil y de las propiedades dinmicas del suelo, en base a los ensayos de
campo.
e. Predimensionamiento de la cimentacin, segn recomendaciones generales o experiencia,
definicin del nivel del fondo de la cimentacin y condiciones de restriccin lateral. Debe
mencionarse los casos en que el fabricante de los equipos cuenta con planos con la forma del
bloque de cimentacin; en tal caso, los estudios deben tener fines verificatorios de las
dimensiones preestablecidas y sus condiciones de restriccin lateral.
f. Ejecucin del anlisis dinmico bajo las condiciones propuestas.
g. Comparacin de las respuestas mximas con las admisibles.
h. En caso de resultados no satisfactorios, modificacin de las dimensiones o de la condicin
lateral de la cimentacin, realizacin de un nuevo anlisis y comparacin de respuestas
mximas, hasta definir dimensiones de diseo.
161
SECCIN DE POST-GRADO - FIC
Adems, debe tenerse presente que las simplificaciones asumidas en los anlisis y la
heterogeneidad del suelo hacen altamente recomendable la medicin posterior, in-situ, de las
vibraciones presentes en la maquinaria en operacin, de tal forma que ser prevista la ocurrencia
de posibles daos, y se tomarn medidas finales para limitar las vibraciones excesivas.
4. CASOS ESTUDIADOS
El suelo est formado por arena fina con cierto contenido de finos no plsticos. El suelo de la
zona a la profundidad de cimentacin es de naturaleza compacta. El basamento rocoso se
encuentra aproximadamente a 30 m de profundidad; el nivel fretico no fue encontrado en las
excavaciones (CISMID, 1992). Se advirti la presencia de sales solubles. Debido a la constitucin
arenosa del suelo fue posible realizar el ensayo SPT. Sus resultados (Olcese y Zegarra, 1992)
dan un valor promedio de N igual a 70, indicando la gran rigidez del suelo a la profundidad de
desplante. Empleando expresiones empricas (CISMID, 1992) se tiene un rango de valores del
mdulo de corte de 35,900 a 22,100 t/m2. Adicionalmente, fueron realizados dos ensayos de placa
cclicos (CISMID, 1992); del anlisis de la parte cclica de la prueba, se defini para el mdulo de
corte un valor de 20,000 t/m2.
Datos de la Maquinaria:
Tipo Centrfuga
Peso de la mquina 331.37 t
Altura total del centro de gravedad a la base 4.22 m.
Solicitaciones:
162
UNIVERSIDAD NACIONAL DE INGENIERA
Los parmetros dinmicos del terreno y los resultados del anlisis se presentan en la Tabla 1.
El desplazamiento mximo considerando el movimiento no acoplado es del orden de 1.38 x 10-5
m, en tanto que para el movimiento acoplado resulta del orden de 1.60 x 10-5 m. Dado que el
desplazamiento lmite a la frecuencia de operacin es 7.6 x 10-5 m, de acuerdo a las
recomendaciones del fabricante, se concluye que la cimentacin cumple el requisito de diseo.
Se aprecia que, en este caso, aunque se produce un lgico incremento de las frecuencias
naturales de vibracin, en orden proporcional a la raz cuadrada del mdulo, las amplitudes
modales no crecen apreciablemente (Figura 5). Sin embargo, para tomar en cuenta la relacin
entre los desplazamientos y la frecuencia en que se producen, con el fin de evaluar la severidad
vibracional, se ha elaborado el grfico de desplazamientos espectrales presentado en la Figura 6.
El anlisis mediante la compilacin de Whitman presenta magnitudes mayores que aquellos
resultantes del mtodo de Gazetas, debido a que el modo rotacional crtico tiene un
amortiguamiento menor en el primer mtodo que en el segundo. Conociendo adems que se
presentarn acciones a frecuencias bajas en el estado de inicio o trmino de la operacin, es
posible que se presenten problemas de vibracin excesiva. Una alternativa para reducir los niveles
vibracionales es el aumento del ancho de la cimentacin, lo que ampliara la capacidad de
amortiguamiento y reducira las frecuencias naturales, alejndola de las frecuencias crticas167
(Salinas, 1995).
163
SECCIN DE POST-GRADO - FIC
de placa cclico. La forma del bloque de cimentacin se presenta en la Figura 7. Los datos
empleados en el anlisis fueron los siguientes:
Datos de la Maquinaria:
Tipo Reciprocante
Peso de la mquina 105.00 t
Altura del centro de gravedad a la base 1.80 m
Solicitaciones:
Dimensiones de la Cimentacin
La Tabla 2 presenta un resumen de los parmetros dinmicos y la respuesta vibratoria del sistema
analizado, que muestra un desplazamiento mximo de 1.42 x 10-6 m., superponiendo las
amplitudes de las dos componentes de frecuencia presentes. La amplitud mxima permisible es
del orden de 3.7 x 10-5 m y por lo tanto, se cumple satisfactoriamente el criterio de diseo
asumido.
Una serie de anlisis dinmicos considerando el rango de valores del mdulo de corte definido por
los ensayos, especficamente para 25000, 50000, 75000 y 100000 t/m2, y de acuerdo a los dos
mtodos implementados, han permitido tener una idea de la variacin de la respuesta con este
parmetro del suelo. Sus resultados, presentados en la Tabla 4 y grficamente en las Figuras 8 y
9, indican que el valor ms bajo del mdulo de corte conduce al resultado ms conservador.
Aunque el nivel de las fuerzas es pequeo, si stas aumentaran en diez veces su magnitud
nominal, el sistema an se encontrara en condiciones seguras, pues los desplazamientos
puntuales, cercanos en este caso a 1.5 x 10-6 m, seran menores al admisible.
164
UNIVERSIDAD NACIONAL DE INGENIERA
Dada la importancia de muchas de las obras de Ingeniera que cuentan con maquinarias
vibratorias, se recomienda considerar la uniformizacin de los criterios de diseo en un acpite de
la norma de diseo de cimentaciones, para el caso de acciones vibratorias. Se recomienda el
empleo de la metodologa presentada para el anlisis de la cimentacin de maquinarias en el
Per.
REFERENCIAS
165
SECCIN DE POST-GRADO - FIC
9. Moore P.J. (editor) (1985); Analysis and Design of Foundations for Vibrations. University of
Melbourne, A.A. Balkema Publishers; Rotterdam, Holanda.
10. Nieto J.A. y Resndiz D. (1967); Criterios de Diseo para Cimentaciones de Maquinaria.
Ingeniera, Vol.57, N3, pp.257-270; Mxico.
11. Olcese M. y Zegarra J. (1992); Estudio de Suelos con fines de Cimentacin Central Trmica
de Ventanilla. Informe Tcnico; Lima, Per.
12. Prakash S. (1981); Soil Dynamics. Mc.Graw-Hill Publishing Co.; New York, N.Y.
13. Richart F.E. (1962); Foundation Vibrations. Transactions ASCE, Vol.127(I), pp.863-898.
14. Richart F.E., Hall J. y Woods R. (1970); Vibrations of Soils and Foundations. Prentice-Hall
Inc.; Englewood Cliffs, N.J.
15. Salinas, R. (1995); Anlisis y Diseo de Cimentaciones Superficiales ante Acciones
Vibratorias. Tesis Profesional, Facultad de Ingeniera Civil, Universidad Nacional de Ingeniera;
Lima., Per.
16. Shrinivasulu P. y Vaidyanathan C.V. (1993); Handbook of Machine Foundations. Tata
Mc.Graw-Hill Publishing Co.; New Delhi, India.
17. Whitman R.V. (1972); Analysis of Soil-Structure Interaction. A State of the Art Review. MIT
Soils Publication N 300, Department of Civil Engineering, Massachussetts Institute of
Technology.
166
UNIVERSIDAD NACIONAL DE INGENIERA
Traslacin Rotacin
Modo de Vibracin Vertical
Transversal Longitudinal
Mtodo: Compilacin de Whitman
Parmetros Dinmicos
Frecuencia Natural (rpm) 777 864 1454
Inercia de Masas 142.38 142.38 1273.47
Rigidez 942300 1165000 29530000
Relacin de Amortiguamiento 0.83 1.15 0.38
Frecuencia de Resonancia (rpm) 1731
Parmetros Dinmicos
Frecuencia Natural (rpm) 1011 701 1191
Inercia de Masas 142.38 142.38 1273.47
Rigidez 1596000 767600 19820000
Relacin de Amortiguamiento 0.66 1.46 0.73
Frecuencia de Resonancia (rpm) 2704
167
SECCIN DE POST-GRADO - FIC
Rotacin
Modo de Vibracin Torsin
Transversal
Mtodo: Compilacin de Whitman
Parmetros Dinmicos
Frecuencia Natural (rpm) 1356 1685
Inercia de Masas 493.087 270.654
Rigidez 9947000 8429000
Relacin de Amortiguamiento 0.23 0.15
Frecuencia de Resonancia (rpm) 1281 1646
168
UNIVERSIDAD NACIONAL DE INGENIERA
TABLA 3. ANALISIS DE LA C.T. VENTANILLA, LIMA. RESULTADOS PRINCIPALES DE LA EVALUACION PARAMETRICA RESPECTO AL
MODULO DE CORTE DEL SUELO
TABLA 4. ANALISIS DE LA C.T. CALANA, TACNA. RESULTADOS PRINCIPALES DE LA EVALUACION PARAMETRICA RESPECTO AL
MODULO DE CORTE DEL SUELO
* El modo rotacional rx tiene acciones a la frecuencia de 360 rpm, para la cual se presentan los resultados en esta Tabla
169
SECCIN DE POST-GRADO - FIC
FZ
MZ
MX MY
FX FY
LZ2
LZ1
LX Y
X LY
0.05 L Pe
m Pr l
ite ec igr
, In
pa au o
c Es
ra i tru
0.02 M n
q Es ctu
Se
ui tr r
na u c al
ve
s tu
ra
0.01 ra
sp
ar
al
0.005
Amplitud de desplazamiento
as
Pr
pe
ob
rs
le
on
m
Cim
0.002 t
as
ic
as
+
en
F pa
tac
0.001 ci ra
lm
in
la
en s
t pe
de
0.0005 Po e rs
co pe on
M
rc as
q
N p e e p
o rc t ib +
uin
pe ep le
as
0.0002 rc tib s
ep le +
t ib s
le +
s
0.0001
100 200 500 2000 5000 10000
Frecuencia rpm
170
UNIVERSIDAD NACIONAL DE INGENIERA
INICIO
Clculo de Rigideces
Clculo de Rigideces y
Amortiguamientos
I= 1, N de frecuencias
Clculo de Amortiguamientos
N = 1,6
SI
Indicador = N
NO Clculo de Frecuencias
Naturales Segn el Procedimeinto
de Solucin Elegido
I = 1, N de frecuencias
Pueden existir SI
Modos acoplados
NO Se seleccion la SI
opcin de clculo
I = 1, N de frecuencias
Clculo de Amplitudes
lmites
N PUNTO DE
INTERES = 1, NP
NF = N de
frecuencias
N = 1,6
SI
indicador = N
NO Impresin de
Desplazamientos
Direcciones x, y, z
Impresin de Desplazamientos
Superpuestos
Impresion de Amplitudes
Lmites
FIN
171
SECCIN DE POST-GRADO - FIC
2.70
10.6
0
C.G.
GENERADOR
A A
18.40
13.2
3
C.G.
TURBINA
8.65
5.91
8.70
9 15
PLANTA
Z
C.G.
GENERADOR
Fz
TURBINA
Fy PTO. 1
1.77
PTO. 0.95
2
1.30 1.50
2.45 2.45
8.70
CORTE A - A
172
UNIVERSIDAD NACIONAL DE INGENIERA
2000
COMPILACION DE WHITMAN
1800
COMPILACION DE GAZETAS
FRECUENCIAS NATURALES (R.P.M.)
1600
x
1400 modo rotacional
eje longitudinal
x
1200 modo horizontal x
1000 x
x modo vertical
+
x
800
x
+
+
600 + +
+
+
+
400
200
0
8000 10000 12000 14000 16000 18000 20000
MODULO DE CORTE (Tn/m2)
173
SECCIN DE POST-GRADO - FIC
10
COMPILACION DE WHITMAN
COMPILACION DE GAZETAS
DESPLAZAMIENTOS MAXIMOS (x 0.1 mm)
G = 20000 Tn/m2
Vibraciones
severas
para las
personas
Vibraciones
problemticas
para las
personas
DESPLAZAMIENTOS HORIZONTALES
DESPLAZAMIENTOS VERTICALES
Vibraciones
fcilmente
apreciables
0.1
1000 10000
100
10
Vibraciones
COMPILACION DE WHITMAN
severas COMPILACION DE GAZETAS
para las
personas
G = 8000 Tn/m2
DESPLAZAMIENTOS MAXIMOS (x 0.1 mm)
Vibraciones
problemticas
para las
personas
DESPLAZAMIENTOS VERTICALES
1
DESPLAZAMIENTOS HORIZONTALES
Vibraciones
fcilmente
apreciables
0.1
1000 10000
1.00
174
UNIVERSIDAD NACIONAL DE INGENIERA
3.60
1.80 1.80
5.50
C.G. DE LA
MAQUINARIA
PTO.1 PTO.2
11.00
5.50
PTO.3
X
PLANTA
Mz C.G. EJE DE LA
MAQUINARIA
My
PTO.1
1.80
PTO.2
PTO.3
2.50
3.60
PERFIL
175
UNIVERSIDAD NACIONAL DE INGENIERA
COMPILACION DE WHITMAN
4500
FRECUENCIAS NATURALES (R.P.M.)
COMPILACION DE GAZETAS
4000
x
3500 x
modo torsional
x
3000 x
+
+
modo rotacional x
2500 eje transversal x +
+
2000 +
+
xx
1500 modo rotacional
+
+ eje longitudinal
1000
500
0
20000 30000 40000 50000 60000 70000 80000 90000 100000
COMPILACION DE WHITMAN
COMPILACION DE GAZETAS
2E-07
+
ROTACIONES MAXIMAS
+
1.5E-07
(rad)
ROTACION TORSIONAL
+ (720 rpm)
+
x
5E-08 x
+
+
+
+
x
x
x
x x
x
0
20000 30000 40000 50000 60000 70000 80000 90000 100000
176
SECCIN DE POST-GRADO - FIC
1.6E-06
COMPILACION DE WHITMAN
1.4E-06 x COMPILACION DE GAZETAS
DESPLAZAMIENTOS MAXIMOS (m)
x
1.2E-06
1E-06
DESPLAZAMIENTOS VERTICALES
8E-07
DESPLAZAMIENTOS HORIZONTALES
6E-07
x
x
4E-07
xx
x
x
2E-07
0
20000 30000 40000 50000 60000 70000 80000 90000 100000
177
SECCIN DE POST-GRADO - FIC
Captulo XI
MTODOS SIMPLIFICADO PARA LA EVALUACIN DEL
POTENCIAL DE LICUACIN DE SUELOS
Denys Parra Murrugarra (3)
Jorge E. Alva Hurtado (1)
RESUMEN
Este trabajo presenta una revisin de los mtodos disponibles para evaluar el potencial de
licuacin de un suelo, utilizando metodologas simplificadas en base a observaciones de campo
del comportamiento de depsitos de arena en terremotos pasados, es decir, mediante
correlaciones empricas entre el potencial a sufrir licuacin de un suelo y algunas caractersticas
fsicas y/o mecnicas del mismo, tales como: la resistencia a la penetracin, tipo de suelo, etc.,
adems de conocer la magnitud y aceleracin mxima superficial del sismo de diseo.
INTRODUCCIN
(1) Laboratorio Geotcnico del CISMID, Facultad de Ingeniera Civil, Universidad Nacional de Ingeniera, Lima, Per.
Ponencia presentada en el VI Congreso Nacional de Mecnica de Suelos e Ingeniera de Cimentaciones, Lima, del 11 al
13 de Noviembre de 1991.
178
UNIVERSIDAD NACIONAL DE INGENIERA
licuacin se ha debido a una serie de grandes terremotos que han ocurrido en todas las regiones
ssmicamente activas del mundo.
Existen bsicamente dos mtodos disponibles para evaluar el potencial de licuacin de una
arena saturada sometida a sacudimiento ssmico.
Debido al hecho que la mayora de los depsitos naturales de arena son muy poco
uniformes y a la extrema dificultad de obtener muestras inalteradas de suelos no cohesivos para
ensayos de laboratorio, se desarrollarn bsicamente los mtodos llamados empricos o
simplificados.
Por las razones indicadas, el uso de ensayos in-situ para evaluar la resistencia a la carga
cclica o caractersticas de licuacin de arenas, ha llegado a ser el mtodo preferido para la
evaluacin de estas propiedades en la prctica ingenieril en los pasados 15 aos. Puesto que no
existe ningn mtodo capaz de medir la resistencia a la licuacin de una arena in-situ, ha sido
necesario desarrollar correlaciones entre los valores de resistencia a la carga cclica determinados
de estudios de comportamiento de campo durante sacudimientos ssmicos y parmetros de
ensayos in-situ, los cuales sern probablemente indicativos de estas propiedades (Seed y De
Alba, 1986).
Debido a las dificultades para la realizacin de este ensayo, algunos parmetros de los
suelos pueden tal vez ser medidos con mayor exactitud sobre un amplio rango de profundidades,
179
SECCIN DE POST-GRADO - FIC
y en condiciones ambientales ms difciles. Sin embargo, debido 1a que el SPT ha sido tan
ampliamente usado en el pasado, la mayora de los datos de comportamiento de campo estn
generalmente correlacionados con este ndice de resistencia del suelo.
Seed e Idriss propusieron que la relacin de esfuerzo cclico para un suelo a una
profundidad determinada durante un terremoto puede ser evaluada mediante:
d = 0.65 a max o r
g o
d
o
N1 = Cn N
donde : Cn = factor de correccin que puede ser evaluado del grfico de la Fig. 1
180
UNIVERSIDAD NACIONAL DE INGENIERA
correlaciones estn dadas para arenas con diversos contenidos de finos y para un sismo de 7
de magnitud. La relacin mostrada puede ser fcilmente extendida a terremotos de diferentes
magnitudes, multiplicando la relacin de esfuerzos cclicos calculados por los factores de
correccin mostrados en la Tabla 1, la cual fue deducida a partir de estudios estadsticos basados
en el nmero de ciclos representativos que inducen terremotos de diferente magnitud, tomando
como base de comparacin una magnitud de 7 , segn se muestra en la Tabla 2 (Seed et al,
1975).
donde :
donde :
z = profundidad en metros
M = magnitud del terremoto.
Na = N1 + Nf
N1 = Cn N
181
SECCIN DE POST-GRADO - FIC
1
Cn = ; o en kg / cm2
o
donde :
La relacin de esfuerzo de corte (l / o) est dada por la Fig. 4. En este grfico se evala la
relacin de resistencia a la licuacin en funcin de los valores de N corregidos y para tres niveles
de deformacin por corte 2%, 5% y 10%.
En base al trabajo realizado por Seed e Idriss (1971), la carga dinmica inducida en el
elemento de suelo por un movimiento ssmico puede ser estimada mediante:
d = max = a max v r
v v g v d
donde :
Basados en un gran nmero de anlisis de repuesta ssmica del terreno, Iwasaki et al (1978)
propusieron la siguiente relacin para rd.
donde :
rd = 1-0.015 z
z = profundidad en metros
Por otro lado, Iwasaki et al (1978) condujeron numerosos ensayos triaxiales cclicos sobre
muestras de arena no disturbada. Como resultado concluyeron que la resistencia cclica no
drenada o relacin de esfuerzos cclicos que causar licuacin en un elemento de suelo sometido
a cargas dinmicas durante un terremoto puede ser evaluada mediante:
- para 0.02 D50 0.6 mm.
182
UNIVERSIDAD NACIONAL DE INGENIERA
l N
= 0.0882 + 0.225 log(0.35/ D50)
v v + 0.7
1 = 0.0882 N
- 0.05
v v + 0.7
Yasuda, Saito, Iwasaki y Morimoto, en 1985, propusieron un nuevo mtodo simplificado para
predecir el potencial de licuacin del ensayo de penetracin cnica con medicin de presin de
poros (CPTU). Este mtodo de prediccin es ms preciso que los mtodos convencionales,
debido a que no slo se toma en consideracin la resistencia por punta, qc, sino tambin la friccin
lateral, fS, y el exceso de presin de poros, u. Este mtodo es esencialmente vlido para predecir
el potencial de licuacin del terreno que soporta tuberas enterradas a travs de una gran rea.
183
SECCIN DE POST-GRADO - FIC
Cada uno de los mtodos simplificados descritos, basados en las experiencias de campo y
de laboratorio de sus autores, proponen diversas expresiones para evaluar, tanto la relacin de
esfuerzos que se requiere para causar licuacin en un suelo dado, cuyas caractersticas son
conocidas( l / v ),como la relacin de esfuerzos que induce un movimiento ssmico cuyas
caractersticas tambin son conocidas ( d / v ).
Por lo tanto, se puede definir el factor de seguridad contra la ocurrencia del fenmeno de licuacin
(FL), mediante la siguiente expresin.
( l / v )
FL =
( d / v)
Por lo tanto :
Si FL > 1 no se producir licuacin
Si FL 1 se producir licuacin
Se presenta en el Anexo I un ejemplo de la corrida del programa DLICUA para tres sondajes de
exploracin de campo, uno realizado en la ciudad de Chimbote y los otros dos en la ciudad de
Moyobamba, en las localidades de Tahuisco y Azungue, considerando diferentes sismos de
184
UNIVERSIDAD NACIONAL DE INGENIERA
diseo. El fenmeno de licuacin ocurri en Chimbote durante el terremoto de 1970 y en los dos
ltimos lugares durante los sismos de 1990 y 1991.
CONCLUSIONES Y RECOMENDACIONES
1. Licuacin es un fenmeno por el cual un suelo arenoso y saturado pierde toda su resistencia
al corte por el efecto de un sacudimiento ssmico severo, comportndose virtualmente como
un fluido viscoso.
2. Debido a que la mayora de los depsitos naturales de arena son poco uniformes y a la
extrema dificultad de obtener muestras inalteradas de suelos no cohesivos para realizar
ensayos de laboratorio, el uso de ensayos in-situ ha llegado a ser el mtodo ms aceptado
para la evaluacin de la resistencia a la carga cclica caracterstica de licuacin de arenas,
durante los pasados 15 aos.
4. Los mtodos simplificados implican el conocimiento de algunos parmetros del suelo, tales
como: la resistencia a la penetracin (valores de N qc) o la velocidad de ondas de corte,
entre otros. Sin embargo, debido a la gran cantidad de datos del ensayo SPT disponibles que
han servido para correlacionar el valor N con la resistencia a la licuacin del suelo, es que se
hace ms confiable utilizar este parmetro del suelo para evaluar su potencial de licuacin.
5. Se define como el factor de seguridad contra la ocurrencia del licuacin (FL), a la relacin
entre los esfuerzos inducidos por el terremoto de diseo y los esfuerzos requeridos para
causar licuacin en el suelo. No se producir el fenmeno si FL>1.
185
SECCIN DE POST-GRADO - FIC
REFERENCIAS
1. De Beer E., Goelen E., Heynen W. y Joustra K. (1988), Cone Penetration Test (CPT):
International Reference Test Procedure, Penetration Testing, ISSMFE, pp. 27-51.
2. Decourt L., Muromachi T., Nixon I.K., Schmertmann J. H., Thorbum S. y Zolkov E. (1988),
Standard Penetration Tests (SPT): International Reference Test Procedure, Penetration
Testing, ISSMFE, pp. 3-26.
3. Iwasaki T. (1986), Soil Liquefaction Studies in Japan. State of the Art, International Journal of
Soil Dynamics and Earthquake Engineering.
4. Iwasaki T., Tatsuoka F., Tokida K. y Yasuda S. (1978), A Practical Method for Assessing Soil
Liquefaction Potential Based on Case Studies at Various Sities in Japan, 5th Japan
Symposium on Earthquake Engineering, November.
5. Liao S. y Whitman R.V. (1986), Overburden Correction Factors for SPT in Sand, Journal of
Geotechnical Engineering, ASCE, March, pp. 373-377.
7. Robertson P.K. y Campanella R.G. (1983), Interpretation of Cone Penetration Test Part I
Sand, Canadian Geotechnical Journal, Vol 20, pp. 718-733.
8. Robertson P.K. y Campanella R.G. (1985), Liquefaction Potential of Sands Using the CPT,
Journal of the Geotechnical Engineering Division, ASCE, vol 111, N 3, pp. 384-403.
9. Seed H.B., Idriss I.M. y Arango I. (1983), Evaluation of Liquefaction Potential Using Field
Performance Data, Journal of Geotechnical Engineering, ASCE, vol. 109, N 3, pp. 458-481.
10. Seed H.B., Tokimatsu K., Harder L. y Chung R. (1985), Influence of SPT Procedures in Soil
Liquefaction Resistance Evaluations, Journal of Geotechnical Engineering, ASCE, vol. 111,
N 12, pp. 1425-1445.
11. Seed H.B. y De Alba P. (1986), Use of SPT and CPT Test for Evaluating the Liquefaction
Resistance of Sands, Use of Insitu Test in Geotechnical Engineering, ASCE, pp. 281-302.
12. Tokimatsu K., Uchida A. y Oh-oka H. (1989), Evaluation of Liquefaction Potential Based on
Shear Wave Velocity, 12th International Conference on Soil Mechanics and Foundation
Engineering, Rio de Janeiro, Brazil, pp. 151-156.
13. Tokimatsu K. y Yoshimi Y. (1981), Field Correlation of Soil Liquefaction with SPT and Grain
Size, International Conference on Recent Advances in Geotechnical Earthquake Engineering
and Soil Dynamics, University of Missouri-Rolla, pp. 203-208.
186
UNIVERSIDAD NACIONAL DE INGENIERA
14. Tokimatsu K. y Yoshimi Y. (1983), Empirical Correlation of Soil Liquefaction Based on SPT N-
value and Fines Content, Soils and Foundations, Vol. 23, N4, pp. 56-74.
15. Yasuda S., Saito K., Iwasaki K. y Morimoto I. (1989), Prediction of Liquefaction Potential
Using CPTU, 12th International Conference on Soil Mechanics and Fondation Engineering,
Rio de Janeiro, Brazil , pp. 131-134.
187
SECCIN DE POST-GRADO - FIC
8 0.89
7 1.00
6 1.13
6 1.32
5 1.50
8 26
7 15
6 10
6 56
5 23
188
UNIVERSIDAD NACIONAL DE INGENIERA
ANEXO I
(Ejemplo 1)
Cont. Esf.
Profundidad Clasif. Densidad D50 Valor N Esf. Efectivo
3 Finos Total 2 FRL 1 FRL 2 FRL 3
(mt.) (sucs) (Ton/m ) (mm) (gol/pie) (Kg/cm )
(%) (Kg/cm2)
1.45 SP-SM 1.80 5.0 0.215 5 0.261 0.261 0.577 0.675 0.848
2.55 SP-SM 1.80 5.7 0.211 11 0.459 0.364 0.839 0.921 0.915
3.65 SM 1.75 15.2 0.184 4 0.655 0.450 0.368 0.561 0.551
4.45 SP-SM 1.85 8.5 0.307 24 0.799 0.514 4.118 5.303 0.931
5.55 SM 1.85 16.2 0.165 21 1.003 0.608 2.298 5.086 0.941
6.45 SM 1.85 16.2 0.172 26 1.169 0.684 4.707 4.948 0.975
7.45 SM 1.85 21.1 0.197 21 1.354 0.769 1.671 4.919 0.831
8.45 SM 1.85 21.1 0.197 57 1.539 0.854 4.637 4.889 1.250
9.45 SM 1.85 27.4 0.145 63 1.724 0.939 4.657 4.882 1.339
189
SECCIN DE POST-GRADO - FIC
(Ejemplo 2)
SONDAJE : CPT 1
NIVEL FREATICO (mt) : 2.40
Cont. Esf.
Profundidad Clasif. Densidad D50 Valor N Esf. Efectivo
Finos Total FRL 1 FRL 2 FRL 3
(Ton/m3) 2
(mt.) (sucs) (mm) (gol/pie) (Kg/cm )
(%) (Kg/cm2)
1.00 SM 1.85 19.0 0.203 5 0.185 0.185 1.414 1.666 1.334
2.00 SC 1.80 37.0 0.174 7 0.366 0.366 1.741 2.337 1.517
3.00 SP-SM 1.90 6.0 0.271 2 0.555 0.495 0.174 0.642 0.650
4.00 SC 1.80 45.0 0.084 5 0.739 0.579 0.809 1.301 1.308
5.00 SC 1.80 45.0 0.084 3 0.919 0.659 0.437 1.031 1.049
6.00 SC 1.80 45.0 0.084 3 1.099 0.759 0.418 0.972 0.986
7.00 SM 1.85 26.0 0.200 11 1.281 0.821 1.128 1.529 1.045
8.00 SM 1.85 26.0 0.200 17 1.466 0.906 1.522 3.468 1.200
190
UNIVERSIDAD NACIONAL DE INGENIERA
(Ejemplo 3)
SONDAJE : CPT3
NIVEL FREATICO (mt) : 0.30
Cont. Esf.
Profundidad Clasif. Densidad D50 Valor N Esf. Efectivo
Finos Total FRL 1 FRL 2 FRL 3
(Ton/m3) 2
(mt.) (sucs) (mm) (gol/pie) 2 (Kg/cm )
(%) (Kg/cm )
1.00 SM 1.85 26.0 0.150 3 0.183 0.113 0.591 0.851 0.791
2.00 SP-SM 1.90 6.0 0.250 12 0.370 0.2900 1.264 2.217 0.989
3.00 SP-SM 1.90 6.0 0.250 8 0.560 0.290 0.782 0.795 0.769
4.00 SP-SM 1.85 6.0 0.250 10 0.749 0.379 0.820 0.862 0.811
5.00 SM 1.85 30.0 0.145 1 0.934 0.464 0.363 0.597 0.451
6.00 SM 1.85 30.0 0.145 3 1.119 0.549 0.361 0.678 0.601
7.00 SM 1.85 30.0 0.145 2 1.304 0.634 0.333 0.636 0.527
7.50 SP 1.95 .0 0.300 30 1.402 0.682 4.799 6.850 1.167
191
SECCIN DE POST-GRADO - FIC
2
ESFUERZO EFECTIVO VERTICAL Kg/cm
1
Dr= 40 - 60%
Dr= 60 - 80%
3
5
CN
0.6
35 15 5 Porcentaje
de Finos
0.5
o
l
0.4
RELACION DE ESFUERZO CICLICO
0.3
0.2
0.1
0
10 20 30 40 50
(N1)60
192
UNIVERSIDAD NACIONAL DE INGENIERA
12
0
0 10 20 30 40 50
0.5
= 10% 5% 2%
0.6
RELACION DE RESISTENCIA A LA LICUACIN
(RELACION DE ESFUERZOS DE CORTE l/)
0.4
0.3
= Deformacin
por corte
0.2
0.1
0
0 10 20 30 40 50
VALORES N CORREGIDOS
193
SECCIN DE POST-GRADO - FIC
1) Arya S.C., O'Neill M. y Pincus G. (1979), "Design of Structures and Foundations for
Vibrating Machines", Gulf Publishing Company, Houston..
2) Barkan D.D. (1962), "Dynamics of Bases and Foundations", Mc Graw Hill, New York.
10) Richart F.E., Woods R.D. y Hall J.R. (1970), "Vibrations of Soils and Foundations", Prentice-
Hall, New Jersey.
11) Whitman R.V. (1975), "Notas de Clase de Dinmica de Suelos", Instituto Tecnolgico de
Massachusetts, Cambridge.
194