Análisis de Interacción Suelo-Estructura
Análisis de Interacción Suelo-Estructura
Facultad de ingenierías
PROYECTO DE INVESTIGACION
“Software FLUSH”
“Software FLUSH”
I. ÍNDICE
BREVE HISTORIA .................................................................................... 4
FLUSH ................................................................................................ 5
1
AGRADECIMIENTO
El presente trabajo de investigación agradecemos a Dios por bendecirnos
porque hiciste realidad este sueño anhelado. A nuestros padres que siempre
estuvieron presentes en cada momento de nuestras vidas y con su esfuerzo
lograron que cumplamos nuestras metas y/o sueños de estudiar nuestra carrera
universitaria. a la universidad tecnología de los andes por darnos la oportunidad
de estudiar y ser profesionales.
También nos gustaría agradecer a nuestros docentes que durante toda
nuestra carrera profesional aportaron sus conocimientos obtenidos en su
experiencia profesional para poder lograr en nosotros unos buenos profesionales
a futuro. Son muchas las personas que han formado parte de nuestras vidas en
el ámbito laboral y profesional a las que nos encantaría agradecerles su amistad,
consejos, apoyo, ánimo y compañía en los momentos más difíciles y más alegres
de nuestra vida. Algunas están aquí con nosotros y otras en nuestros recuerdos
y en nuestro corazón, sin importar en donde estén queremos darles las gracias
por formar parte de nuestras vidas, por todo lo que nos han brindado y por todas
sus bendiciones.
2
IV. INTRODUCCIÓN
3
BREVE HISTORIA
4
Estructura, para lo cual se analizaron diferentes condiciones de suelo,
trabajando con perfiles de suelos homogéneos y perfiles de suelos que
presentan un incremento de la rigidez con la profundidad, y edificios en pórticos
y muros de cortante representativos de los sistemas estructurales empleados en
la ciudad de Medellín. En total se analizaron más de 250 casos con el método
completo, empleando el programo FLUSH, para lo cual fue necesario desarrollar
interfaces para la entrada de datos ( INPUT FLUSH 2005) y para el
procesamiento de los resultados (OUTPUT FLUSH 2005). Como resultado de
los análisis se encontró que los fenómenos de interacción más notorios no
necesariamente ocurren en edificios esbeltos sobre suelos blandos, como
considera la NSR-98, sino cuando ocurren coincidencias entre el periodo de
vibración del primero o algún modo superior de vibración de la estructura y el del
suelo.
FLUSH
FLUSH: Es un programa de computadora para el análisis trimendimensional
aproximado de problemas de interacción suelos estructura. FLUSH, emplea
modelos con elementos finitos que permiten representar el suelo mediante
elementos sólidos y la cimentación y la estructura por medio de elementos viga,
permite además considerar fronteras viscosas que evitan la concentración de
ondas en los bordes de los modelos y simula el efecto tridimensional. Este
programa realiza los análisis solucionando las ecuaciones del sistema en el
dominio de las frecuencias modelando la rigidez del suelo y el amortiguamiento
mediante el módulo de rigidez complejo. La no linealidad del suelo la simula
mediante el método lineal equivalente.
5
4. Una capacidad aproximada de 3D que hace posible realizar ciertos análisis
de interacción estructura-suelo-estructura esencialmente al mismo costo que el
análisis 2D.
5. Límites de transmisión.
6. Un solucionador de ecuaciones fuera del núcleo que reduce los problemas
de tamaño del núcleo.
7. Desconvolución interna. Esta función elimina la necesidad de realizar un
análisis de respuesta del sitio independiente para determinar los movimientos de
la base rígida.
8. Funciones convenientes de manejo de archivos que brindan capacidad de
reinicio y un registro permanente para la recuperación posterior de los detalles
de la solución.
9. Historiales de tiempo de salida impresos, trazados o almacenados de
momentos de aceleración y flexión.
10. Cálculo de fuerzas cortantes máximas en elementos viga.
11. Espectros de respuesta de velocidad de aceleración impresos o
almacenados.
12. Trazado de funciones de amplificación de Fourier.
13. Corrección de línea de base opcional de movimientos de entrada y salida
1.2 Alcances sobre el FLUSH
Actualmente durante las fases de elaboración y ejecución de los proyectos
en edificaciones, Los edificios en muros estructurales tienen una altura por piso
de 2.5 m para todos los pisos y constan de 7 ejes de carga en una dirección y 9
ejes de carga en la otra. Para estas edificaciones se tiene una alta
densidad de muros de diversas longitudes con un espesor de 0.13 m. Las losas
de piso se definieron como losas macizas armadas en dos direcciones con un
espesor de 0.10 m.
Para simplificar los modelos a utilizar en el programa FLUSH, se realizó una
condensación estructural que consiste en transformar la estructura
tridimensional en un vástago equivalente con masas concentradas, con unas
propiedades mecánicas y geométricas tales que el vástago equivalente tiene
idéntico comportamiento ante cargas horizontales y similar respuesta inercial
Considerando estos resultados, se propuso la ejecución de un estudio
paramétrico destinado a identificar las variables que tienen un mayor efecto en
6
los resultados de la Interacción Sísmica Suelo Estructura, para lo cual se
analizaron las siguientes variables en los perfiles de suelo:
7
Rigidez del suelo, se consideraron suelos homogéneos y suelos cuya
rigidez varía con la profundidad.
Espesor del perfil de suelo.
Impedancia con la roca sobre la cual está apoyado el perfil de suelo.
En cuanto a las estructuras se consideraron edificios representativos de las
soluciones estructurales normalmente empleadas en la Ciudad de Medellín. Se
analizaron edificios con alturas entre 5 y 20 pisos en pórticos o en muros de
cizalladora; apoyados en cimentaciones superficiales.
El método idóneo para realizar estos análisis es el método Completo,
mediante modelos de elementos finitos analizados con el programa FLUSH, para
lo cual fue necesario desarrollar dos interfaces que permitieran hacer eficiente la
utilización del programa FLUSH, como son INPUT FLUSH 2005, para los datos
de entrada, y OUTPUT FLUSH 2005, para el pos procesamiento de los
resultados.
Este análisis paramétrico se realizó como parte de la línea de investigación
en Interacción Suelo Estructura que se desarrolla en la Escuela de Ingeniería a
través de dos Trabajos Dirigidos de Grado.
Pero a veces, como con las bases de datos, desea asegurarse de que los
datos lleguen al almacenamiento. En estos casos, puede vaciar los búferes para
8
que sea menos probable que un corte de energía deje los datos en un estado
cuestionable.
Ilustración 1 Coeficiente de interacción sísmica máxima en función de la rigidez del suelo para los
diferentes modelos de suelo (MERINO, 2009).
9
Ilustración 2 Axial circunferencial m’aximo según diferentes códigos de la rigidez del suelo para los
diferentes modelos de suelo. (GUARÍN, 2009)
El código FLUSH
También permite el estudio dinámico de problemas de ISE en deformación
plana, como es el caso tratado, permitiendo además la resolución del problema
de Campo-Libre, así como estudios de deconvolución y demás disponibilidades,
presenta una serie de limitaciones derivadas de las capacidades
computacionales propias de la época en la que se desarrolló (años 70). Por ello,
con el desarrollo y puesta a punto del código ESES, se propone una alternativa
a los estudios de ISE en condiciones sísmicas, incorporando un nuevo Borde
Transmisor, así como nuevas capacidades que el programa posee, como son:
• Permite leer las mallas generadas mediante cualquier pre-procesador de
elementos finitos, siguiendo tan sólo unas pautas concretas, lo cual le
proporciona una gran versatilidad a la hora de generar la geometría, así como
una menor inversión de tiempo en el desarrollo de la misma.
• Los elementos finitos para el mallado de la zona de estudio, son elementos
cuadriláteros isoparamétricos de 8 nodos, en lugar de los tradicionales de 4
nodos como ocurre en FLUSH, con lo cual se consigue una mayor precisión con
menor esfuerzo computacional.
10
• Permite modelar la estructura con este mismo tipo de elementos, mucho
más precisos que los elementos viga empleados tradicionalmente en otros
programas de ISE como era el caso de FLUSH, dando la opción de calcular,
mediante una subrutina en esta Tesis desarrollada (subrutina ANGEL), los
esfuerzos axiles, cortantes y flectores máximos en los elementos del
revestimiento del túnel mediante la integración numérica de las tensiones a lo
largo de la sección. En otras palabras, evita la inconsistencia que, por ejemplo,
se produce en el programa FLUSH y en otros códigos comerciales (ABAQUS y
FLAC) en los cuales los elementos viga que modelan la estructura poseen,
lógicamente, grados de libertad de giro en un intento burdo de unir elementos 2-
D con elementos 1- D.
• Analiza todo el espectro de frecuencias del desarrollo en serie de Fourier
sin interpolar los resultados entre determinadas frecuencias como, por ejemplo,
hace el programa FLUSH, lo cual proporciona una mejor aproximación a los
resultados por no despreciar ni interpolar para frecuencias intermedias
previamente elegidas. Asimismo tiene las siguientes desventajas:
• El código se reduce al estudio de problemas en deformación plana, y no
permite el análisis de problemas axisimétricos, aunque para la presente
investigación no significaba un inconveniente.
• No permite un análisis lineal equivalente (EQLM) como hace el programa
FLUSH, aunque se podría implementar con una cierta relativa facilidad esta
forma de proceder. Lo anterior es una mera técnica de aproximación sin base
matemática, basada solamente tratar de representar la pérdida de rigidez que
experimenta el suelo con la deformación angular, por lo que no ha resultado de
interés su implementación en el código desarrollado.
• No realiza estudios de deconvolución, en los cuales se obtiene el
acelerograma en la base del modelo, a partir de uno definido en algún punto del
mismo en la condición de Campo-Libre; sin embargo, el procedimiento de
deconvolución es más útil en situaciones donde se conoce el perfil de suelo, así
como el registro de aceleraciones, algo no siempre disponible.
11
Ilustración 3 Información básica y extensiones de archivos Fuente
[Link]
12
1.3.1 Comparación tenso-deformacional con el código
FLUSH
13
Borde Consistente desarrollado por Waas [96] y además permite emplear
modelos de suelo viscoelásticos, así como aplicar el Método Lineal Equivalente
(EQLM). Para la modelización del sistema suelo-estructura posee las siguientes
características: • Elementos cuadriláteros isoparamétricos de 4 nodos para la
modelización del suelo.
Ilustración 5 Coeficiente de iteración sísmica máximo según diferentes códigos para un modelo de
suelo VEL
14
Ilustración 6 Diferencia (%) de los axiales circunferenciales máximos con los modelos de suelos NL y
EQLM respecto del modelo de suelo VEL
La diferencia entre los valores de los esfuerzos máximos que predicen los
códigos utilizados ayuda a ver la diferente estimación de cada código en función
de sus características. Así por ejemplo, la diferencia entre los esfuerzos
estimados con FLUSH y con ESES, llega a ser de casi del 20 %, lo cual viene a
confirmar la mejora introducida en ESES con respecto a FLUSH, pues aun
empleando ambos códigos un Borde Transmisor exacto, y habiéndose
demostrado que estiman las mismas deformaciones al modelar el problema
utilizando el código ESES (la estructura se modela con elementos
isoparamétricos de 8 nodos), se obtienen unos resultados más conservadores,
15
que a la par podrían ser asumidos como más exactos por las mejoras
introducidas en el código desarrollado. De todo lo expuesto se llega a la
conclusión que son más fiables los códigos que trabajan en el dominio de la
frecuencia por emplear un borde más exacto y por el tipo de amortiguamiento
que permiten introducir.
Para simplificar los modelos a utilizar en el programa FLUSH, se realizó una
condensación estructural que consiste en transformar la estructura
tridimensional en un vástago equivalente con masas concentradas, con unas
propiedades mecánicas y geométricas tales que el vástago equivalente tiene
idéntico comportamiento ante cargas horizontales y similar respuesta inercial
ante excitaciones dinámicas que la estructura original.
En este proceso se está suponiendo que se puede concentrar la masa
del edificio en los nodos que definen los elementos “viga”, esto sería
estrictamente valido para el caso de edificios en pórticos, puesto que la masa de
los pisos, dada por las losas y las vigas de amarre, es mucho mayor que la
dada por las columnas. En el caso de los edificios en muros estructurales de
cortante, la masa de los muros es comparable con la masa total por piso; por
tanto es necesario distribuir la masa entre los nudos y los elementos viga.
Ilustración 7Axial circunferencial máximo según diferentes códigos para un modelo de suelo VEL
16
pasos de tiempo muy pequeños o en muchas frecuencias, como es el caso del
código FLUSH, que sólamente calcula la solución del problema para un
determinado número de frecuencias, interpolando los resultados para
frecuencias intermedias; esto tiene el inconveniente de no considerar
componentes frecuenciales que pudieren resultar de importancia.
La mejora introducida en el código ESES, respecto de la limitación citada, es
que es capaz de estudiar todo el intervalo discretizado de frecuencias dadas por
la FFT, por ello permite extraer la influencia de introducir una frecuencia de corte,
pues cuanto menor es el número de frecuencias a analizar, menor será el tiempo
de simulación, y por tanto se podrán realizar modelos más grandes o de menor
tamaño de malla. Para ver esta influencia, se compararán los resultados de
calcular con todo el espectro (fmax=50 Hz), con los resultantes de tomar
frecuencias de corte de 10, 15 y 25 Hz, para un modelo de túneles a una
profundidad media en un estrato de propiedades uniformes según las dadas en
la Tabla 6.1, de radio R=4 m y e/R=10 % y para un amortiguamiento interno del
5 %. se recogen los resultados obtenidos en términos de axiles circunferenciales
máximos, considerando una frecuencia de corte en el análisis, respecto al
análisis de todo el espectro de frecuencias de Fourier.
17
al estudio de todo el espectro de frecuencias de Fourier, cuanto más rígido es el
estrato, mayor es el error que se comete al despreciar las frecuencias altas, dado
que su frecuencia natural es mayor.
1.3.3 La deconvolución
Es una operación matemática. Se usa en restauración de señales para
recuperar datos que han sido degradados por un proceso físico que puede
describirse mediante la operación inversa, una convolución.
Tabla 1
18
V. CONCLUSIONES
Los espectros de la microzonificación sísmica tienen ordenadas
espectrales muy inferiores a las que se encuentran al realizar un análisis
de la amplificación del movimiento sísmico en la masa de suelo por
medio de un programa como el SHAKE o el FLUSH. Esta observación
es válida para los periodos cortos del espectro, es decir, los espectros
normativos menosprecian el contenido de alta frecuencia del movimiento
que llega a la cimentación.
La No Linealidad de los suelos donde la velocidad de onda de corte
aumenta con la profundidad toma un papel secundario en cuanto a la
respuesta de estos suelos ante excitaciones sísmicas, para los cuales su
comportamiento estará determinado principalmente por la
configuración y magnitud de las rigideces de sus estratos. Sin embargo
esta particularidad no se cumple o no es tan significativa para suelos
blandos y/o de poco espesor.
En lo referente al modelado numérico del problema, los tipos de Bordes
Transmisores disponibles, aquellos actualmente incorporados en los
programas comerciales, no son exactos; sin embargo, si se sitúan lo
suficientemente alejados de la zona en estudio, aportan resultados
aproximados.
Mientras que el código FLUSH estudiaba el problema en ciertas
frecuencias, interpolando en las frecuencias intermedias, se ha
demostrado que el hecho de despreciar las frecuencias altas puede
afectar a los resultados, tanto más cuanto más rígido sea el estrato debido
a que su frecuencia natural es mayor. Por ello, el código ESES permite el
estudio del problema en todo el espectro de frecuencias de Fourier.
19