Trabajo Final Analisis No Lineal Estatico Pushover - Copia

UNIVERSIDAD NACIONAL DE COLOMBIA SEDE MANIZALES

FACULTAD DE INGENIERIA Y ARQUITECTURA ESPECIALIZACION EN ESTRUCTURAS

VULNERABILIDAD SISMICA

ANALISIS NO LINEAL ESTATICO "PUSHOVER" USANDO SAP2000

Trabajo Presentado por: CRISTIAN DAVID BOLIVAR BUITRAGO

Trabajo Presentado a:

Doctor CARLOS ALBERTO BERMUDEZ MEJIA

MANIZALES, (CALDAS), JUNIO DE 2012

ANALISIS NO LINEAL ESTATICO “PUSHOVER”

2

INTRODUCCION

A continuación se desarrolla un informe detallado de la valoración sísmica de un edificio modelado para elaborar el presente trabajo usando el método no lineal estático de carga incremental más conocido como “Pushover”, el cual es una técnica eficiente para estudiar la capacidad, resistencia – deformación de una estructura bajo una distribución esperada de fuerzas inerciales además de determinar la respuesta máxima por medio del “Método del Espectro de Capacidad”, evaluando el máximo desplazamiento, por medio de la intersección entre la curva de capacidad y el espectro de demanda reducido. No obstante un análisis elástico da una buena indicación de la capacidad elástica de la estructura, las cargas siempre son directamente proporcionales a las deformaciones lo cual no permite pronosticar la respuesta inelástica que puede ocurrir durante una solicitación sísmica, tampoco puede predecir el mecanismo de falla y tomar en cuenta la redistribución de fuerzas durante la fluencia continua de los elementos. El análisis inelástico ayuda a demostrar como realmente trabaja el edificio, por medio de la identificación de los modos de falla y el colapso progresivo. El único objetivo de los dos métodos es determinar el punto de desempeño “Performed Point” o máxima respuesta de desplazamiento. Con este punto se pueden obtener las Probabilidades de daño en la Estructura usando las curvas de fragilidad, siendo una herramienta útil para evaluar edificios existentes y lograr objetivos específicos de rehabilitación. La demanda es la representación del movimiento del suelo, mientras que la capacidad es la habilidad que posee la estructura para resistir la demanda sísmica.


3

OBJETIVOS

Obtener la curva de capacidad lateral más allá del rango elástico.

Obtener la formación secuencial de mecanismos y fallas en los elementos.

Determinar la capacidad lateral de la estructura.

Cuales elementos serán mas susceptibles a fallar primero.

Determinar la ductilidad local de los elementos y global de la estructura.

Verificar el concepto de vigas débiles y columnas fuertes.

Verificar la Degradación global de la resistencia.

Verificar los desplazamientos relativos inelásticos.

Verificar los criterios de aceptación a nivel local de cada elemento.


4

1. DESARROLLO

1.1 Procedemos a modelar un edificio el cual a de cumplir lo exigido en la norma sismo resistente 2010 (NSR-10), además de incluir lo que se presenta a continuación.

Modelo Estructural.

Aplicación de los patrones incrementales de carga lateral

Determinación de la Capacidad Lateral Global de la Estructura.

Determinación de la demanda inelástica que se espera en la respuesta máxima de la estructura.

Evaluación de los elementos estructurales cuando ocurre la máxima respuesta en la estructura.

Evaluación de la Respuesta Global de la Estructura. 2 MODELO ESTRUCTURAL. a) Programa Utilizado para el Análisis Para el presente trabajo utilizamos el programa SAP2000 V14. El modelo consiste en una estructura compuesta por un sistema de pórticos de concreto. El análisis fue efectuado en tres dimensiones. Con relación a las vigas y columnas hemos desarrollado un modelo de rotulas plásticas según los criterios vistos en clase que consisten en inducir rotulas mediante Sap 2000 al 5% y al 95% de las longitudes de los elementos. El patrón incremental de cargas laterales parte de las cargas estáticas según las características geológicas y de zonificación asumidas para este ejercicio. A su vez, la geometría del modelo se muestra en la figura 1. Se asume que el modelo ha sido analizado y diseñado previamente antes de iniciar este ejemplo. Este programa permite considerar la degradación de la resistencia en los elementos luego de que estos pasan su rango de fluencia y tienen una alta demanda de ductilidad. Tampoco limita el comportamiento hasta la pos fluencia si no que considera la degradación en ciclos sucesivos en la resistencia flexional del elemento lo cual es una gran ventaja que nos permite obtener un comportamiento más real de la estructura. b) Materiales

