DISENO SISMICO COMPLETO -- EDIFICIO DE 2 PISOS

[ ]

mostrar "======================================================"
mostrar "  DISENO SISMICO COMPLETO -- EDIFICIO DE 2 PISOS"
mostrar "  Integra 8 piezas construidas hoy en una sola tuberia"
mostrar "======================================================"

E_c = 25000000.0
A_c = 0.16
I_c = 0.00213333
L_col = 3.2
L_viga = 5.0

mostrar ""
mostrar "--- PASO 1: Ensamblaje de rigidez ---"
crear_matriz_dinamica(18, 18, K)
rigidez_global_elem(E_c, A_c, I_c, L_col, 1.5707963, k_c1)
ensamblar_elemento_portico(K, k_c1, 1, 3)
rigidez_global_elem(E_c, A_c, I_c, L_col, 1.5707963, k_c2)
ensamblar_elemento_portico(K, k_c2, 2, 4)
rigidez_global_elem(E_c, A_c, I_c, L_viga, 0.0, k_v1)
ensamblar_elemento_portico(K, k_v1, 3, 4)
rigidez_global_elem(E_c, A_c, I_c, L_col, 1.5707963, k_c3)
ensamblar_elemento_portico(K, k_c3, 3, 5)
rigidez_global_elem(E_c, A_c, I_c, L_col, 1.5707963, k_c4)
ensamblar_elemento_portico(K, k_c4, 4, 6)
rigidez_global_elem(E_c, A_c, I_c, L_viga, 0.0, k_v2)
ensamblar_elemento_portico(K, k_v2, 5, 6)
mostrar "Portico de 2 pisos ensamblado (12 GDL libres)."

crear_matriz_dinamica(12, 12, K_red)
para i = 1 hasta 12
    para j = 1 hasta 12
        K_red[i,j] = K[i+6, j+6]
    finpara
finpara

mostrar ""
mostrar "--- PASO 2: Deriva bajo carga lateral ---"
F_carga = [80.0,0.0,0.0,0.0,0.0,0.0,50.0,0.0,0.0,0.0,0.0,0.0]
resolver_sistema_nxn(K_red, F_carga, desplaz)
desplaz_pisos = [desplaz[1], desplaz[7]]
graficar_barras(desplaz_pisos, "Desplazamiento lateral por piso (m)")

mostrar ""
mostrar "--- PASO 3: Periodos naturales (modelo de cortante) ---"
kxi_col = 12.0 * E_c * I_c / (L_col*L_col*L_col)
k_piso1 = 2.0 * kxi_col
k_piso2 = 2.0 * kxi_col
K_cort = [[k_piso1+k_piso2, 0.0-k_piso2],[0.0-k_piso2, k_piso2]]
M_cort = [[100.0,0.0],[0.0,100.0]]
analisis_modal(K_cort, M_cort, 100, frec_cort, modos_cort)
periodo1 = 6.28318530718 / frec_cort[1]
periodo2 = 6.28318530718 / frec_cort[2]
mostrar "T1 (periodo fundamental, s):"
mostrar periodo2
pgrafico = [periodo2, periodo1]
graficar_barras(pgrafico, "Periodos naturales del edificio (s)")

mostrar ""
mostrar "--- PASO 4: Espectro de diseno y respuesta sismica ---"
Sa1 = espectro_diseno_sa(periodo1, 1.0, 0.5)
Sa2 = espectro_diseno_sa(periodo2, 1.0, 0.5)
resp_modales = [[Sa1],[Sa2]]
srss_combinar(resp_modales, resp_total)
mostrar "Respuesta sismica combinada SRSS:"
mostrar resp_total[1]

pcurva = [0.05,0.1,0.2,0.3,0.4,0.5,0.6,0.8,1.0,1.5,2.0]
sa_curva = [0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0]
para i = 1 hasta 11
    sa_curva[i] = espectro_diseno_sa(pcurva[i], 1.0, 0.5)
finpara
graficar_barras(sa_curva, "Espectro de diseno sismico Sa(T)")

