Degtra
Degtra
Degtra
us ( z z , t ) = As e i ( t + k s *zs ) + Bs e i ( t k s *zs )
(1) (2)
ur ( zr , t ) = Ar e i ( t + kr *zr ) + Br e i ( t kr *zr )
donde T es la frecuencia angular de la onda armnica y k* es el nmero de onda complejo. La superficie del suelo est libre de esfuerzos de cizalla (zs = 0), de forma que:
An. VI - P7 - Pgina 1
donde Gs* = G ( 1 + 2i>) es el mdulo de cizalla complejo del suelo. Sustituyendo la ecuacin (1) en la ecuacin (3) y diferenciando se obtiene que:
Gs*ik s ( As e ik s ( 0) Bs e ik s ( 0) )e i t = Gs*ik s ( As Bs )e i t = 0
(4)
que se satisface cuando As = Bs . La compatibilidad de los desplazamientos y la continuidad de los esfuerzos en la frontera suelo/roca requiere: us (zs = H) = ur (zr = 0) Js (zs = H) = Jr (zr = 0) Sustituyendo las ecuaciones (1) y (2) en la ecuacin (5): (5) (6)
As (e ik s H + e ik s H ) = Ar + Br
De la ecuacin (6) y la definicin de esfuerzos de cizalla ( = Gs*
(7)
u ) se tiene: z
As iGs* k s* (e ik s H e ik s H ) = iGr* k r* ( Ar Br )
o
Gs* k s* * * ik s H e ik s H ) = Ar Br * * = As ( e Gr k r
La razn:
(8)
Gs* k s* * * * = z Gr k r
donde "z* se conoce como la razn de impedancia compleja. Resolviendo las ecuaciones (7) y (8) simultneamente, se obtiene que:
1 * * * * As (1 + z )e ik s H + (1 z )e ik s H 2 1 * * * * Br = As (1 z )e ik s H + (1 + z )e ik s H 2 Ar =
[ [
] ]
(9a) (9b)
Si una onda de cizalla de amplitud A, propagndose verticalmente, viajara hacia arriba a travs de la roca y el suelo no estuviera presente, el efecto de la superficie libre en
An. VI - P7 - Pgina 2
2 As =
4A (1 + )e
* z
* ik s H
* + (1 z )e ik s H
La funcin de transferencia, F(T) , definida como la razn de la amplitud de la superficie del suelo a la del sustrato rocoso, est dada por:
F ( ) = 2 (1 + )e
* z
* ik s H
* + (1 z )e ik s H
Obviamente, la funcin de transferencia es una funcin compleja que puede reescribirse usando la ley de Euler como:
F ( ) =
1 * cos( k H ) + i z sin( k s* H )
* s
(10)
2.2 Caso de mltiples capas de suelo La nica diferencia con respecto al caso anterior reside en que la funcin de transferencia no es la misma para un depsito de diferentes capas de suelo. La funcin de transferencia para un depsito de este tipo debe tener en cuenta la reflexin y transmisin de ondas en las interfases entre capas adyacentes. Consideremos el depsito de suelo mostrado en la Figura 2. Dentro de una lmina dada, por ejemplo la j, los desplazamientos horizontales vendrn dados Figura 2. Caso de mltiples capas de suelo. por:
u j (z j , t ) = ( Aj e
ik * z j j
+ Bje
ik * z j j
)e it
(11)
En la frontera entre la capa j y la capa j+1, la compatibilidad de los desplazamientos requiere que:
An. VI - P7 - Pgina 3
+ Bje
ik l*h j
(12)
* j +1
( Aj e
ik *h j j
+ Bje
ik *h j j
(13)
Observa como las ecuaciones (12) y (13) son anlogas a las ecuaciones (7) y (8), respectivamente. Definiendo "j * como la razn de impedancia compleja en la frontera entre las capas j y j+1, las amplitudes de onda para la capa j+1 pueden obtenerse a partir de las amplitudes de la capa j resolviendo las ecuaciones (12) y (13).
A j +1 =
1 1 ik *h ik *h A j (1 + * )e j j + B j (1 * )e j j j j 2 2
(14a) (14b)
B j +1 =
1 1 ik *h ik *h A j (1 * )e j j + B j (1 + * )e j j j j 2 2
En la superficie del suelo (z1 = 0), la condicin de que el esfuerzo de cizalla sea cero implica que A1 = B1 . Aplicando las ecuaciones (14) para j = 1, 2, 3,..., N, los coeficientes Aj+1 y Bj+1 pueden relacionarse con Aj y Bj por: Aj+1 = aj+1 (T) A1 (15a) Bj+1 = bj+1 (T) B1 (15b) donde las funciones aj+1 (T) y bj+1 (T) representan los efectos de las interacciones de las ondas que tienen lugar en todas las interfases de las capas por encima de la capa j+1. Entonces, una funcin de transferencia que relaciona los movimientos en los techos dos capas cualquieras, por ejemplo i y j, pueden expresarse como:
Fij ( ) = ai ( ) + bi ( ) a j ( ) + b j ( )
(16)
2.3 Anlisis lineal equivalente. El comportamiento no lineal y anelstico del suelo est bien establecido en la ingeniera geotcnica. La no linealidad del comportamiento esfuerzo-deformacin del suelo implica que el mdulo de cizalla del suelo est cambiando continuamente. La anelasticidad implica que el suelo se descarga siguiendo un camino diferente al camino de carga,
An. VI - P7 - Pgina 4
Para aproximar la respuesta real del suelo no lineal y anelstica, se puede utilizar una aproximacin lineal equivalente. En esta aproximacin, se llevan a cabo anlisis lineales con propiedades de suelo que se ajustan iterativamente para ser consistentes con un nivel efectivo de deformacin de cizalla inducido en el suelo. En la aproximacin lineal equivalente, el mdulo de cizalla secante, como se muestra en la Figura 3, aproxima un modulo de cizalla promedio durante el ciclo completo de carga. Conforme aumenta el nivel de la deformacin de cizalla, el modulo de cizalla secante decrece. La relacin entre el mdulo de cizalla secante y la amplitud del esfuerzo de cizalla puede caracterizarse por medio de la curva de reduccin de modulo. La naturaleza de esta curva, que tiene una ordenada de razn de modulo (= G/Gmax ) y una abscisa de log (deformacin de cizalla), ha sido muy bien establecida para muchos suelos. El algoritmo de solucin usado en SHAKE supone amortiguamiento viscoso del suelo que se representa usando un modulo complejo de cizalla. El amortiguamiento viscoso implica un comportamiento que estara caracterizado por bucles esfuerzo-deformacin elpticos. Puesto que los bucles esfuerzodeformacin reales son rara vez elpticos, se usa una razn de amortiguamiento equivalente. La razn de amortiguamiento equivalente es igual a la razn de amortiguamiento que se obtendra basada en el rea dentro del bucle de histresis, el mdulo de cizalla secante y la deformacin de cizalla mxima. La relacin entre esta razn de amortiguamiento equivalente y la deformacin de cizalla se caracteriza por medio de una curva de amortiguamiento. La naturaleza de esta curva, que tiene una ordenada de razn de amortiguamiento y una abscisa de log (deformacin de cizalla) ha sido bien establecida para muchos suelos. En un anlisis lineal equivalente, la primera iteracin se realiza usando el mdulo de cizalla y las razones de amortiguamiento que corresponden a niveles de deformacin de cizalla estimados inicialmente. En el software que utilizaremos, la primera iteracin est basada en una deformacin de cizalla asumida de 0.001%. Siguiendo la primera iteracin,
An. VI - P7 - Pgina 5
donde R( es un factor de reduccin de la deformacin que a menudo se toma como: R( = (M -1 ) / 10 El mdulo de cizalla y la razn de amortiguamiento correspondientes a (eff se usan, entonces, para una nueva iteracin. Este proceso se repite hasta que la deformacin efectiva calculada no cambia mucho de una iteracin a otra. En este momento, se dice que el proceso lineal equivalente ha convergido. Aunque la aproximacin lineal equivalente permite aproximar los efectos ms importantes de comportamiento no lineal y anelstico de los suelo, debemos llamar la atencin a que sigue siendo un mtodo lineal de anlisis. El modulo de cizalla compatible con la deformacin y la razn de amortiguamiento permanecen constantes durante toda la duracin de un terremoto - cuando las deformaciones inducidas en el suelo son pequeas y cuando son grandes. Las deformaciones permanentes no se pueden calcular y tampoco las presiones intersticiales. Sin embargo, la aproximacin lineal equivalente se ha mostrado como una metodologa que proporciona estimaciones razonables de la respuesta ssmica del suelo bajo muchas condiciones de importancia prctica. 3. METODOLOGA 3.1 Introduccin El primer paso en esta experiencia ser iniciar el programa EduShake. Haz click sobre el cono EduShake para ejecutar el programa. La primera pantalla muestra el smbolo de Edushake y una barra de herramientas con seis botones, que de izquierda a derecha seran: Input Manager (Entrada de datos), Solution Manager (Ejecucin del problema), Output Manager (salida de datos), Utilities (Utilidades), Help (Ayuda), y Exit (Salir del programa) como se muestra en la Figura 4. EduShake incluye una utilidad que permite aadir nuevas curvas de reduccin de mdulo y de amortiguamiento o modificar las incluidas en las libreras del programa. EduShake posee adems un sistema completo de ayuda que opera como una ayuda tpica de Windows y a la que puede accederse de dos formas diferentes. Primero, haciendo click en el botn Help de la barra de herramientas y la segunda pulsando la tecla F1 una
An. VI - P7 - Pgina 6
3.2 Inicio de un anlisis Para llevar a cabo un anlisis en EduShake es necesario definir las propiedades de todas las capas de suelo dentro del perfil que se va a analizar, especificando las caractersticas del movimiento de entrada que ser aplicado al perfil de suelo, definiendo los parmetros que se van a calcular y documentando los datos de entrada. Sigue los pasos que se describen a continuacin para llevar a cabo el anlisis propuesto en esta experiencia. Definicin de las propiedades de las capas 1. Haz click sobre el botn Input Managerde la barra de herramientas. Esto abrir una nueva ventana que te permitir definir un perfil de suelo (Profile), seleccionar un movimiento de entrada (Input Motion) y guardar un registro de tus datos de entrada en un documento (Figura 5). Tambin es posible abrir un fichero de entrada existente sin ms que pinchar en el botn File y Open de la barra de herramientas de esta ventana.
An. VI - P7 - Pgina 7
2.
Escribe un ttulo para el perfil del suelo en la casilla Profile Title. Este ttulo debe identificar tu perfil de forma que puedas distinguirlo del de otros introducidos por tus compaeros. Puedes escribir Perfil del Alumno ..., por ejemplo. Pulsa la tecla Tab de tu teclado para pasar automticamente al siguiente campo que has de rellenar. El curso se situar en el campo No. of Layers. Introduce el nmero de capas para el anlisis. Para esta experiencia, introduce 16' y pulsa la tecla Tab (as consideramos un perfil de suelo compuesto por 15 capas ms el semiespacio infinito sobre el que yacen). Observa como EduShake proporciona 16 campos cuando se introduce el nmero de capas - cada campo te permitir introducir los datos apropiados para cada capa (Figura 6). Introduce la profundidad del nivel fretico. Puesto que EduShake lleva a cabo un anlisis de esfuerzos totales, la profundidad del nivel fretico no es necesaria por un anlisis de respuesta del suelo. Para esta experiencia, introduce el valor 0.0 en el campo correspondiente a pies (ft) o metros (m). EduShake permite introducir los datos de entrada tanto en unidades del SI como en US. A continuacin hemos de introducir los datos correspondientes a cada una de las capas. Asegrate de estar en el campo correspondiente a la Capa 1 (la etiqueta 1 estar en color azul). Introduce el nombre del material correspondiente a esta capa dentro del campo Material Name. Este puede ser cualquier cadena alfanumrica, por ejemplo, Arcillas limosas poco consolidadas.
An. VI - P7 - Pgina 8
3.
4.
5.
Figura 6. Ventana para la caracterizacin de cada una de las capas del perfil de suelo de anlisis.
6.
El siguiente paso ser la seleccin de una curva de degradacin de los mdulos G/Gmax. EduShake tiene implementada un conjunto de estas curvas extradas de la literatura de la ingeniera geotcnica. Para ver la lista de estas curvas, haz click sobre el botn al lado derecho del campo correspondiente a Modulus Reduction Curve. En esta experiencia seleccionaremos la curva Vucetic-Dobry (Figura 7). El modelo Vucetic-Dobry describe el comportamiento de degradacin de los mdulos (y del amortiguamiento) como una funcin del ndice de plasticidad (introduce un valor 10 para el ndice de plasticidad). EduShake interpolar entre la curvas de degradacin de los mdulo (y amortiguamiento) presentadas por Vucetic y Dobry para obtener curvas que correspondan al ndice de plasticidad que se ha introducido. Observa que EduShake inicialmente ajusta el modelo de amortiguamiento (Damping Curve) exactamente igual que el modelo de degradacin de los mdulos (Figura 7). Si deseramos cambiarlo lo haramos de forma similar a la explicada en el paso 6 para la curva de degradacin de los mdulos. En nuestro caso mantendremos los valores por defecto. Comprueba las curvas de degradacin de los mdulos y amortiguamiento pinchando en el botn Plot Modulus and Damping Curves. Esto mostrar la curva seleccionada en negro y las curvas del modelo Vucetic-Dobry en verde (Figura 8). Pulsa Return para volver a la pantalla principal.
An. VI - P7 - Pgina 9
7.
8.
9.
Ahora introduce en el campo correspondientes al espesor (Thickness), bien en pies (ft) o en metros (m), un valor de 1.5 m para la Capa 1 (5 ft).
Figura 8. Ventana correspondiente a las curvas de modulo de cizalla y amortiguamiento seleccionadas. An. VI - P7 - Pgina 10
12.
13.
14. 15.
An. VI - P7 - Pgina 12
16.
Ahora debemos especificar que informacin se va a calcular durante el anlisis de respuesta del suelo. Ve al campo correspondiente a la Capa 1 y haz click en el botn Select Output (seleccionar salida). Sobre el formulario que aparece, seala las cajas correspondientes a velocidad y aceleracin. Tambin marca la caja correspondiente al espectro de respuesta en aceleraciones e introduce las razones de amortiguamiento de 5, 10 y 20 % (Figura 13). Repite este proceso para las Capas 2 y 11, tambin marcando los grficos de deformacin de cizalla y esfuerzo de cizalla para dichas capas. Ahora que todas las capas estn definidas, revisa que no existan errores haciendo click sobre el botn View Profile (ver perfil). El perfil muestra la variacin del peso por unidad de volumen y la velocidad de ondas de cizalla con la profundidadlos valores escogidos para esta experiencia muestra un gran contraste de impedancias (producto de la densidad y la velocidad de ondas) en la frontera entre las arcillas limosas y las arcillas consolidadas (fondo de la Capa 10) y entre las arcillas consolidadas y el substrato rocoso (fondo de la Capa 15). Como veremos los contrastes de impedancias tendrn una gran influencia sobre la respuesta ssmica de nuestro perfil de suelo. Las localizaciones para la cuales se va a calcular la salida aparecen indicadas por ovalos verdes(Figura 14).
17.
An. VI - P7 - Pgina 13
An. VI - P7 - Pgina 14
Figura 14. Resumen del perfil de suelo a analizar mediante el programa EduShake.
Especificacin del movimiento del suelo. 18. El siguiente paso consiste en la especificacin del movimiento del suelo. Haz click sobre el botn Input Motion (Movimiento de entrada). Aparecer un formulario que permitir seleccionar un movimiento del suelo, definir sus caractersticas, verlo grficamente y calcular varios parmetros del movimiento del suelo asociados con el movimiento (Figura 15). Introduce 1' para el nmero de movimientos. EduShake permite usar slo dos movimientos del suelo pero se puede analizar un perfil con mltiples movimientos de entrada cambiando varias caractersticas (e.g. la aceleracin pico) o los registros disponibles de movimiento del suelo. Para esta experiencia, deja la razn de deformacin, el mximo nmero de iteraciones y el error de tolerancia en sus valores por defecto. Haz click sobre el botn Open en el campo correspondiente a Motion 1'. Del menu Open File que aparece seleccionar el fichero denominado yerba.eq. Este fichero contiene un movimiento del suelo obtenido en el sustrato rocoso de las Isla Yerba Buena en la baha de San Francisco durante el terremoto de Loma Prieta de 1989. Las caractersticas del movimiento del suelo registrado aparecern en la caja Object Motion (los valores de aceleracin pico, paso temporal y frecuencia de
An. VI - P7 - Pgina 15
19.
20.
21
Asigna el movimiento al substrato rocoso introduciendo 16 en el campo Layer dentro del campo Object Motion Location. Especifica este movimiento como un movimiento aflorante marcando la caja apropiada (Outcrop). Selecciona Yes en la caja Animation (Figura 16). Esto inducir a EduShake a calcular la respuesta en las localizaciones necesarias para llevar a cabo una animacin de la respuesta del suelo que podr ser vista en el Output Manager. Con esto se finaliza la introduccin de las caractersticas del suelo y del movimiento de entrada. Seguidamente archiva los datos introducidos en un fichero .dat Documentacin de la entrada de datos
22.
23.
24
Haz click en el botn Report. Se abrir un procesador de textos con una plantilla que muestra todos los datos de entrada en forma tabulada (Figura 16). Este informe puede imprimirse directamente desde EduShake, o archivarse (en formato Rich Text) para ser posteriormente editado con un procesador de textos ms potente.
An. VI - P7 - Pgina 16
Figura 16. Informe sobre los datos de entrada en el procesador de textos de EduShake.
25.
Los grficos tambin se pueden copiar a este informe. Vuelve a la pantalla Input Motion y dibuja el grfico de aceleracin del movimiento de entrada frente al tiempo. Haz click en el botn Copy to Report y vuelve a la pantalla Report. Observars como dicho grfico se ha copiado en la posicin en la que estaba el cursor. Haz lo mismo con el resto de grficos de la pantalla Input Motion. Ahora archiva este informe como un fichero Rich Text, con el nombre informe.rtf.
26.
3.3 Ejecucin del anlisis 27. 28. Para ejecutar el anlisis haz click en el botn Solution Manager En este momento se analizar el fichero de datos que nosotros deseemos. La extensin por defecto de los ficheros de datos EduShake es .dat aunque se pueden usar ficheros con otras extensiones. En esta experiencia seleccionaremos el fichero de datos que grabamos en el paso 23 de esta experiencia. EduShake llevar a cabo, entonces, el anlisis seleccionado. Aparecer una pantalla indicativa del proceso iterativo que se est haciendo, en cuya esquina inferior izquierda podr leerse Running mientras dura el proceso de ejecucin y Finished cuando se complete el anlisis. Cuando finalice el anlisis, haz click en Si y se crear un fichero con el mismo nombre del fichero de datos y extensin .lyr, en el que se hallarn archivados todos los resultados (Figura 17).
An. VI - P7 - Pgina 17
29.
3.4 Obtencin y representacin de los resultados 30. 31. Para visualizar los resultados obtenidos, pincha en Output Manager y abre el fichero .lyr creado en el paso 29. Hay numerosas opciones para ver los resultados obtenidos en el Output Manager. Los tipos de grficos primarios estn organizado en un conjunto de seis campos: Ground Motion Plots, Stress and Strain Plots, Response Spectrum Plots, Depth Plots, Other Parameters, y Animation (Figura 18). Los grficos de movimiento del suelo permiten dibujar historias temporales y espectros de Fourier en la parte superior de las capas seleccionadas en el paso 16 de esta experiencia (Figura 19). Es posible dibujar uno o ms movimientos en el mismo grfico sin ms que seleccionar las cajas correspondientes. Para esta experiencia haz un par de grficos con estos movimientos y copialos al informe que habas creado en el paso 26. Los grficos de esfuerzo y deformacin permiten dibujar historias temporal y espectros de esfuerzos de cizalla y deformacin de cizalla - en nuestra experiencia podemos representar, por ejemplo, los correspondientes a las capas 2 y 11 (Figura 20).
32.
33.
An. VI - P7 - Pgina 18
Figura 19. Representacin de la historia temporal de la aceleracin para la parte superior de la Capa 1 (amarillo), la Capa 11 (negro) y la Capa 16 (rojo).
34.
Los grficos de espectros de respuesta permiten dibujar uno o ms espectros de respuesta segn las capas definidas en el paso 16. En esta experiencia, obtn los
An. VI - P7 - Pgina 19
Figura 20. Representacin de la historia temporal de la deformacin de cizalla para la Capas 2 (amarillo) y la Capa 11 (negro).
35.
Los dibujos de profundidad pueden usarse para dibujar la variacin de diversos parmetros con la profundidad para uno o ms movimientos del suelo. Puesto que nuestro anlisi usa slo un movimiento de entrada, marca la caja correspondiente a Motion 1' e intenta dibujar algunas de estos parmetros con la profundidad (Figura 22). La opcin correspondientes a Other Parameters ofrece la oportunidad de calcular una variedad de parmetros del movimientos del suelo tiles, en la parte superior de cualquiera de la capas definidas. Resalta la capa de inters y haz click en el botn Calculate para ver los valores numricos de estos parmetros. Copia dicha informacin al informe (Figura 24). La opcin Animation permite ver la variacin de los desplazamientos horizontales tanto con la profundidad como con el tiempo para cada movimiento de entrada. Haz click en el botn Plot.
36.
37.
An. VI - P7 - Pgina 20
Figura 21. Espectro de respuesta en aceleracin para la Capa 1 con un amortiguamiento del 5% (amarillo), 10% (negro) y 20 % (azul).
An. VI - P7 - Pgina 21
Figura 24. Animacin del desplazamiento del perfil de suelo con la profundidad debido al movimiento de entrada.
An. VI - P7 - Pgina 22
39.
3.5 Realizacin de la memoria de la prctica 40. Revisa el informe que has recopilado, utiliza un procesador ms potente para escribir un informe detallado del proceso y las conclusiones que se pueden sacar de los resultados obtenidos.
An. VI - P7 - Pgina 23