Resistencia del concreto: 3000 psi -21 Mpa

Resistencia del Acero: 60 Ksi-4200 Mpa


5

c) Sistema Estructural Hemos modelado la estructura en tres dimensiones. La base de la estructura ha sido restringida asumiendo un empotramiento perfecto. d) Espectro de diseño DATOS DE ENTRADA ANALISIS

GRUPO DE USO = II (NSR-10 A.2.5) T (s) Sa (g) Sv (m/s) Sd (m)

COEFICIENTE DE IMPORTANCIA I = 1,10 (NSR-10 A.2.5) 0 0,00 0,894 0,000 0,000

CIUDAD O MUNICIPIO Manizales (NSR-10 A.2.3) 1 0,03 0,894 0,021 0,000

ZONA DE AMENAZA SISMICA Alta (NSR-10 Tabla A.2.3-2) 2 0,06 0,894 0,052 0,000

VALORES DE ACELERACION PICO EFECTIVA ==> 3 0,09 0,894 0,093 0,001

Aa= 0,25 (NSR-10 Tabla A.2.3-2) 4 0,12 0,894 0,143 0,003

Av= 0,25 (NSR-10 Tabla A.2.3-2) 5 0,15 0,894 0,204 0,005

CLASIFICACION DEL SUELO TIPO= D (NSR-10 A.2.4) 6 0,25 0,894 0,349 0,014

VALORES COEFICIENTES EFECTOS LOCALES ==> 7 0,28 0,894 0,391 0,017

Fa= 1,30 (NSR-10 Tabla A.2.4-3) 8 0,42 0,894 0,587 0,039

Fv= 1,90 (NSR-10 Tabla A.2.4-4) 9 0,56 0,894 0,782 0,070

10 0,70 0,894 0,977 0,110

T0 TC TL 11 0,71 0,883 0,977 0,111

0,146 0,702 4,560 12 0,79 0,794 0,977 0,124

13 1,05 0,599 0,977 0,164

14 1,30 0,481 0,977 0,204

15 1,56 0,401 0,977 0,245

16 1,82 0,345 0,977 0,285

17 2,08 0,302 0,977 0,325

18 2,33 0,269 0,977 0,366

19 2,59 0,242 0,977 0,406

20 2,85 0,220 0,977 0,446

21 3,11 0,202 0,977 0,487

22 3,36 0,186 0,977 0,527

23 3,62 0,173 0,977 0,567

24 3,88 0,162 0,977 0,608

25 4,56 0,138 0,977 0,715

26 4,81 0,124 0,926 0,715

27 5,06 0,112 0,881 0,715

28 5,31 0,101 0,839 0,715

29 5,56 0,092 0,801 0,715

30 5,81 0,085 0,767 0,715

31 6,06 0,078 0,735 0,715


6

e) Cargas

ELEMENTO Peso Unidad

Muros 300 kg/m2

Pisos 90 kg/m2

Mortero Pega 66 kg/m2

Cielo raso 10 kg/m2

Caseton 40 kg/m2

Vigueteria 188 kg/m2

Torta de concreto 120 kg/m2

TOTAL 814 Kg/m2

ESTIMACIÓN CARGA MUERTA

f) Dimensiones edificio.

DESCRIPCIÓN DIMENSIÓN Un

Numero de pisos 4 Pisos

Altura entre pisos 2,8 m

Distancia entre ejes

en plano ZX5 m

Numero de vanos

en plano ZX4 un

Distancia entre ejes

en plano ZY6 m

Numero de vanos

en plano ZY3 un

g) Derivas.

TABLE: Joint Displacements - Absolute

Joint OutputCase CaseType StepType U1 U2 U3 R1 R2 R3

Text Text Text Text m m m Radians Radians Radians Derivas

100 SISMX LinRespSpec Max 0,023065 0,007939 0,000372 0,000225 0,000613 0,000283 0,003305 CUMPLE





7

Figura 1 modelo tridimensional

Figura 3 Vista frontal Figura 2 Vista lateral izquierda


8

2.1 MODELACIÓN UTILIZANDO SAP2000 Ya que no es un curso de modelación no enfatizaremos en la forma en que el programa hace los cálculos puesto que su algoritmo utiliza los parámetros de ATC 40 y FEMA para el análisis no lineal de estructuras, pero si mostramos gráficos y datos de resultados y explicaremos su contenido.