mostrar ""
mostrar "--- PASO 5: Chequeo P-Delta ---"
deriva_piso1 = desplaz[1] / L_col
suma_P = 3000.0
V_piso = 130.0
B2 = pdelta_b2(suma_P, deriva_piso1, V_piso, L_col)
mostrar "Factor de amplificacion B2:"
mostrar B2

mostrar ""
mostrar "--- PASO 6: Simetria en planta ---"
columnas_planta = [[-2.5,-2.5,0.0],[2.5,-2.5,0.0],[2.5,2.5,0.0],[-2.5,2.5,0.0]]
tiene_c4 = simetria_cn(columnas_planta, 4)
mostrar "Tiene simetria C4 en planta:"
mostrar tiene_c4

mostrar ""
mostrar "--- PASO 7: Clasificar seguridad con arbol de decision ---"
X_entren = [[0.5,0.005,1.1],[0.8,0.015,1.2],[1.5,0.025,1.5],[0.3,0.008,1.05],
            [2.0,0.03,1.6],[1.0,0.018,1.35],[2.5,0.04,2.0],[0.6,0.012,1.15]]
y_entren = [1,1,0,1,0,1,0,1]
crear_matriz_dinamica(16, 6, arbol_seguridad)
masc_entren = [1,1,1,1,1,1,1,1]
construir_nodo2(X_entren, y_entren, masc_entren, 8, 3, 2, 0, 4, 2, 1, arbol_seguridad)

T_fundamental = periodo2
punto_edificio = [T_fundamental, deriva_piso1, B2]
clasificacion = predecir_arbol(punto_edificio, arbol_seguridad, 1)

punto_malo = [2.2, 0.035, 1.7]
clasif_malo = predecir_arbol(punto_malo, arbol_seguridad, 1)

mostrar ""
mostrar "--- CLASIFICACION FINAL (1=seguro, 0=revisar diseno) ---"
mostrar "Nuestro edificio:"
mostrar clasificacion
mostrar "Caso deliberadamente malo (verificacion):"
mostrar clasif_malo

deriva_normalizada = deriva_piso1 / 0.02
B2_normalizado = (B2 - 1.0) / 0.4
comparacion = [deriva_normalizada, B2_normalizado, clasificacion]
graficar_barras(comparacion, "Resumen: deriva/limite, B2 normalizado, clasificacion")

Salida Jota: ====================================================== -Salida Jota: DISENO SISMICO COMPLETO -- EDIFICIO DE 2 PISOS -Salida Jota: Integra 8 piezas construidas hoy en una sola tuberia -Salida Jota: ====================================================== -Salida Jota: -Salida Jota: --- PASO 1: Ensamblaje de rigidez --- -Salida Jota: Portico de 2 pisos ensamblado (12 GDL libres). -Salida Jota: -Salida Jota: --- PASO 2: Deriva bajo carga lateral --- -Salida Jota: -Salida Jota: --- PASO 3: Periodos naturales (modelo de cortante) --- -Salida Jota: T1 (periodo fundamental, s): -Salida Jota: 0.514 -Salida Jota: -Salida Jota: --- PASO 4: Espectro de diseno y respuesta sismica --- -Salida Jota: Respuesta sismica combinada SRSS: -Salida Jota: 1.395 -Salida Jota: -Salida Jota: --- PASO 5: Chequeo P-Delta --- -Salida Jota: Factor de amplificacion B2: -Salida Jota: 1.014 -Salida Jota: -Salida Jota: --- PASO 6: Simetria en planta --- -Salida Jota: Tiene simetria C4 en planta: -Salida Jota: 1.000 -Salida Jota: -Salida Jota: --- PASO 7: Clasificar seguridad con arbol de decision --- -Salida Jota: -Salida Jota: --- CLASIFICACION FINAL (1=seguro, 0=revisar diseno) --- -Salida Jota: Nuestro edificio: -Salida Jota: 1.000 -Salida Jota: Caso deliberadamente malo (verificacion): -Salida Jota: 0.000 -

Shift + Enter ejecuta la celda actual

Ayuda de Jota