2.1.1 PATRON INCREMENTAL DE CARGAS LATERALES. A fin de tomar en consideración el efecto de la carga de gravedad en la aplicación de la carga lateral hemos incluido un primer análisis controlado directamente por las fuerzas asumidas, donde se espera los elementos permanezcan en su rango lineal a fin de capturar la perdida de rigidez inicial dado a la aplicación de la cargas de gravedad gradualmente. En dicho análisis hemos tomado el 100% de las cargas de gravedad combinado con el 25% de la carga viva. El programa determina automáticamente la distribución de la carga lateral por piso, la misma se aplicó en las dos direcciones dado que la rigidez en cada dirección no es la misma.

2. OBTENCIÓN CURVA PUSHOVER MEDIANTE SAP 2000

ANALISIS DE RESULTADOS

A. DEFORMACION EN CADA INCREMENTO DE CARGA LATERAL Y SECUENCIA DE FORMACION DE ROTULAS

Figura 4. Figura donde se evidencian las primeras rótulas

PRIMERAS

ROTULAS


9

Figura 4. Cortante en la base vs desplazamiento (Capacidad lateral en dirección X) Curva

Pushover

El primer evento del comportamiento descrito en la grafica nos puede indicar que una columna del primer nivel alcanza la fluencia, indicando un comportamiento dúctil y favorable. El primer cambio de pendiente presente en la grafica nos da a entender un cambio en el comportamiento de la estructura, posiblemente continua la fluencia en los elementos, antes de darle paso a las primeras rótulas las cuales han de generarse como lo muestra el modelo en las columnas del segundo nivel, dándole paso al colapso de estos elementos, fenómeno evidenciado en la grafica, en la zona donde hay un cambio brusco con desplazamiento constante y cortante en la base en aumento.


10

Step Displacement BaseForce AtoB BtoIO IOtoLS LStoCP CPtoC CtoD DtoE BeyondE Total

m Kgf

0 -0,000052 0 408 0 0 0 0 0 0 0 408

1 0,023826 419,90103 406 2 0 0 0 0 0 0 408

2 0,033738 567,75536 354 54 0 0 0 0 0 0 408

3 0,039579 614,60314 312 96 0 0 0 0 0 0 408

4 0,078157 732,99108 232 118 58 0 0 0 0 0 408

5 0,125138 826,69269 200 108 98 0 0 0 2 0 408

6 0,125143 1038,89117 172 136 98 0 0 0 2 0 408

7 0,177433 1345,6711 150 138 96 14 6 0 4 0 408

8 0,177438 1480,43614 148 130 98 20 8 0 4 0 408

9 0,180487 1499,63046 146 126 102 20 8 0 6 0 408

10 0,180492 1541,60286 144 120 110 20 8 0 6 0 408

11 0,232283 1812,4432 128 114 110 34 14 0 8 0 408

12 0,232283 1812,4432 128 114 110 34 14 0 8 0 408

13 0,232283 1812,44321 128 114 110 34 14 0 8 0 408

14 0,233109 1815,6029 126 116 108 36 14 0 8 0 408

15 0,19969 1714,57887 122 118 110 36 14 0 8 0 408

TABLE: Pushover Curve - EMPUJE_X

Los criterios de aceptación se usan para demostrar el estado de las rotulas asignadas, cuando se muestran los resultados del análisis, estos criterios no afectan el comportamiento de la estructura. Estos valores de aceptación se refieren a deformaciones, rotaciones y desplazamientos.

3. ESPECTRO DE CAPACIDAD MAXIMA RESPUESTA ESPERADA EN LA ESTRUCTURA BAJO LA DEMANDA SISMICA SELECCIONADA. ESPECTRO DE CAPACIDAD OBTENIDO CON SAP 2000

PROCESO NSP (Análisis No Lineal Estático Pushover) La Técnica Pushover es Apropiada para:


11

Obtener la curva de capacidad lateral más allá del rango elástico. Obtener la formación secuencial de mecanismos y fallas en los elementos concepto General de la Técnica Pushover:

Consiste en un proceso sucesivo de análisis estático incremental que toma en cuenta la variación de la rigidez en los elementos en cada elemento.

El análisis se efectúa incrementando la carga lateral hasta que la estructura alcanza ciertos limites de desplazamientos o se vuelva inestable

Existen diferentes técnicas para estimar la máxima respuesta en una estructura, en este caso se trabajado el “METODO DEL ESPECTRO DE CAPACIDAD”, método que radica en convertir las curvas de capacidad lateral global y la curva de demanda (espectro reducido) a coordenadas espectrales, superponerlas en una misma grafica, reducir la demanda tomando en consideración la capacidad adicional de disipar energía que tiene la estructura por medio de los ciclos de histéresis y determinar la intersección entre dichas curvas, que va a representar la máxima respuesta de desplazamiento que se espera en la estructura. La demanda será representada por medio del espectro de respuesta creado para Manizales que es donde hipotéticamente se encuentra la estructura

En la aplicación del método el programa requiere los factores Ca y Cv para crear el espectro de demanda, estos componentes constituyen la aceleración como función de la gravedad según el código UBC-97 (En otras palabras son coeficientes sísmicos) usado en EE.UU. En nuestro caso (para Colombia) debemos determinar los factores como sigue:

Aa = 0.25 Av = 0.25 Fa = 1.30 Fv= 1.90 Según ASCE-07 Sms = Ca(UBC-97) y Sm1= Cv Entonces Sms = Fa*Aa = 0.32 = Ca Sm1 = Fv*Av = 0.47 = Cv Entre otros términos usados seleccionados para la aplicación de este método podemos mencionar:

Según la Tabla 8-4 (ATC-40) Asumiremos una Estructura Tipo B (Promedio entre una estructura nueva y vieja).

Desde la Tabla 8-3 (ATC-40) El βeff Máximo se puede tomar como 40% (Dado que no estamos seguro cuan estables serán los ciclos de histéresis en una estructura existente debemos desarrollar una familia de espectros de demanda para visualizar que ocurre si el amortiguamiento efectivo varia). Para estos fines vamos a usar 5%, 10%, 15%, 20%.


12

El amortiguamiento inherente usado será igual 0.05, asumido para el acero estructural.

El factor Kappa según la Tabla 8-1 para un comportamiento tipo B, es calculado automáticamente por el programa según el amortiguamiento histeretico dado por la curva de capacidad lateral global de la estructura.

A continuación se presenta la tabla de información para calcular la curva pushover.


13

Para visualizar en que paso del análisis la estructura alcanza su respuesta máxima solo hay que ir a la tabla a través del menú File

Step Teff Beff SdCapacity SaCapacity SdDemand SaDemand Alpha PFPhi

m m

0 0,524272 0,050000 0,000000 0,000000 0,006224 0,091163 1.000.000 1.000.000

1 0,524272 0,050000 0,016696 0,244526 0,006224 0,091163 0,836354 1.430.183

2 0,541848 0,070391 0,023822 0,326638 0,006084 0,083415 0,846570 1.418.420

3 0,571719 0,106586 0,028379 0,349519 0,006010 0,074018 0,856433 1.396.472

4 0,762829 0,228073 0,059344 0,410548 0,007546 0,052206 0,869568 1.317.885

5 0,902808 0,243184 0,095361 0,470996 0,008706 0,043001 0,854861 1.312.798

7 0,796759 0,118948 0,131110 0,775080 0,010532 0,062263 0,849307 1.332.587

9 0,825209 0,110205 0,131809 0,892416 0,012243 0,082893 0,845593 1.353.710

11 0,770171 0,050000 0,163237 1.083.716 0,011911 0,079075 0,818857 1.368.004

12 0,771096 0,050000 0,163237 1.083.716 0,011911 0,079075 0,818438 1.369.697

13 0,760394 0,050000 0,163237 1.083.716 0,011911 0,079075 0,823145 1.378.090

14 0,778702 0,058040 0,163717 1.085.862 0,011897 0,078910 0,814550 1.423.292

TABLE: Pushover Curve Demand Capacity - ATC40 - EMPUJE_X

Usado para buscar el paso correspondiente a la respuesta máxima

Respuesta Max


14

De la tabla anterior podemos deducir que la máxima respuesta se encuentra entre los pasos 5 y 7, debido a que se evidencia un salto muy significativo entre estos dos pasos.

4. PUNTO DE DESEMPEÑO OBTENIDO CON LOS MÉTODOS VISTOS EN CLASE (ATC40 Y RISK UE)

Definición

Desempeño: Una vez se ha definido la curva de capacidad y los desplazamientos de

demanda, se puede verificar el punto de desempeño o punto de respuesta máxima en la estructura. Este punto es aproximadamente igual al obtenido por el método de los coeficientes el cual no se explica en este trabajo. MÉTODO RISK UE: A continuación se expone tabla donde se encuentran los valores a partir de los cuales se grafica la curva que intercepta la curva de capacidad.

Sdi µ Saei Ti Rui Sai

0,01 1 0,894 0,2121669 1 0,894

0,015 1,5 0,894 0,25985032 1,21654193 0,73486986

0,02 2 0,894 0,3000493 1,50008217 0,59596735

0,025 2,5 0,894 0,33546532 1,83866329 0,48622279

0,03 3 0,894 0,36748384 2,22494614 0,40180748

0,035 3,5 0,894 0,39692792 2,65386631 0,33686701

0,04 4 0,894 0,42433379 3,12166895 0,28638527

0,045 4,5 0,894 0,45007395 3,62543138 0,24659134

0,05 5 0,894 0,4744196 4,16279733 0,21475943

0,055 5,5 0,894 0,49757547 4,73181606 0,1889338

0,06 6 0,894 0,51970063 5,33083861 0,16770345

0,065 6,5 0,894 0,54092157 5,95844772 0,15003908

0,07 7 0,894 0,56134084 6,61340841 0,13517992

0,075 7,5 0,894 0,58104297 7,2946322 0,12255587

0,08 8 0,894 0,6000986 8,00115035 0,11173393

0,081 8,1 0,894 0,60383757 8,14541124 0,10975505


15

0,000

0,100

0,200

0,300

0,400

0,500

0,600

0,700

0,800

0,900

1,000

0,000 0,020 0,040 0,060 0,080 0,100 0,120 0,140

Sa

Sd

RISK UE

ESPECTRO

CURVA-CAPACIDAD

Sdi_VS_Sai

Grafica vs Curva Risk UE

Al graficar los datos obtenidos en la tabla anterior obtenemos la curva con la cual debemos interceptar la curva de capacidad para así obtener el punto de desempeño para el método RISK UE. A continuación exponemos los valores obtenidos para este método tanto de manera grafica como de manera analítica.

Sdi= 0,035

Sai= 0,3368

PUNTO DE DESEMPEÑO METODO RISK - UE (VALORES

DE INTERCEPCIÓN DE LAS GRÁFICAS)

Formulas utilizadas


16

MÉTODO ATC40 Proceso A Esta es la aplicación más directa, es verdaderamente iterativo pero, se basa en formulas que pueden ser fácilmente programadas en hojas de cálculo. Este es más bien un método analítico que grafico siendo el método con la aplicación más directa, es el más recomendado.

Calculo Usando el Proceso A. Este proceso de iteración puede ser realizado a mano o

sobre una hoja de Excel para hacer converger el punto de desempeño. 1. Desarrollar el espectro elástico de un 5%, apropiado para la localización.

2. Transformar la curva de capacidad a espectro de capacidad.

3. Graficar las dos curvas en un mismo grafico.

4. Seleccionar un punto asumido inicial de desempeño por el método “Aproximación de igual desplazamiento”, ver figura.


17

5. Desarrollar la representación bilineal del espectro de capacidad

6. Calcular el factor los factores de reducción espectral.

7. Desarrolle el espectro de demanda reducido y graficar en el mismo gráfico.


18

8. Determine si el espectro de demanda intercepta el espectro de capacidad en el punto, api, dpi (estos son los puntos asumidos), de lo contrario verifique si el desplazamiento en el punto de intersección di, está dentro de la tolerancia aceptable del dpi (0.95dpi < di < 1.05dpi). (La siguiente grafica muestra la Tolerancia Entre el Punto Asumido y el Punto de

Intersección).

9. Si el espectro de demanda no intercepta el espectro de capacidad dentro de la tolerancia, se selecciona un nuevo valor api, dpi y se regresa al paso 5.

10. Si el espectro de demanda intercepta el espectro de capacidad dentro de la tolerancia aceptable, entonces el punto asumido api, dpi será el punto de desempeño (ap, dp), y el

desplazamiento dp representa el máximo desplazamiento que se espera en el terremoto.


19

GRAFICA OBTENIDA SIGUIENDO EL PROCEDIMIENTO ANTERIORMENTE MENCIONADO:

0,000

0,100

0,200

0,300

0,400

0,500

0,600

0,700

0,800

0,900

1,000

0,000 0,020 0,040 0,060 0,080 0,100 0,120 0,140

Sa

Sd

METODO ATC - 40

ESPECTRO

CURVA-CAPACIDAD

Sdi_VS_Sai

ESPECTRO REDUCIDO

Comparación resultados obtenidos con los métodos desarrollados y el programa SAP

Dpi Api Sdi Sai Sd Sa

0.026 0.35 0.035 0.33 0.081 0.47

METODO ATC - 40 METODO RISK - UE SAP 2000

El comportamiento sísmico de un edificio sometido a un terremoto se puede cuantificar mediante su respuesta máxima, que se define mediante el punto de capacidad por demanda, también conocido como punto de desempeño (SDpd, SAPD). Este punto se obtiene mediante la intersección entre el espectro de capacidad y el espectro de


20

demanda. El espectro de demanda se obtiene a partir del espectro de respuesta elástico con 5% de amortiguamiento, reducido para niveles mayores de amortiguamiento efectivo.

El movimiento del suelo durante un terremoto produce un patrón complejo de desplazamientos en la estructura que puede variar con el tiempo. Para una estructura dada y un movimiento del suelo, la demanda de desplazamiento es un estimado de la máxima respuesta esperada del edificio durante el movimiento del terreno.

Los espectros de capacidad obtenidos reflejan el comportamiento frágil de este sistema estructural, pudiéndose observar que, para el modelo aquí estudiado, el desplazamiento último obtenido es bajo para este tipo de edificios, pero difieren de un método a otro llevándonos a pensar cual de todos es el correcto o que influyó en uno o en otro método para que los resultados no fueran los esperados o al menos aproximados entre sí. A partir del desplazamiento último la estructura falla debido al mecanismo plástico desarrollado durante el análisis.

Los métodos RISK UE Y ATC – 40 que para nuestro caso muestran valores algo semejantes entre sí, son el resultado de corregir mediante iteraciones valores que en primera medida no se acomodaban a los criterios de construcción grafica implementados y obtenidos a partir de formulas y que a su vez nacen de la información obtenida durante la modelación en el programa SAP 2000, y que a su vez no son datos similares a los obtenidos en el programa pudiendo este fenómeno obedecer a procedimientos de modelación de la estructura no regulares desde el punto de vista de la teoría estructural, pero que a manera de instrucción a la hora de comprender la esencia de los métodos puede ser un instructivo a la hora de proporcionarnos una comprensión de los tres métodos utilizados, máxime no haber obtenido unos resultados equiparables entre si para el propósito de la vida profesional pero brindándonos el entrenamiento necesario para comprender la finalidad de los procedimientos.

La metodología aquí aplicada se basa en el análisis estático no lineal y en el método del espectro de capacidad y ha permitido, de una manera fiable y claro está, en el campo probabilista, evaluar el desempeño de una estructura mediante los métodos descritos anteriormente.


21

5. COMPARACIÓN DEL NIVEL DE SOLICITACIÓN SÍSMICA QUE INDICA EL PUNTO DE DESEMPEÑO CON EL SUPUESTO AL DISEÑAR DE ACUERDO CON LA NSR-10 QUE PERMITE EL USO DE R, FACTOR DE DISIPACIÓN SÍSMICA PARA EL CÁLCULO DE LAS FUERZAS SÍSMICAS.

Definición

En los últimos años se han realizado importantes cambios o ajustes en las técnicas de diseño y evaluación de estructuras de edificios. Nuevas teorías en relación con el nivel de riesgo sísmico aceptable se han promovido al introducir en las normas el que ha sido llamado el diseño por comportamiento. En este enfoque se define un objetivo de seguridad o desempeño de la estructura, que va desde un estado previo al colapso con un daño severo hasta un estado sin daño estructural y no-estructural y un funcionamiento sin interrupciones del edificio, según varios niveles de solicitación sísmica. Dichas solicitaciones corresponden a sismos de diferente intensidad que varían en su período de retorno, y por lo tanto implican una estimación de la amenaza sísmica para varios sucesos con diferente probabilidad. De esta forma se diseña o se evalúa una estructura no sólo para una demanda sísmica, como ha sido tradicional, sino a varias demandas con probabilidades de excedencia diferentes.

Calculo del factor de disipación de energía.

R=Sa1/Sa2

Sa1= Punto desempeño espectro de diseño

Sa2= Punto desempeño espectro reducido

R obtenido = 1,035

R NSR -10= 7 (Tabla A.3.3)

Los valores obtenidos para el coeficiente de disipación de energía, R, son en cualquiera de los casos, mucho menores que los presentados en las NSR-10.

El coeficiente de disipación de energía depende de la demanda sísmica y no se entiende por qué en las NSR-10 fijan un valor que depende básicamente del sistema estructural usado.

Por otro lado y según estudios, para periodos largos (a), se supone que si la estructura se comportase elásticamente bajo la acción de una fuerza F, sufriría un desplazamiento ΔE, mientras que si su comportamiento fuera plástico idealizado, se alcanzaría el mismo desplazamiento para una fuerza reducida F/R. Para periodos cortos, la propuesta se basa en la igualdad de energías, suponiendo que la energía del sistema elástico es igual a la energía del sistema plástico

De las anteriores expresiones parten las que definen el R en varios de los códigos de diseño, incluyendo las NSR-10:

Periodos largos, T > 0.5 s R = μ

Periodos cortos, T ≤ 0.5 s R = 2μ −1


22

Los valores obtenidos para el coeficiente de disipación de energía varían notablemente, dependiendo del método de cálculo usado.

Literaturas existentes en la red sugieren mediante estudios y artículos que la NSR–10 incluya una aclaración en cuanto a la aplicación del R, pues en la práctica del diseño estructural se siguen al pie de la letra las indicaciones del código y se diseñan todos los elementos, incluyendo los de soporte, para las fuerzas reducidas, las cuales están muy por debajo de las que harán plastificar la estructura y por lo tanto se generarán fallas frágiles que pueden ocasionar el colapso.

Los artículos encontrados en la red también sugieren el uso de programas de uso común en ingeniería, los cuales traen herramientas que permiten evaluar en forma relativamente fácil la curva de capacidad de la estructura y su desempeño ante determinada demanda sísmica. Se sugiere que se usen estas herramientas para evaluar las solicitaciones que liberan los pórticos, y que se usen esos resultados no solo para evaluar el desempeño, sino para el diseño de los elementos de soporte y estructuras ligadas al pórtico.

Adicionalmente a los requisitos de resistencia que involucran el coeficiente de disipación de energía, los códigos de diseño, incluyendo la NSR-10, han adoptado criterios de rigidez cada vez más exigentes; tales requisitos impiden definir de antemano el nivel de carga sísmica reducida donde el comportamiento inelástico se presenta y por lo tanto los valores propuestos para el R, predefinidos en los códigos, entrarán en contradicción con los niveles de ductilidad que se pueden alcanzar con las estructuras diseñadas para cumplir esos requisitos de rigidez.

Por lo tanto y siendo consecuentes con lo anteriormente expuesto, los valores de R de la NSR – 10 como los calculados a partir de lo arrojado por el programa SAP 2000 tienen una diferencia abismal que no encuentra razón de ser, sabiendo de antemano que el tema ha sido ampliamente documentado y se han llegado a datos satisfactorios que no dan lugar a duda.


23

6. CURVA DE FRAGILIDAD HAZUS 99

-2,000E-01

0,000E+00

2,000E-01

4,000E-01

6,000E-01

8,000E-01

1,000E+00

1,200E+00

0 5 10 15 20 25 30 35 40

Curvas de Fragilidad

LEVE

MODERADO

SEVERO

COLAPSO


24

Mediante la obtención de las curvas de fragilidad se puede obtener el desempeño de la estructura.

Para obtener las curvas de fragilidad es necesario definir una medidas objetivas que indiquen, cuando el daño de una estructura pasa de un estado a otro; esas medidas se conocen como niveles de daño o estados limites de daño, los cuales son leve, moderado, severo y colapso

La fragilidad de un edificio está relacionada con su vulnerabilidad sísmica y se cuantifica mediante curvas de fragilidad.

En función del punto de desempeño. HAZUS 99 obtiene las curvas de fragilidad para encontrar el nivel de daño estructural y no estructural.

7. GRADO DE DAÑO QUE EXPERIMENTARÁ EL EDIFICIO


25

Según los datos obtenido anteriormente, y siguiendo los paso para obtener la información correspondiente al grado de daño vemos que el edificio se encuentra ante al probabilidad del 2% de que ocurra colapso, 18% de que ocurra daño severo, 68% de que ocurra daño moderado y un 90% de que ocurran daños leves, algo que puede sorprender ya que hablamos de un edificio cuyas características pueden reflejar un mal comportamiento ante eventos sísmicos, pues los edificios más bajos tienen un mayor índice de daño, lo cual puede ser debido a que los edificios de poca altura poseen mayor rigidez y menos masa que los edificios altos, por lo cual, el período es más bajo y las aceleraciones espectrales aumentan. Además los desplazamientos espectrales decrecen al disminuir el número de niveles de los edificios. También se puede concluir cómo el amortiguamiento efectivo aumenta a medida que disminuye el número de pisos, lo que puede estar relacionado con la disipación de energía de los edificios. Para una misma acción sísmica definida por medio de un espectro de demanda, las estructuras bajas se acercan más a los períodos cortos que los edificios altos.


26

CONCLUCIONES

La metodología aquí aplicada se basa en el análisis estático no lineal y en el

método del espectro de capacidad y ha permitido, de una manera fiable en el

campo probabilista, evaluar la vulnerabilidad sísmica mediante curvas de fragilidad

y matrices de probabilidad de daño.

Literaturas existentes en la red sugieren mediante estudios y artículos que la NSR–10 incluya una aclaración en cuanto a la aplicación del R, pues en la práctica del diseño estructural se siguen al pie de la letra las indicaciones del código y se diseñan todos los elementos, incluyendo los de soporte, para las fuerzas reducidas, las cuales están muy por debajo de las que harán plastificar la estructura y por lo tanto se generarán fallas frágiles que pueden ocasionar el colapso.

Los artículos encontrados en la red también sugieren el uso de programas de uso común en ingeniería, los cuales traen herramientas que permiten evaluar en forma relativamente fácil la curva de capacidad de la estructura y su desempeño ante determinada demanda sísmica. Se sugiere que se usen estas herramientas para evaluar las solicitaciones que liberan los pórticos, y que se usen esos resultados no solo para evaluar el desempeño, sino para el diseño de los elementos de soporte y estructuras ligadas al pórtico.

Adicionalmente a los requisitos de resistencia que involucran el coeficiente de disipación de energía, los códigos de diseño, incluyendo la NSR-10, han adoptado criterios de rigidez cada vez más exigentes; tales requisitos impiden definir de antemano el nivel de carga sísmica reducida donde el comportamiento inelástico se presenta y por lo tanto los valores propuestos para el R, predefinidos en los códigos, entrarán en contradicción con los niveles de ductilidad que se pueden alcanzar con las estructuras diseñadas para cumplir esos requisitos de rigidez.

Por lo tanto y siendo consecuentes con lo anteriormente expuesto, los valores de R de la NSR – 10 como los calculados a partir de lo arrojado por el programa SAP 2000 tienen una diferencia abismal que no encuentra razón de ser, sabiendo de antemano que el tema ha sido ampliamente documentado y se han llegado a datos satisfactorios que no dan lugar a duda.

Los métodos RISK UE Y ATC – 40 que para nuestro caso muestran valores algo semejantes entre sí, son el resultado de corregir mediante iteraciones valores que en primera medida no se acomodaban a los criterios de construcción grafica implementados y obtenidos a partir de formulas y que a su vez nacen de la información obtenida durante la modelación en el programa SAP 2000, y que a su


27

vez no son datos similares a los obtenidos en el programa pudiendo este fenómeno obedecer a procedimientos de modelación de la estructura no regulares desde el punto de vista de la teoría estructural, pero que a manera de instrucción a la hora de comprender la esencia de los métodos puede ser un instructivo a la hora de proporcionarnos una comprensión de los tres métodos utilizados, máxime no haber obtenido unos resultados equiparables entre si para el propósito de la vida profesional pero brindándonos el entrenamiento necesario para comprender la finalidad de los procedimientos.

El comportamiento calculado puede no ser el esperado por la Norma ya que como se vio por ejemplo en el factor de disipación sísmica R, el cual independientemente de los datos arrojados por el programa siempre es considerado mucho mayor por la norma que por los resultados obtenidos mediante los procedimientos utilizados para el desarrollo de este trabajo.

Puede ser posible el cumplir derivas aun entrando en el campo plástico siempre que no se halla experimentado un fenómeno elástico bastante severo y prolongado que halla inducido en la estructura deformaciones remanentes que no son susceptibles de recuperar.


28

BIBLIOGRAFIA

http://www.scielo.org.co/pdf/iei/v28n1/v28n1a05.pdf http://www.santaeufemiadelarroyo.es/extras/archivos/Caracter%EDsticas%20s%EDsmicas%20de%20las%20construcciones%20de%20tierra%20en%20el%20Per%FA.pdf http://www.manizales.unal.edu.co/gestion_riesgos/refuerzo3.php http://www.scielo.org.co/pdf/iei/v28n1/v28n1a05.pdf MANUAL DE ANALISIS NO LINEAL ESTATICO “PUSHOVER” USANDO SAP 2000 V14”. MORRISON INGENIEROS. ATC – 40 FEMA – 356 TESIS DOCTORAL PROFESOR CARLOS BERMUDEZ

http://www.scielo.org.co/pdf/iei/v28n1/v28n1a05.pdf

http://www.santaeufemiadelarroyo.es/extras/archivos/Caracter%EDsticas%20s%EDsmicas%20de%20las%20construcciones%20de%20tierra%20en%20el%20Per%FA.pdf



http://www.manizales.unal.edu.co/gestion_riesgos/refuerzo3.php

http://www.scielo.org.co/pdf/iei/v28n1/v28n1a05.pdf

Trabajo Final Analisis No Lineal Estatico Pushover - Copia

Documents

Transcript of Trabajo Final Analisis No Lineal Estatico Pushover - Copia